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Zur Phasenumwandlung zwischen Supra- und Normalleiter 
im kritischen Magnetfeld 


Von F. Beck 
(Mit 2 Abbildungen) 


Inhaltsiibersicht 


Fiir eine langsame, gleichférmige Verschiebung der Grenze zwischen Supra- 
und Normalleiter werden im Fall des unendlichen supraleitenden Halbraumes 
die Gleichungen der London-Laueschen Theorie gelést. Das Ergebnis recht- 
fertigt eine Herleitung der zum Zusammenbruch der Supraleitung fiihrenden 
kritischen Feldstärke aus dem Energiesatz der Elektrodynamik’). 


1. Einleitung 


Bis sich die Auffassung des „Zweiphasensystems“ Supra-Normalleiter durch- 
zusetzen begann, konnte man über den Übergang aus dem supra- in den normal- 
leitenden Zustand wenig aussagen. Die Zerstörung der Supraleitung durch ein 
starkes äußeres Magnetfeld war empirisch bekannt und man glaubte, daß dabei 
die das Feld abschirmenden Supraströme in Ohmsche Ströme verwandelt, deren 
Energie also irreversibel in Joulesche Wärme umgesetzt würde. Das war vor 
Entdeckung des Meißner-Effektes (1933) und diese führte dann auch zu jener 
Auffassung des Zweiphasensystems. Es setzten nun die bekannten Anwendungen 
der Thermodynamik reversibler Vorgänge ein (Rutgers, Casimir und Gorter), 
die unter Benutzung der Temperaturabhängigkeit des magnetischen Schwellen- 
wertes als empirische Tatsache den Supraleiter thermodynamisch als vollkommenes 
Diamagneticum (u = 0) behandeln. Diese Betrachtungen erfordern nicht die Lon- 
don-Lauesche elektrodynamische Theorie; sie gehen auf den Mechanismus der 
Supraleitung nicht ein und können daher auch nicht die bekannten Abweichungen 
von der Silsbeeschen Hypothese für dünne Supraleiter liefern. Es sei ausdrücklich 
betont, daß bereits das Hinschreiben der zwei Londonschen Ergänzungsglei- 
chungen ohne jede Zusatzannahme notwendig zur gesamten Thermodynamik 
des Zusammenbruches führt, wie dies M. v. Laue?) in voller Allgemeinheit zuerst 
gezeigt hat. Die v. Lauesche Betrachtung geht von den Gleichgewichtsbedin- 
gungen aus, die sich durch das Auftreten elektrodynamischer Spannungen an der 
Grenze zwischen Supra- und Normalleiter ergeben. Freilich sucht H. London 
bereits in seiner oben zitierten Arbeit den Schwellenwert, ebenfalls unter Zugrunde- 
legung der Londonschen Gleichungen aus dem Energiesatz der Elektrodynamik 


1) H. London, Proc. Roy. Soc. London 152, 650 (1935). 
2) M.v.Laue, ‘Ann. Physik 82, 71 (1938) (s. auch ,, Theorie der Supraleitung*, 2. Aufl. 
Berlin, Springer 1949). 
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abzuleiten (die Méglichkeit der Ableitung aus dem Spannungstensor ist nur ange- 
deutet), jedoch beschränkt sich seine Betrachtung auf den Fall einer ebenen, sich 
mit gleichférmiger Geschwindigkeit verschiebenden Grenzfläche. Sie bedarf ferner 
folgender zwei Annahmen: Erstens dürfen die bei der Bewegung der Grenze auf- 
tretenden elektrischen Feldstärken keine merkliche Joulesche Wärme hervor- 
bringen und zweitens muß der Poyntingsche Energiestromvektor an der Grenze 
stetig, d.h. es müssen die üblichen Grenzbedingungen für © und § auch an der 
bewegten Grenze erfüllbar sein. In der Londonschen Arbeit ist nur die erste An- 
nahme näher begründet). Die zweite Annahme zu rechtfertigen, gab den Anstoß, 
die Feldgleichungen für den vorliegenden Fall zu lösen. Da dies unseres Wissens 
bisher noch nicht geschehen ist, erscheint es gerechtfertigt, das Ergebnis kurz 
mitzuteilen. 


2. Lösung der Feldgleichungen für eine langsame, gleichförmige Verschiebung 
der Grenze zwischen Supra- und Normalleiter 


Wir beschränken uns auf den einfachsten Fall des unendlich ausgedehnten, 
isotropen supraleitenden Halbraumes und wenden die lineare Theorie an (4 = 
const. für 7’ = const.). Die Supra- und Normalleitung tragende Materie grenze 
in z = a an Vakuum an. 


Die zu erfüllenden Gleichungen sind: 


19 1 (a 
div = 0 (III), div © = 0 (ladungsfrei) (IV), { liter, 
° =o (V) ‘ 
zh o, im Supraleiter 
s Wz=V o = (a, im Normalleiter } (Va) 
0 im Vakuum 
aS 6 (vn), 
N Z*2Z9 Nur im 


c-Arot 3! = — (VII), { Supraleiter. 


V Z=a Wegen der unendlichen Ausdehnung 
des Halbraumes hängen alle in Be- 
Abb. 1. Die Grenze zwischen Supra- und tracht kommenden Größen nur von z 
Normalleiter genüge der Gleichung z, = vt und tab: (z,t); E= E(z, i); 


3) Das merkliche Auftreten Joulescher Wärme würde nicht der Reversibilitat des 
Phasenüberganges Supra-Normalleiter widersprechen, da die hierbei irreversibel umge- 
setzte Energie aus dem äußeren Felde, dessen Gesamtenergie in die thermodynamischen 
Rechnungen nicht eingeht, stammen könnte. Dies ist z. B., wie G. U. Schubert) ge- 
zeigt hat, beim Eintritt des Meißner-Effektes der Fall, wo die während der Feldver- 
drängung auftretenden Ohmschen Ströme die gleiche Größenordnung haben wie die ent- 
stehenden Supraströme. 


4) G. U. Schubert, Ann. Physik 5, 213 (1950). 
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Y= (zt); (st). Zunächst folgt dann aus den beiden Divergenz- 
gleichungen (III) und (IV) $,= 0 und €,=0. Wir versuchen nun, die Feld- 
gleichungen für die einzelnen Bereiche zu lösen. 


a) Supraleitender Bereich 


Es soll angenommen werden, daß die Feldverteilung im Supraleiter in bezug 
auf die bewegte Grenze unverändert bleibt. Dies trägt dem Umstande Rechnung, 
daß das Feld nur in einer dünnen Schutzschicht hinter der Grenze vorhanden ist 
und von dieser während ihrer Bewegung mitgeschleppt wird. Ferner führt es zu 
den von der Thermodynamik geforderten konstanten Werten von $ und § an 
der Grenze. Dann kann die Abhängigkeit von z und ¢ nur in der Kombination 
z— vt auftreten. Dies führen wir als neue Variable ein: 


C=2-vt. (1) 
Dann wird 


Wir wählen nun noch die Achsenrichtungen so, daß allein $, von Null verschieden 
ist, und setzen an (£ >0): 


=1' 0). (2) 
Dann folgt aus (I) durch Integration 
Ga-ivth (3) 


Die Konstanten k, und k, verschwinden, weil in größerem Abstand von der Grenze 


das Innere des Supraleiters feldfrei sein muß. Gemäß (V) ist damit auch der 
Ohmsche Strom bestimmt. Für den — folgt aus (VI) 


27 (0) dl + kz. (4) 


k, ist wieder im Hinblick auf das cae Innere zu bestimmen. SchlieBlich er- 
gibt (II), zusammen mit (VIII), und Integro-Differentialgleichung 


0) 
aus der durch Differentiation die Differentialgleichung 
[23 1 


hervorgeht. Die noch verbleibende Gl. (VII) liefert nichts Neues und ist mit 
dem gemachten Ansatz verträglich. 
Die für  -> co verschwindende Lösung von (5a) ist leicht anzugeben: 


7? = He (6) 


+4 1 (3 1 0,% 


mit 


yi-s B(1 5) V 


Da die Konstante k, in (4) gerade so zu bestimmen ist, daß (5) das unbestimmte 
Integral der homogenen Gl. (5a) ist, stellt (6) auch die Lösung von (5) dar. Dies 
läßt sich auch durch Einsetzen leicht bestätigen. 
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b) Normalleitender Bereich (a — vt < ¢ < 0) 


Auch hier wollen wir die alleinige Abhängigkeit von £ annehmen; wie sich 
nachher ergibt, ist dies zur Erfüllung der Grenzbedingungen notwendig: 


92 = 7" (~). (7) 
Aus (I) folgt, ebenso wie im Supraleiter, 
Ü 
und (II) führt mit (V) und (Va) zu der Differentialgleichung 
+ =0 
Als Differentialgleichung 1. Ordnung hat sie nur die eine unabhängige Lösung 
aft = 


mit 


Die Verbindung der beiden Lösungen (6) und (9) ergibt sich aus den Grenz- 
bedingungen. Im vorliegenden Fall, wo D = E und $ = @ ist, lassen sich diese 
auch für die bewegte Grenze leicht aus (I) und (II) herleiten. Dazu legen wir ein 
im Raum festes Rechteck der langen Seite 1 und der kurzen Seite d so, daß jene 
mit Ctang parallel wird und diese normal zur Grenzfläche verläuft. Integrieren 
wir (I) über die Fläche des Rechtecks, im Augenblick, wo die Grenze über dieses 
hinwegstreicht, so ergibt sich aus 


_ 
(3), Ao=d 


und 


die Gleichung 
— = lim dv —= (1 — Dra). 


d>0 
Entsprechend ergibt sich aus (II) 


- Da = Jim Gus d+ a} — G2), 


© 
da % endlich bleibt. Die zwei letzten Gleichungen sind nur miteinander verträglich, 


wenn 
Dia = Hr und = Cn 
ist. 

Nach (2), (6), (7) und (9) bedeutet dies 


H* — Hn = Ho. 


Wie ein Blick auf (3) und (8) zeigt, sind damit die Grenzbedingungen für € eben- 
falls erfüllt. 


ec) Vakuum (z < a) 
In diesem Bereich können wir nicht mehr die Abhängigkeit von z und ¢ in der 
Kombination ¢ = z— vt vorschreiben (dies würde auf 9, = const., €, = const. 
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führen, was mit den Randbedingungen in z = a nicht verträglich ist). Wir setzen 


daher an 
9.=n(2t); €,=e(,d); 
(I) und (II) führt auf die Differentialgleichungen 
_ 1a, 1 


cat’ ca 
mit den fiir den vorliegenden Fall brauchbaren Lésungen 


n= A, erhtz-et) 4 dhs ehlzitc), gg = — A, erhte-et) + A,ehtetet), 


Die Konstanten A,, A, und h ergeben sich aus den Randbedingungen in z = a. 
Die Grenzbedingungen fordern hier wiederum Stetigkeit von Stang und Ctang, 
so daß jene nach (7), (8) und (9) lauten 


Dies führt auf folgende Werte 
h v 


Damit ist die Lösung vollständig bestimmt. Für die einzelnen Feldgrößen er- 
geben sich in den drei Bereichen die folgenden Gleichungen. Wir geben neben der 
exakten Lösung auch die für kleine Geschwindigkeiten gültige Näherung unter 
Vernachlässigung von Gliedern zweiter Ordnung in v/c an. 


¢= 0 (supraleitend): 


1 0% 
= H! e- IS He ~(1+ 2 
1 
= — yi BE — 
c c (10) 
1 Osv 
1 o,v\ aq) 
1 Ogv 
= — HZ enßt — Ho ov 
c 
2 
n= y! +7 :G 1— 1 
a—vt<f€< 0 (normalleitend): 
= He rn: w He 
(11) 
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z <a (Vakuum): 
= + 83} 
104 =— Fern {s, — Sj} 


v onr® 


Diese Gleichungen zeigen folgendes: Im supraleitenden Bereich ist die Lösung 
ganz ähnlich der für ruhende Grenzflächen®). Die Felder klingen mit wachsendem £ 
sehr rasch exponentiell ab. Das Auftreten des Ohmschen Stromes bedingt lediglich 
eine geringe Änderung der Eindringtiefe. Im stationären Fall B-1, ist sie hier 
gegeben durch 


(12) 


(y By BO (1 


Für positives v, d.h. Wachsen des normalleitenden Bereichs, ist sie 
im umgekehrten Fall, für negatives v, größer. Allerdings ist das Gliea 2 $ I von 
der Größenordnung 10? v/c (a ~ 10, 6 ~ 108), also für kleine Geschwindigkeiten 
(v/ce ~ 10-10) kaum ins Gewicht fallend. 

Im normalleitenden Bereich stellt die Lösung kein homogenes Feld mehr dar, 
wie im stationären Fall, die Feldgrößemwachsen vielmehr bei positivem v mit der 
Entfernung von der Grenzfläche (negative ¢-Richtung) exponentiell an. Dieses 
Auftreten der Exponentialfunktion in (11) ist eine Folge der enegieverzehrenden 
Ohmschen Ströme. Denn die insgesamt erzeugte Joulesche Wärme wächst 
mit der Dicke der normalleitenden Schicht, und damit bei positivem v mit der 
Zeit an. Diese im Laufe der Zeit immer größer werdende in Wärme umgesetzte 
Energie muß, da das Feld an der Begrenzung des supraleitenden Bereichs seinen 
Wert dauernd beibehalten soll, als Energiestrom aus dem Vakuum nachgeliefert 
werden. Da einer Vergrößerung des Poyntingschen Vektors eine Vergrößerung 
der Feldstärken entspricht, müssen diese exponentiell anwachsen. Die Strecke, 
auf der das Feld nicht wesentlich von seinem konstanten Wert bei £ = 0 abweicht, 
kann jedoch mit abnehmender Geschwindigkeit beliebig groß gemacht werden, 
wie das lineare Auftreten von v in der e-Funktion zeigt. 

Ebenso ist es im Vakuum. Dort hat die Lösung natürlich nur bis in die Nähe 
der Quellen des Feldes Gültigkeit. Wenn diese nicht zu weit von der Grenze 
z =a entfernt sind, kann man wegen des Faktors (v/c)? in der Exponentialfunk- 
tion [Gl. (12)] das Feld in erster Näherung als räumlich konstant und für nicht 
zu große Zeiten, d.h. nicht zu große Dicke des normalleitenden Bereichs, auch 
als zeitlich konstant ansehen. 

Der Lösung des vorliegenden Problems ist der Einfachheit halber die skalare 
Theorie zugrunde gelegt worden. Es macht jedoch keine prinzipiellen Schwierig- 
keiten, die von M. v.Laue®) entwickelte, allen nichtkubischen Supraleitern 


5) 8. 2. B. das zitierte Buch von M. v. Laue, § 7c. 
6) M. v. Laue, Ann. Physik 8, 31 (1948). 


F.B 
gere 
lich 
Erg 
Fel 
Gre 
erfü 
tin 
pro 
im 
iibe 
den 
gesi 
klei 
The 
geä 
mit 
tris 
vol 
doc 
die 
8. 
sei 
sti 
no 
En 
un 
gre 
un 
un 
un 
mi 
su 
0, 
bz 
un 
Eı 
eil 


(12) 


F. Beck: Phasenumwandlung zwischen Supra- und Normalleiter im kritischen Magnetfeld 323 


gerecht werdende tensorielle Theorie heranzuziehen. Die Lösung wird dann ähn- 
lich modifiziert wie im stationären Fall’). 

Im Hinblick auf die Londonsche Ableitung des kritischen Feldes zeigt das 
Ergebnis, daß es in der Tat möglich ist, für die bewegte Grenzfläche Lösungen der 
Feldgleichungen für Supra- und Normalleiter zu finden, mit denen sich die üblichen 
Grenzbedingungen (Stetigkeit der tangentiellen Komponenten von € und $) 
erfüllen lassen. Damit ist die für jene Ableitung notwendige Stetigkeit des Poyn- 
tingschen Vektors gewährleistet. Ferner zeigen Gin. (10) und (11), daß die 
pro Zeit- und Volumeneinheit erzeugte Joulesche Wärme sowohl im supra- wie 
im normalleitenden Bereich der zweiten Potenz von v/c proportional ist. Gegen- 
über der an der Grenze geleisteten, ihrer Verschiebung proportionalen Arbeit kann 
demnach die Joulesche Wärme beliebig klein gemacht werden. Damit wird der 
gesamte Vorgang der Grenzverschiebung durch das kritische Feld für unendlich 
kleine Geschwindigkeiten reversibel, ganz entsprechend den Vorgängen, die die 
Thermodynamik betrachtet. 

Weiter widerlegt das Ergebnis die von H. London in seiner zitierten Arbeit 
geäußerte Ansicht, das Entstehen bzw. Verschwinden des Suprastromes un- 
mittelbar an der bewegten Grenze ginge spontan, ohne jede Hilfe eines elek- 
trischen Feldes vor sich (I. c., 8. 655). Dann nämlich wäre die Ableitung des 
vorstehenden Ergebnisses allein aus der Elektrodynamik nicht möglich, bei der 
doch die zeitliche Änderung des Suprastromes und elektrische Feldstärke durch 
die Grundgleichung (VI) verknüpft sind. 


3. Ableitung der kritischen Feldstärke aus dem Energiesatz der Elektrodynamik 


Da die eingangs zitierte Arbeit H. Londons wenig bekannt geworden zu 
sein scheint, sei es hier gestattet, die bekannte Beziehung zwischen kritischer Feld- 
stärke und Differenz der freien Energien von normal- und supraleitender Phase 
nochmals aus dem Energiesatz herzuleiten. 

H. London geht aus von der für reversible, isotherme Prozesse gültigen 
Energiebeziehung 

F: Gesamte freie Energie 
öF—-dA=0 4: An dem System geleistete Arbeit (13) 


und wendet sie auf den in Abb. 2 dargestellten Bereich an. In einem beliebig abge- 
grenzten, festen Volumen V seien supra- 
und normalleitende Phase anwesend, 
und ihre ebene Grenzfläche b bewege sich 
unter dem Einfluß eines MagnetfeldesH, 
mit der Geschwindigkeit v in den 
supraleitenden Bereich hinein. o, bzw. 
o, seien die Oberflächen der normal- 
bzw. supraleitenden Volumenteile V,, 
und V,. 


Die Änderung der gesamten freien 


Energie in der Zeit dt setzt sich aus App,2, Zur Ableitung der kritischen Feld- 
einem Feld- und einem Phasenanteil stärke aus dem Energiesatz 


1. ce, § 7i. 
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zusammen: 


OF =(f,—f,)b-vdt+d [ Wav 
( nt ) 
f„ ist die aufs Volumen bezogene freie Energie der Normalphase, /, entsprechend 
für die supraleitende Phase. W ist die freie Feldenergie pro Volumeneinheit. Die 
an dem System geleistete Arbeit ist 


6A = — Ot f S,, do ©: Poyntingscher Vektor, 
(9) 
so daß sich ergibt 


f Wdv+6t fS,do=0. (14) 
(Vn+Vs) (9) 

Die Summe der beiden letzten Terme läßt sich nun berechnen, ohne die Feld- 
größen im einzelnen zu kennen unter Voraussetzung der Gültigkeit der Maxwell- 
schen Gleichungen in Supra- und Normalleiter, der Londonschen Gleichungen 
nur im Supraleiter und der Stetigkeit von €, und §, an der Grenzfläche. London 
geht aus vom Energiesatz, der für beide Phasen gleichermaßen lautet ®) 


+ dvS+0®=0. (15) 
Dabei gilt für W und © 
im Supraleiter: 


G 
W = + + J (i, dg)®) 


im Normalleiter: 

W= 9} 

S=clG, 9]. 


Wie aus Abschnitt 2) hervorgeht und wie man auch direkt aus (I) folgern kann, 


ist E mv v/c, alle in € quadratischen Glieder können daher für kleine Geschwindig- - 


keiten vernachlässigt werden. Dann ergibt Integration von (15) über das in Abb. 2 
angedeutete Volumen V 


av + --| [ + (18) 
(Vn) (Vs) (On) (os) (0) 

da ©, an der Grenzfläche b stetig ist. Beim Übergang zur totalen zeitlichen 

Änderung der Feldenergie ist zu bedenken, daß die Volumina, über die auf der 

linken Seite zu integrieren ist, selbst zeitabhängig sind. Es ergibt sich 


d ow ov ow 


(Vn) (Vn) (Vn) 
8) 8. z. B. M. v. Laue, l.c., § 5. 
*) Wir formulieren die Supraleitungsenergiedichte W' gleich allgemein (1. c., § 20). 
Fiir den linear skalaren Fall wird 


ia dg) = 443". 
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W. ist die Energiedichte an der sich zeitlich ändernden Begrenzung des Volumens. 
Entsprechend wird 


d aw 
[var = [er - an (20) 
(Vs) (Vs) 
Setzt man (19) und (20) in (18) ein, so erhält man wegen 9, = 9) (s. Abschnitt 2) 
Go 
wav+4 fwav=4 WdV=— | G,do— [ (i,dg)-bv 

at dt J (i’, dg) 
oder 


Go 
6 f war=-ä [ S,do—bvot f (it,dg). 
(Vn+Vs) (9) 0 


Zusammen mit (14) ergibt sich 
J = Wo. (21) 


Dies ist die bekannte!) Beziehung fiir die Grenzenergiedichte, die zum Zusammen- 
bruch der Supraleitung führt. Für die linear skalare Theorie geht (21) über in 


(21a) 
Mit Hilfe der in 2) gegebenen Lösung läßt sich dieses für den dortigen Spezial- 
fall des „‚dicken‘‘ Supraleiters explizit herleiten. Dazu nehmen wir als Volumen V 


einen Zylinder in z-Richtung, dessen eine Grundfläche von der Größe der Flächen- 
einheit im feldfreien Teil des Supraleiters, die andere im Normalleiter in unmittel- 
barer Nähe der Grenze liegt, so daß ¢ <; Zn ist. Dann kann man die Exponen- 
tialfunktion in (11) entwickeln und erhält aus (11) bis auf Größen zweiter Ord- 
nung in v/c 
= —c EH, = Hv. 
(Zylinder) 
Die Änderung der Feldenergie wird, da sich das Feld unverändert mit der Grenze 
verschiebt, 
6 WdV=Wvdt =} Hv dt. 
(Zylinder) 
Aus (14) folgt dann 
4. 


Dies ist die für dicke Supraleiter gültige Beziehung für die kritische Feldstärke. 
In der Schutzschicht des Supraleiters wird der Energiestromvektor nach (10) 


1 Ogv 
&= Hoye Ctr 
Dies gibt der Tatsache Ausdruck, daß die gesamte in der Schutzschicht enthaltene 
Energie eine feste Fläche als Poyntingscher Energiestrom durchsetzt. 


10) In dieser, auch für die nichtlineare Theorie gültigen Formulierung erstmals bei 
M. v. Laue, Ann. Physik 5, 197 (1950). (Dort bedeutet f die freie Energie pro Mol.) 


_| 
7) 
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Für v < 0, d. h. Vergrößerung des supraleitenden Bereichs, ergibt sich aus der 
Lösung (11) auch im Normalleiter ein exponentiell abklingendes Feld. Läßt man 
die Begrenzung der Materie nach —oo rücken, dehnt man also den Normal- 
leiter nach negativen z-Werten hin unendlich aus, so bedeutet das abklingende 
Feld physikalisch, daß die bei der Umwandlung Supra-Normalleiter gewonnene 
Energie in diesem vollständig in Joulesche Wärme umgesetzt wird. Die Bewegung 
ist dann ,,selbsterhaltend“, d. h. sie findet, einmal eingeleitet, statt, ohne daß das 
Feld von außen beeinflußt zu werden braucht. 

Sämtliche hier angestellten elektrodynamischen Betrachtungen machen keine 
Annahme darüber, welcher Mechanismus zur Bewegung der Grenze zwischen 
Supra- und Normalleiter führt. Sie gelten in ihrer Näherung für kleine Geschwin- 
digkeiten auch dann, wenn der Supraleiter an einen ihm chemisch verschiedenen 
Normal- oder Nichtleiter angrenzt, das Feld H® an der Grenze einen beliebigen 
Wert <H, hat, und jener sich als Ganzes mit der Geschwindigkeit v verschiebt. 
Denn die Maxwellschen Gleichungen (I) bis (IV) gelten unverändert auch für 
bewegte Materie, während Ohmsches Gesetz und Londonsche Gleichungen von 
der angegebenen, im Ruhsystem gültigen Form, erst in zweiter Ordnung in v/c 
abweichen. Dann liefert die Energiebilanz die Arbeit, die man leisten muß, um 
den Supraleiter gegen die in der Grenzschicht wirkenden Spannungen zu ver- 
schieben. 

Zum Abschluß sei darauf hingewiesen, daß man gewöhnlich die vom Felde an 
der Materie geleistete Arbeit aus der Änderung des elektrodynamischen Potentials 
erhält!) und dessen Beziehung zur freien Feldenergie von den Quellen des Feldes 
(nur permanente Magnete, nur Ohmsche Ströme oder beides zusammen) abhängt. 
Daß wir hier trotzdem die Arbeit aus der Änderung der Feldenergie berechnen 
konnten, liegt daran, daß wir ein abgeschlossenes Volumen ohne Einschluß von 
Quellen, aber unter Berücksichtigung der ein- bzw. ausströmenden Energie be- 
trachtet haben. Die Energiebilanz muß für einen solchen Bereich natürlich 
stimmen. Für den ganzen, Feld und Quellen einschließenden Raum ist jedoch 
allein die Änderung des elektrodynamischen Potentials maßgebend. 


Herrn Professor Dr. Max v. Laue möchte ich für klärende Diskussionen meinen 
herzlichsten Dank aussprechen. 


1) Bei Anwesenheit von Supraleitern im Feld angegeben von M. v. Laue, Z. Physik 
125, 571 (1949). 


Berlin-Dahlem, Kaiser-Wilhelm-Institut für physikalische Chemie und 
Elektrochemie. 


(Bei der Redaktion eingegangen am 2. November 1951). 
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Die Bestimmung sehr hoher Verteilungstemperaturen') 
Von Ernst Palm 


(Mit 4 Abbildungen) 


Inhaltsübersicht 


Unter der Voraussetzung, daß zwei Strahlern zwei voneinander verschiedene 
Verteilungstemperaturen zugeordnet werden können, lassen sich ohne weitere An- 
nahmen beide Temperaturen ledighich durch Messung der Intensitätsverhältnisse 
der Spektren beider Strahler bestimmen. Die Unsicherheit der Bestimmung 
hängt nicht nur von der Genauigkeit der Messung, sondern auch von den beiden 
Temperaturen selbst ab„deren relative Differenz groß und die auch selbst möglichst 
hoch sein müssen. Die Methode erlangt ihren besonderen Wert deshalb erst bei 
der Bestimmung extrem hoher Verteilungstemperaturen, wie sie gelegentlich 
in Gasentladungen oder bei astrophysikalischen Objekten vorkommen. 


$ 1. Die Voraussetzungen für die Anwendung des Planckschen Gesetzes. 
Abkürzungen 


Wenn sich die relative Intensitätsverteilung irgendeines kontinuierlichen 
Spektrums in der Form des Planckschen Strahlungsgesetzes 


3 


(1) 


E,, 7 = const - 


darstellen läßt, sei der „„Parameter‘‘ 7’ entsprechend einem Vorschlag von Reeb 
und Richter?) als ,,Verteilungstemperatur der untersuchten Strahlung defi- 
niert. Allein in diesem Sinne®) wird dieser Begriff im folgenden benutzt werden. 
Wir setzen dabei voraus, daß es überhaupt Strahler mit einem derartigen Spek- 
trum gibt. Streng genommen trifft das nur für Hohlraumstrahler zu, und auch 
da sind bestimmte Bedingungen einzuhalten, damit die den Hohlraum verlassende 
Strahlung eine ungestörte Plancksche Strahlung darstellt. Ist nämlich bei- 
spielsweise der Hohlraum mit einem Medium hoher Dichte gefüllt, dann treffen 
die für die Ableitung des Planckschen Gesetzes geltenden Annahmen nicht mehr 
zu, und es können bei hohen Temperaturen Abweichungen auftreten, die ihren 


1) Mitgeteilt auf der Tagung des Verbandes Deutscher Physikalischer Gesellschaften 
in Bad Nauheim im Oktober 1950. 

2) O. Reeb u. M. Richter, Lichttechnik 2, 186 (1950). 

3) Gelegentlich wird statt dessen der Ausdruck „Farbtemperatur‘“ benutzt; jedoch 
sollte der Begriff Farbtemperatur nur in Verbindung mit dem sinnlichen Eindruck der 
Strahlung verwendet werden, wobei dann freilich unter Umständen auch einem Selektiv- 
strahler eine Farbtemperatur zugeordnet werden kann. 
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Ursprung in der angeregten Linienstrahlung der Hohlraumfüllung haben. Bei 
anderen als Hohlraumstrahlern ist das Plancksche Gesetz niemals streng erfüllt. 
Bei der Strahlung fester glühender Körper wäre Voraussetzung, daß der Absorp- 
tionsgrad *) des Körpers im gesamten Spektralgebiet konstant ist oder — wenig- 
stens im Bereich der Wienschen Näherung — dem Gesetz 


Oy = 
folgt. 

Die Anwendbarkeit des Strahlungsgesetzes wiirde durch diese Umstinde 
stark eingeschränkt werden, wenn nicht eine Reihe von kontinuierlichen Strah- 
lungen durch die Plancksche Formel überraschend gut angenähert werden könn- 
ten, wenigstens in einem größeren Spektralgebiet, das in der Regel auch den sicht- 
baren Bereich enthält. Das trifft nicht nur für eine große Zahl glühender fester 
Körper zu (den positiven Krater des Kohlelichtbogens, glühenden Wolframdraht 
usw.), sondern oft auch für die kontinuierliche Strahlung glühender Gase (Edelgas- 
hochdruckentladungen, astrophysikalische Objekte). Bei letzteren ist die große 
Schichtdicke von Sternatmosphären von besonderer Bedeutung, da in diesem 
Falle auffallende Strahlung schließlich nahezu völlig absorbiert wird®) und damit 
die Grundvoraussetzung der Hohlraumstrahlung erfüllt fst. 

Nachdem es also als erwiesen gelten kann, daß kontinuierliche Strahlung sehr 
häufig durch eine Plancksche Verteilung nach (1) so gut angenähert werden kann, 
daß man ihr wenigstens in beschränkten Spektralgebieten ®) eine Verteilungs- 
temperatur zuordnen kann, besteht das Problem der Messung und Berechnung 
dieser Temperatur, das auf mannigfache Weise gelöst werden kann (siehe § 2). 
Hier soll ein Weg angegeben werden, dessen besondere Eigenart darin besteht, 
daß, obwohl ein bekannter Vergleichsstrahler nicht benutzt wird, dennoch auf 
Absolutmessungen von Strahlungsenergien völlig verzichtet wird. Die besonderen 
Voraussetzungen, unter denen dieser Weg zum Erfolg führen kann, werden ein- 
gehend diskutiert werden. 


An Abkürzungen führen wir ein: 


die u-bezogene Wellenzahl 
_ 


A 
= Anzahl der Wellen je 1 Länge 


das 7-Maß 
w— 1 entspricht der Wellenlänge A = 1 u (Ultrarot) 
t = 1 entspricht der Temperatur T = 14380° K. 


4) Statt „Absorptionsvermögen‘‘, siehe O. Reeb, Lichttechnik 2, 206 (1950). 

5) A. Unsöld, Physik der Sternatmosphären. Berlin 1938; Pannekoek, Astrophys. 
J. 84, 481 (1936). 

6) Zu diesem Gebiet möge außer dem Sichtbaren das mit normalem Aufwand zu- 
gängliche Ultrarot und Ultraviolett gehören, also } etwa von 220 bis 3000 nm. 

7) Zugrunde liegt der auf der 9. Generalkonferenz für Maß und Gewicht in Paris im 


Oktober 1948 für c, = = festgelegte Wert: c, = 1,438 cm grad. 
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Damit wird: 
ETF 
w3 
Ey. = const - (2) 


Im folgenden beziehen sich die Indices 1 und 2 auf zwei zu vergleichende 
Strahlungen mit den Temperaturen 7, und 7. Wir werden noch abkürzen: 


=3(U +7) 


t_= Im — 7). 


§ 2. Zusammenstellung bekannter Methoden 


Einige der üblichen Methoden zur Bestimmung von Farb- bzw. Verteilungs- 
temperaturen seien hier kurz erwähnt. 


1. Die Methode der ausgewählten Ordinaten (,,selected ordinate method‘) °) °) 


Es wird auf dem Weg über die drei Farbenindices die dem visuellen Eindruck 
entsprechende ähnlichste Farbtemperatur nach einem von Bowditch und Null 
angegebenen Rechenverfahren gefunden, wobei man von der gemessenen absoluten 
Intensitätsverteilung ausgeht. 


2. Die Farbtemperaturgerade !°) 


Trägt man In => über w auf, so erhält man im Gültigkeitsbereich der Wien- 
schen Näherung des Planckschen Gesetzes 
= const - w* - (3) 


eine Gerade, deren Steigung ein Maß für die Verteilungstemperatur ist. Auch hier 
muß die absolute Intensitätsverteilung zuvor gemessen werden. 


8. Die Methode der Isochromaten !!) 


Auch diese Methode wird gewöhnlich im Bereich der Wienschen Näherung 
benutzt. Man vergleicht den unbekannten Strahler mit einer Standard-Licht- 
quelle, die auf eine Reihe verschiedener Temperaturen eingestellt wird. Es muß 
also bei dieser Lichtquelle die Abhängigkeit der Temperatur (7',) von der Belastung 
bekannt sein. Man mißt und trägt für einige w-Werte als Parameter über r, auf: 

In = — w (t,—T,) + const. (4) 
Diese Geradenschar hat einen gemeinsamen Schnittpunkt, dessen Abszisse die ge- 
suchte Verteilungstemperatur 7’, liefert. 


4. Die Methode der Isothermen 


Sie entspricht der Isochromaten-Methode; nur bleibt die Vergleichslichtquelle 
auf einer bestimmten und bekannten Temperatur 7,. Unabhängig Veränderliche 


8) F.T. Bowditch u. M. R. Null, J. Opt. Soc. Amer. 28, 500 (1938). 
®) D. B. Judd, J. Opt. Soc. Amer. 26, 421 (1936). 

10) H. G. MacPherson, J. Opt. Soc. Amer. 30, 189 (1940). 

11) K. Langbein, Diss. Basel 1918. 
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ist nicht die Temperatur, sondern die Wellenzahl. Auch dann wird Gl. (4) durch 
eine Gerade wiedergegeben, deren Steigung 2 t_ ist: 


Liefert die Messung keine Gerade, dann wird durch diese Gleichung die ,,Grada- 
tions- oder Gradiententemperatur“ 12) als eine Funktion von w definiert. Im andern 
Falle kann der Differentialquotient durch den Differenzenquotient ersetzt werden 
und ist dann u. a. auch die Grundlage der Farbtemperaturmessung mit dem 
Pulfrich-Photometer. 

5. Auf die Möglichkeit, ohne bekannten Vergleichsstrahler nur aus Relativ- 
messungen die Verteilungstemperatur bestimmen zu können, ist zum erstenmal 
von E. Warburg") hingewiesen worden. Den damals — 1917 — gegebenen ex- 
perimentellen Möglichkeiten entsprechend, wurde das Strahlungsintensitäts- 
verhältnis für die Gesamtstrahlung thermoelektrisch und für eine geeignete Wellen- 
länge im sichtbaren Gebiet photometrisch gemessen. Bezeichnen wir ersteres mit 
x und letzteres mit y und seien die Temperaturen der beiden miteinander ver- 


glichenen schwarzen Strahler 7, und 7',, denen die 2 -Maße 7, und T, ent- 


sprechen, dann wird: 


1 In 
y 
Vx-1 


Obwohl dieses Verfahren zur Erzielung einer größeren Genauigkeit sich noch 
modifizieren läßt, hat es sich nicht durchsetzen können, da die Anforderungen an 
die Meßgenauigkeit weder durch die Meßmethoden selbst noch durch Gewähr- 
leistung einer hinreichenden Konstanz des Strahlers erfüllt werden konnten. 
Erst eine Erweiterung der spektralen Messungen in das sehr langwellige Gebiet 
erschließt bei ausschließlicher Anwendung objektiver Meßmethoden einen Weg, 
die Verteilungstemperaturen beider Strahler mit befriedigender Genauigkeit zu 
bestimmen. Wesentlich ist, daß man in einem Spektralgebiet mißt, in dem die 
Abweichungen von der Wienschen Näherung (3) möglichst groß sind, da es theore- 
tisch unmöglich wäre, Temperaturen aus Relativmessungen zu bestimmen, wenn 
die untersuchten Strahlungen streng dem Wienschen Gesetz folgen würden. 


$ 3. Das Intensitätsverhältnis grauer Strahlungen 


Das im folgenden beschriebene Verfahren beruht zufolge dieser Überlegungen 
darauf, daß man das Intensitätsverhältnis zweier als annähernd grau vorausge- 
setzter Kontinua in einem möglichst ausgedehnten Spektralgebiet mißt. Der 
theoretische Verlauf dieses Verhältnisses, dessen Logarithmus wir mit n bezeichnen 
wollen, ist für die beiden Strahlungen in ähnlicher Weise charakteristisch wie 
die Plancksche Verteilung für die einzelne Strahlung. 


1) H.Kienle, Nullpunkt und Skala der Temperaturen der Sterne, Nachr. Ges. 
Wiss. Göttingen, Fachgr. II, Neue Folge 8, Nr. 17, 201. 
13) E. Warburg, Verh. Dtsch. Physik. Ges. 19, 3 (1917). 
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Für zwei Strahlungen mit den sees T, und r, und den Emissionsgraden ¢, 


und &, ist: 


ne(vo) = In + In 


wt. 


= In? + (w). (5) 
& 
Nach einer elementaren etic. isles wir für 7: 
Sin *(r;, 
g=t.0+ (6) 
Sin * (T+ — t_) 
In Abb. 1 ist 7 (w) für einige t-Wertepaare aufgetragen. 
Für kleine w ist, wenn man - in eine SEGEN nach w entwickelt: 


und also 2 


dn 
= 
Für große w ist asymptotisch 
und also 7 
d 
( =) =2:.. (8) |_ 00 
dw oo 7 
L — 77680 
§ 4. Die Berechnung 
der Verteilungstemperaturen 
Die Gln. (7) und (8) zeigen De, 
uns, wie man T_ aus dem gemes- 
senen Verlauf von 7(w) bestimmen 
kann, wobei man grundsätzlich 7 Sichtbares Geb. 
sowohl das kurzwellige als auch 0 7 2 
das langwellige Ende des Spek- w 


trums benutzen kann. In der Abb. 1. Das spektrale Intensitätsverhältnis der 
Strahlung zweier schwarzer Körper (logarith- 
misch) mit den Temperaturen 7, und 7. 

E 10-4 cm 
dort meist genauer durchgefiihrt „=h2 w= 


E 
werden kann. Sie ist, nebenbei r 


bemerkt, identisch mit der Messung, die bei dem Isothermen-Verfahren (siehe 
oben) erforderlich ist. 


Aus der Gl. (6) folgt nun weiter: 
— wr-)-Tgdwr, = Ig 


_ ighwr_ 

Ty Ar Tg -wr)' (9) 
Da t_, wie wir bereits gesehen haben, sich z. B. aus Messungen im Spektralgebiet 
der Wienschen Näherung leicht bestimmen läßt — siehe (8) —, erhalten wir einen 
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geschlossenen und mit Hilfe der einschlägigen Tabellen‘) leicht auswertbaren 
Ausdruck für r;. Freilich muß dabei beachtet werden, daß tatsächlich nicht n, 
sondern 7, gemessen wird. Nur wenn ¢, = &, ist, ist 7 = 7. Es ist somit im all- 


gemeinen notwendig, In “2 besonders zu bestimmen. Aber auch diese Zahl 


kann aus dem Verlauf von », durch Extrapolation der Wienschen Asymptoten 


ermittelt werden. Es gilt ja: 


(Neo = In + oo © In 2 + 2t_w. 
1 1 


Bringt man die Asymptote (1e)w+o zum Schnitt mit der 7,-Achse, dann 


5 
“ 
f 
2 a : 
fy 
7 
Y CH: 20 fı 
0 1 2 


Abb. 2. Der Einfluß von Meßfehlern auf das 
Ergebnis der Temperaturberechnung. 

Ar_ Av 


Die Kurven a und 6 entsprechen den Kurven- 

scharen der Abb. 1. Der schraffierte Bereich, der 

die f,-Werte begrenzt, entspricht dem Gebiet B 
in Abb. 3 


Eintr, 


ist der Achsenabschnitt In. 


1 
Aus tT, und r_ ergibt sich 
nun in einfacher Weise: 


_ 14880 . 
= 
T, = 14880 °K. 
§ 5. Der Einflu8 von MeBfehlern 


Da infolge der Eigenart des 
Verfahrens ein erheblicher Ein- 
flu8 von MeBfehlern auf das 
Ergebnis zu erwarten ist, ist die 
folgende Fehlerabschätzung von 
besonderer Wichtigkeit. Ein rela- 

Ar_ 


tiver Fehler gr bzw. = bei der 
Bestimmung der Größe von Ar_ 
bzw. v =; wirkt sich auf r, 
und damit auf 7,, in folgender 
Weise aus: 


mit 


) (11a) 


Ww 
Tw + Ga (In v —t_- w) 


Sint, w 1 


0 


(11b) 


14) K. Hayashi, Fünfstellige Tafeln der Kreis- und Hyperbelfunktionen sowie der 
Funktionen exp (x) und exp (—x), Berlin 1944. 
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Da bei Anwendung dieser Gleichungen w stets klein ist, können wir als Näherung 
und obere Grenze für f, und f, angeben: 


(1 + w) (12a) 
t_\2 
_ 
ty I (12b) 


Diese Näherungen sind in Abb. 2 dargestellt. Wie man erwarten muß, steigt der 
Einfluß von Meßfehlern sehr rasch mit größer werdenden r, - w-Werten, d.h. 
mit kleineren Temperaturen 

und mit kleineren Wellen- °K Q,° 


lingen. Im Diagramm der 20000 - 
Abb. 3 werden die Tempe- GY Z 3 N 
raturbereiche gezeigt, die hier- ° 

nach überhaupt für eine An- VB Gh AN 
wendung der Methode in T> 768 N 

Frage kommen. Im Gebiet A Yi Ke 
Größenordnungen schlechter 10000 5 & > 

als die eigentlichen Meßfehler. 
Innerhalb des Gebietes B yy Y 
geht die Größenordnung des b, Yj 
Fehlers des Ergebnisses nicht & 
über die Größenordnung der 

Meßfehler hinaus. Soweit man 


wird die Unsicherheit des 
Auswertungsergebnisses um 
im Übergangsgebiet Messun- 
gen macht, darf man keine 


NN. 


OL LLL 


genauen Ergebnisse erwarten. 10000 20000 °K 
1 
Abb. 3. Der Anwendungsbereich der Methode. 
§6. Über die — A: Anwendbarkeit ausgeschlossen. 
der Methode B: Die Unsicherheit der berechneten Temperatur 
Danach sind alle Versuche, liegt in der Größenordnung der Fehler der 


Strahlungsmessungen. 
Temperaturen von Hohlraum- pas Gebiet B entspricht dem schraffierten Gebiet 
strahlern oder von glühenden in Abb. 2. Die eingezeichneten Punkte entsprechen 
festen Körpern nach dieser den Kurven in Abb. 1 

Methode zu bestimmen, von 

vornherein grundsätzlich zum Scheitern verurteilt, und es ist verständlich, daß auch 
die Warburgsche Methode!) nicht die seinerzeit erhoffte Anwendung finden 
konnte, nämlich die Temperaturskala für hohe Temperaturen ausschließlich an das 
Plancksche Strahlungsgesetz anzuschließen und den mit einer kleinen Unsicher- 
heit behafteten Anschluß an die thermodynamisch definierte Temperaturskala 
für die Temperaturen bis zum — (1063° C = 1336° K)**) zu ver- 
meiden. 


15) E. Warburg, a. a, O. 
16) F, Kohlrausch, Praktische Physik, Bd. 1, S. 148, Leipzig 1950. 
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Selbst bei der Temperatur des positiven Kohlebogenkraters (etwa 4000° K) 
sind die Verhältnisse noch recht ungünstig. Vor allem ist die zeitliche und örtliche 
Konstanz der Kratertemperatur nicht ausreichend, um die Meßfehler hinreichend 
klein zu machen. Bekanntlich'!”) bildet der Krater nach Untersuchungen von 
Weizel, Faßbender, Finkelnburg und Schluge keine gleichförmig leuch- 
tende Fläche. Vielmehr zieht sich der Kohlebogen wie andere Lichtbögen zu einem 
anodischen Brennfleck zusammen, der sich innerhalb des sogenannten Bogen- 
ansatzes mit sehr großer Geschwindigkeit in statistisch unregelmäßiger Weise 
bewegt. Alle Messungen der Leuchtdichte vermitteln also lediglich ein zeitliches 
und örtliches Integral über eine im einzelnen unbekannte Verteilung, die nach 
kinematographischen Untersuchungen beträchtliche Schwankungswerte aufweist. 
In der örtlichen — zeitlich gemittelten — Verteilung hat man auf Grund visueller 
Leuchtdichtemessungen!®) um den 
Docht herum ein Intensitätsmaxi- 
mum und nach dem Rand allmäh- 
lich abfallende Intensität festge- 
stellt. 

Um dies letztere, für die Beur- 
teilung der Leuchtdichtekonstanz 
des Kraters sehr wesentliche Ergeb- 
nis objektiv zu prüfen, wurde die 
Leuchtdichte quer über den Krater 
mit einem Registriergerät!?) ge- 
messen (Dauer der Abtastung des 
.- 10mm Kraters: etwa 2 Sekunden) und 
Dochtkohle MIRA. 108, aufbellung Betzage 
der Ringsdorff-Werke, —-Kohle: Homogen A, der mittleren Intensität gefunden. 
72, von Siemens-Plania, Stromstärke: 15 A Der Docht strahlte gleichfalls etwas 

(Nach einer Registrierkurve) heller als seine Umgebung (s. Abb. 4). 
Verwendet wurde eine Dochtkohle 
von 10 mm Durchmesser der Ringsdorff-Werke, Marke MIRA. Die Belastung 
betrug 15 A. Wir vermuten als Ursache der Randaufhellung, daß der Brennfleck 
bei seiner Bewegung den Kraterrand bevorzugt. Bei anderen Kohlen dürften die 
Verhältnisse ähnlich liegen. In qualitativer Übereinstimmung hiermit stehen 
amerikanische Messungen ®), bei denen für die hellste Stelle des Kraters pyro- 
metrisch die Farbtemperatur 3995° K, dagegen spektrographisch für einen kleinen 
Bezirk des Kraters eine mittlere Verteilungstemperatur von 3860° K gefunden 
wurde. 


Nach neuesten Messungen *!) treten zeitliche Temperaturschwankungen von 
+ 12° auf sowie Stromdichteschwankungen von +15%. Nach all dem ist der 


17), Weizel-Rompe, Theorie elektr. Lichtbögen und Funken, S. 103, Leipzig 1949. 
18) F,Patzelt u. K. Baldewein, Wiss. Veröff. Siemens-Konzern 21, 353 (1943) 
(Bild 1, 2 u. 6), 
. - Beschreibung des Registriergeräts in einer demnächst erscheinenden Arbeit des 
erf. 
20) H. G. MacPherson, a. a. O. 
#1) J. Euler, Z. ang. Phys. 8, 260 (1951). 
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positive Kohlebogenkrater zur Heranziehung der oben beschriebenen Methode 
nicht geeignet. Ihr Anwendungsgebiet beginnt erst bei den hohen Temperaturen 
oberhalb aller bekannten Siedepunkte, also bei Temperaturen, wie sie in Licht- 
bögen und in Sternatmosphären vorkommen. In der Tat kann die spektrale Ver- 
teilung solcher Objekte in vielen Fällen sehr gut durch eine Plancksche Verteilung 
angenähert werden*), so daß dort von einer Verteilungstemperatur gesprochen 
werden kann. Bei astrophysikalischen Objekten entfällt durch geeignete An- 
wendung unserer Methode die Notwendigkeit, die Verteilungstemperatur über 
irgendeinen Anschlußstern an eine irdische Strahlungsquelle anzuschließen und 
damit die Unsicherheit des spektralen Verlaufs der atmosphärischen Extinktion 
in Kauf zu nehmen. 


2) §. a. P. Schulz, Z. Naturforschg. 2a, 583 (1947); Ann. Physik (6) 1, 95 (1947). 


Greifswald, Institut fiir Gasentladungsphysik der Deutschen Akademie der 
Wissenschaften zu Berlin. 


(Bei der Redaktion eingegangen am 15. Februar 1952.) 
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Die Dielektrizitätskonstanten heterogener Mischkörper 
aus isotropen und anisofropen Substanzen 


Von Walter Niesel 


Inhaltsübersicht 


Es wird eine allgemeine Theorie zur Berechnung der D.K. (Dielektrizitäts- 
konstanten) heterogener Gemische aus isotropen und anisotropen Substanzen 
aufgebaut. Die Einzelteile der gemischten Aggregate (Elemente) werden durch 
Ellipsoide angenähert. In den Sonderfällen des Kugel- und Lamellenaggregates 
ergeben sich die bereits von Bruggemann!) abgeleiteten Formeln. 


Nomenklatur 
7 — rel. Dielektrizit. Konst. im Inneren von Ellipsoiden 
Eq = rel. Dielektrizit. Konst. außerhalb von Ellipsoiden 
€’ — resultierende Dielektrizitätskonstante 
a,b,c = Hauptachsen von Ellipsoiden 
E,,, Ey, Ez,\ x, y, Komponenten der Feldvektoren des homogenen, unge- 
D,,, Dya; Dz,| störten Feldes außerhalb der Ellipsoide 
Ze Die entsprechenden Komponenten im Inneren der Ellipsoide 

L, M,N = Entelektrisierungsfaktoren 

6; — Volumenanteil eingesprengter Ellipsoide 
6, = Volumenanteil des Mediums außerhalb der Ellipsoide 
Eu Ey, Ey = resultierende D. K. bei Feld in Richtung der x-, y-, z-Achsen der 


Ellipsoide. 


I. Bedeutung des Problemes 


Bruggemann löste das oben angegebene Problem für die Sonderfälle des La- 
mellen- und des Kugelaggregates. In der Literatur seit 1936 sind nur wenige 
Hinweise auf die Bruggemannsche Veröffentlichung enthalten, und die Er- 
gebnisse der Theorie sind kaum angewendet worden. Der Hauptgrund hierfür 
liegt wohl hauptsächlich darin, daß für D. K.-Unterschiede bis zu einer Zehner- 
potenz die schon lange bekannten empirischen Mischungsformeln von Wiener 


aus isotropen und anisotropen Substanzen“, Ann. Physik 5, 24, 636 (1935). Bruggemann, 
„Dielektrizitätskonstanten und Leitfähigkeiten von Vielkristallen der nichtregulären 
Systeme‘ Ann. Physik 5, 25, 645 (1936). : 


1) Bruggemann, „Dielektrizitätskonstanten und Leitfähigkeiten der Mischkörper 
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und Liehtenecker?) ziemlich gut gelten. Da Bruggemann in seiner Arbeit 
zeigte, daß diese empirischen Formeln eine gute Näherung darstellen, lag keine 
Veranlassung vor, die zum Teil komplizierteren Bruggemannschen Formeln 
zu benutzen. Ein weiterer Grund ist darin zu suchen, daß die Annäherung durch 
Kugeln oder Lamellen für manche Fälle zu grob ist. 

In den letzten Jahren traten jedoch bei der Untersuchung seignette-elektri- 
scher Substanzen Probleme auf, die mit Hilfe der Bruggemannschen Formeln 
hätten gelöst werden können. Rushman?) leitete z. B. für das Mischungssystem 
„Luft-Bariumtitanat‘, das durch mehr oder weniger starke Sinterung von BaTiO, 
entsteht, eine Näherungsformel ab. Die von Bruggemann für diesen Fall abge- 
leiteten Formeln gelten jedoch, solange die Voraussetzungen für die Anwendung 
erfüllt sind, sehr exakt, wie in Abschnitt V, näher ausgeführt wird. 

Andererseits kann man gesintertes BaTiO, ohne Lufteinschlüsse als einen Viel- 
kristall auffassen, dessen Elementarkristalle stark anisotrop sind. Es liegt nahe, 
die entsprechenden Bruggemannschen Formeln zur Berechnung des Vielkri- 
stalles anzuwenden. Es zeigt sich, daß die Bruggemannschen Formeln die 
experimentellen Ergebnisse nur unbefriedigend wiedergeben. Bei Anwendung 
der im folgenden entwickelten allgemeineren Formeln kann man jedoch die ex- 
perimentellen Ergebnisse gut erklären und darüber hinaus eine Reihe weiterer 
Erscheinungen mühelos deuten, z. B. das Altern des BaTiO,, die Auswirkungen 
verschiedener Brenntemperaturen, die Wirkung von Verunreinigungen (siehe 
Abschnitt V). 

Bei Anwendung der Bruggemannschen Formeln und der hier abgeleiteten 
Formeln müssen bestimmte Voraussetzungen erfüllt sein, die im nächsten Ab- 
schnitt angeführt werden. In Abschnitt V soll an Hand einiger Beispiele gezeigt 
werden, wann die einzelnen Formeln angewendet werden dürfen. 


II. Grundlagen 
a) Voraussetzungen 


Bei der Berechnung der D. K. einer Mischung von zwei Substanzen sind drei 
typische Fälle zu unterscheiden: 

1. Die beiden Substanzen bilden eine homogene Mischung, stellen also eine 
einzige Phase dar. 

2. Kleine Teilchen der ersten Substanz (Elemente) sind in die zweite Substanz 
eingesprengt. Die zweite Substanz bildet eine homogene Phase. 

3. Beide Substanzen bestehen aus einzelnen kleinen Teilchen. Diese Elemente 
der beiden Substanzen erfüllen den gesamten Raum. (Das ist dadurch möglich, 
daß die Zwischenräume zwischen den einzelnen Elementen jeweils durch kleinere 
Elemente ausgefüllt werden.) 

Zu 1. Die Berechnung der D. K. ist in diesem Falle ein molekulartheoretisches 
Problem und wird hier nicht behandelt. 

Zu 2. Dies ist der am häufigsten auftretende Mischungsfall. (Z. B. Emul- 
sionen, Kristallpulver in Luft, keramische Stoffe mit Lufteinschlüssen u. a.) 


2) Lichtenecker, „Die D. K. natürlicher und künstlicher Mischkörper‘‘, Physik. Z. 
27, 115 (1926). 

3) Rushman, „The effective Permittivity of two-phase-Systems‘‘, Proc. physic. Soc. 
London 59, 1011 (1947). 
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Zu 3. Dieser Mischungsfall ist selten ideal verwirklicht. (Etwa durch zwei im 
Eutektikum auskristallisierte Substanzen.) Unter 3. fallen jedoch die Vielkristalle 
anisotroper Stoffe. Die anisotropen Vielkristalle stellen sozusagen eine Mischung 
mit sich selbst dar. 

Um die Fälle 2. und 3. berechnen zu können, müssen die Elemente (Einzel- 
bestandteile) der Substanzen bestimmte Voraussetzungen erfüllen, die von Brug- 
gemann angegeben wurden: 

«) Die Elemente der gemischten Körper enthalten sehr viele Atome bzw. 
Moleküle. Jedes Element ist homogen isotrop oder anisotrop. 

ß) Alle Elemente sind sehr klein gegen die Dimensionen des Aggregates. 

y) Die Elemente sind regellos verteilt. 

Jedes den Voraussetzungen «), 8), y) genügende Aggregat wird im folgenden 
als quasihomogen bezeichnet. Wird an dieses Aggregat ein Feld gelegt, so be- 
zeichnen wir dieses Feld als quasihomogenes Feld. Die gemittelten E- und D- 
Komponenten über Raumelemente, groß gegenüber den Struktureinheiten, be- 
sitzen überall im quasihomogenen Aggregat dieselben Werte. 


b) Hypothetische Grundsätze 


Bruggemann stellt vier Grundsätze auf, um seine Formeln abzuleiten. 
Von den vier Grundsätzen wird jedoch im folgenden nur einer benötigt, den ich 
in abgeänderter Fassung wiedergebe: 

Streut man in ein homogenes Medium mit der D. K. e, ein Aggregat von ein- 
ander ähnlichen Elementen mit der D. K. e, ein (z. B. Ellipsoide), die alle dieselbe 
Orientierung besitzen (bei Ellipsoiden parallele Achsen), so besitzt der Misch- 
körper eine neue bestimmte D. K. e’. 

Streut man andererseits genau dasselbe Aggregat in derselben Konzentration 
mit derselben D.K. e, und derselben Orientierung in ein quasihomogenes 
Medium, mit der nach außen makroskopisch gemessenen D.K.e, ein, so besitzt 
der Mischkörper dieselbe resultierende D. K.e’ wie im ersten Fall. 


Auf drei der vier Bruggemannschen hypothetischen Grundsätze kann ver- 
zichtet werden, da bei der Ableitung der Einstreuungsformeln für kleine Einstreu- 
ungen ein neuer Weg eingeschlagen wird. 


IH. D.K.-Änderung bei Einstreuung kleiner Volumenanteile von einander 
ähnlichen Ellipsoiden 


a) Die D. K. eines quasihomogenen Aggregates kann man aus der Kapazität 
eines ebenen Plattenkondensators bestimmen, der als Dielektrikum dieses quasi- 
homogene Aggregat besitzt. In dem unendlich ausgedehnten Kondensator denken 
wir uns ein räumliches orthogonales Achsensystem, dessen z-Achse senkrecht zu 
den Kondensatorplatten steht. Die Kondensatorplatten selbst liegen in der x—y- 
Ebene oder parallel dazu. Es ist bei Benutzung des Gaußschen Maßsystemes 
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An Stelle von ff D,dF kann D,F geschrieben werden. Die Mittelung hat dabei über 
eine Fläche parallel zu den Kondensatorplatten zu erfolgen. Da die Mittelung über 
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jede parallele Fläche denselben Wert ergibt, kann die flächenmäßige auch durch 
eine räumliche Mittelung ersetzt werden. Die räumliche kann unter Umständen 
bedeutend leichter durchzuführen sein als die flächenmäßige Mittelung. 


Für f E, dz kann E,d gesetzt werden. Da es gleich ist, an welcher Stelle 


0 
(x, y) die Mittelung erfolgt, kann die lineare wieder durch die räumliche Mittelung 
ersetzt werden. 
Wir erhalten: 


1 
; 
4x E,d 
Durch Vergleich mit der Kapssitätsformel für einen Plattenkondensator mit 
homogenem Medium erkennt man: 
= 1 
(1) 
Dabei soll sich die Mittelung der D,- und E,-Komponente über das gesamte Konden- 
satorvolumen erstrecken. 

b) Bringt man ein Ellipsoid in ein weit ausgedehntes homogenes elektrisches 
Feld, so wird das Innere des Ellipsoides homogen elektrisiert. Wir wählen ein 
kartesisches Koordinatensystem so, daß die Hauptachsen des Ellipsoides mit den 
Koordinatenachsen zusammenfallen. Die Ellipsoidfläche ist dann gegeben durch: 


2% 


Der Zusammenhang der Feldstärken im Innenraum und im Außenraum der 
Ellipsoide ist *): 
E,, = E,,(1 + D) 


Ey, = Ey, (1 + x%,M) (2) 


E,, = E,,(1+%,N). 


Gleichungen (2) gelten nur dann, wenn die Hauptachsen des Ellipsoides mit den 
Hauptachsen des e-Tensors zusammenfallen. 

L, M, N, sowie x haben hier andere Bedeutungen als im Winkelmannschen 
Handbuch. Die Bezeichnungsweise richtet sich nach Sommerfeld 5), siehe vorn 
unter Nomenklatur. 

In Sonderfällen ist: 

Kugel: L = M = N = }. 

Kreiszylinder; Rotationsellipsoid mit unendlich langer Rotationsachse, 

Zylinderachse in z-Richtung: L = M = 4, N=0, 

Lamelle: Ellipsoid mit einer verschwindenden Achse, Lamellennormale in 
z-Richtung: L = 1, M=N=0. 

c) Es soll jetzt die D. K.-Änderung berechnet werden, die durch Einstreuung 
kleiner Mengen einander ähnlicher, achsenparalleler Ellipsoide entsteht, wenn ein 
Feld in z-Richtung anliegt, die c-Achsen also in z-Richtung weisen. Wir benutzen 


4) Winkelmann, „Handbuch der Physik‘, Bd. V, S. 124, Leipzig 1908. 
5) Sommerfeld, „Lehrbuch der theoretischen Physik“, Bd. III, rs 99, Leipzig 1949. 
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dabei die in IIIa abgeleitete Beziehung (1) 


e = 

Zur Berechnung der mittleren E,- und D,-Komponenten zerlegen wir das Feld 
der Ellipsoide in zwei Teile: 

1. Ungestörtes homogenes Feld der Feldstärke E,, (im folgenden nur mit E, 
bezeichnet), das sowohl im Inneren als auch außerhalb der Ellipsoide anliegt. 

2. Störfeld, das durch Oberflächenladungen zustande kommt. Diese entstehen 
durch die Unterschiede der Polarisationen im Innern und außerhalb der Ellip- 
soide. Das Störfeld ist innerhalb der Ellipsoide homogen. Seine Größe ist dort: 


Est, = Ey— Eg, (Est, = Feldstärke des Störfeldes im Inneren), 
wobei E,, durch Formel (2) gegeben ist, so daß wird 
E, x,N 
Ese, = 1+x,N 


Da die Summe der Oberflächenladungen, die das Störfeld erzeugen, Null ist, muß 
die mittlere E,-Komponente des Störfeldes ebenfalls Null sein. Wenn wir mit 
6, den Volumenanteil der Ellipsoide und mit 6, den Volumenanteil des homogenen 
Mediums bezeichnen (6; + 6, = 1), so gilt daher die Beziehung: 


6; Est, + é, Est, = 0 


6, Est. = — 6; Est,- 
Jetzt lassen sich die mittleren E,- und D,-Komponenten leicht angeben: 


E, ist gleich der ungestörten Feldstärke E, 
D, = + + 6,2, 18, + Estq] 
= 6,¢,E, + 6;¢, Est, + 6,) —€, 6; Est, 
= 6, + 6; Eq + &; Est;] 
7 7 — &4 
¢, E, + 8, (€,—&) Eq [1— |= E, [1 


€ 5, 1 


Verallgemeinerungen: Je nachdem in welcher Achsenrichtung des Ellipsoides 
das Feld angelegt wird, erhält man: 


= e,(1 + =e, (1+ = x): 


Liegt das Feld in einer Richtung an, die die Winkel a, 8, y mit der x, y, z-Achse 
einschlieBt, so wird wegen des tensoriellen Charakters der D. K. 


[1 +6, (cost + cos? B 1 + cos*y 5 y)] (3) 
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Bei regelloser Verteilung der Ellipsoide muß räumlich gemittelt werden: 


ER 
(3a) 
Im folgenden will ich zur Abkiirzung setzen: 


%, 
MM. 


Die Einsprengungsformel nimmt dann die Gestalt an: 
=e, (14+ 46,F (x, L, M, (3c) 


IV. Ableitung der allgemeinen Mischungsformeln 


Formel (3c) gilt nur für sehr kleine Einsprengungen exakt. 

Die Aufstellung der Mischungsformeln für endliche Einsprengungen und 
Mischungen erfolgt jetzt genau nach dem bereits von Bruggemann angegebenen 
Verfahren — nur allgemeiner —, wobei der unter IIb) angegebene Grundsatz 
verwendet wird. 


a) Differentialgleichung für die Einsprengung in ein homogenes Medium 


ö, sei der Volumenanteil der eingesprengten Substanz, 6, der Volumenanteil 
der homogenen Substanz. 


Um die D. K. dieser Mischung zu berechnen, nehmen wir an, die homogene 
Substanz sei eine Fliissigkeit, in die die Elemente der zweiten Substanz in kleinen 
Portionen eingesprengt werden. Um eine Volumeneinheit der Mischung herzu- 
stellen, miissen wir in 6, Volumeneinheiten der homogenen Substanz 6, Volumen- 
einheiten der 2. Substanz einsprengen. 

Angenommen, wir haben zunächst nicht die gesamte Menge von 6; Volumen- 
einheiten eingestreut, sondern nur die Menge von 6 Volumeneinheiten. Dann ist 
das Volumen der Mischung (6, + 6)-Volumeneinheiten. Die D. K. dieser Ein- 
sprengung sei 

Nehmen wir jetzt eine weitere Einsprengung von dö Volumeneinheiten vor, 
so sind das im Verhältnis zum augenblicklichen Volumen bei dieser Einstreuung: 


Volumenanteile. 


Die resultierende D. K. nach dieser kleinen Einsprengung wird nach Gl. (3e): 


€ =e[1 + 


dé 
L, M,N) 
oO F (x, L, M,N). 
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Dies ist eine separierbare Diffgl. für & mit der Anfangsbedingung ¢ = & 6=0, 
6; "Res 

0 &q 
€Res ist die resultierende D. K. Sie wird im folgenden, wenn keine Verwechslung 
möglich ist, mit ¢ bezeichnet. 

Durch die Transformation 6, + 6 = 1— 6’ erhält man: 


oj e 
dé’ de 
&a 

Durchrechnung von drei Sonderfällen bei Einstreuung isotroper Substanzen 
ergibt: 

1. Kugelaggregat: 

L=M=N=} = ty, =, =& 
F (x, L, M, N) = ‘ 

Einsetzen in Formel (4) und Ausintegrieren ergibt: 


1-4, = (5) 
Dies ist die bereits von Brugge mann für diesen Fall berechnete Formel. 
2. Lamellenaggregat: 
L=1,M=N=0, 


F (x, L, M,N) = x[ + 2] )- 


Einsetzen in Formel (4) und Ausintegrieren ergibt: 


+ 2 6; (&; — &a) 


€ 
&—& 2, 


3. Kreiszylinderaggregat: 
L=0,M=N=} 


P(x, L, M,N) = x[1 + ral 
— 4 Ze (e+ 


Die entsprechenden Einsprengformeln fiir anisotrope Substanzen werden sehr 
unhandlich. Fiir ein elektrisch einachsiges Kugelaggregat habe ich die Differential- 
gleichung gelöst. 

Die hiernach berechneten D. K.-Werte erwiesen sich für alle praktisch vor- 
kommenden Fälle als fast identisch mit den D. K.-Werten, die man erhält, wenn 
man in die gewöhnliche Kugeleinsprengformel für e, die mittlere D. K. des massiven 
eingesprengten Stoffes setzt. 
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b) Ableitung der Gleichungen fiir die Mischung zweier inhomogener Substanzen 


Gegeben seien zwei Aggregate, die jedes für sich nur aus einander ähnlichen 
Ellipsoiden bestehen. 


Aggregat 1: &,,, Ey, &2,, Zu M,, Nı- 
Aggregat 2: &,, Ey, &2,, Lu My, Nz. 


Diese Aggregate mische ich mit den Volumenanteilen 6,, 6, (6, + 6, = 1) 
so zusammen, daß nur noch vernachlässigbar kleine Zwischenräume übrig bleiben. 
(Siehe Bemerkung zu Ila, 3). Das-so gebildete quasihomogene Aggregat möge 
die D. K. e besitzen. Dasselbe Aggregat kann ich auch dadurch herstellen, daß ich 
in eine zunächst homogene Substanz mit der D. K. e abwechselnd von der Sub- 
stanz „l‘“ und von der Substanz ‚2‘ kleine Mengen einsprenge, bis die gesamte 
homogene Substanz verdrängt ist. Nach der ersten kleinen Einsprengung 6; der 
Substanz „1“ besitzt das Aggregat nach Gl. (3c) die neue D. K.: 


fl 4 (%4, L. M,, N,)] 
mit 
x,y, 


Jetzt sprenge ich soviel von der zweiten Substanz in das nun quasihomogene 
Feld ein, daß wieder die ursprüngliche D. K. entsteht. Nach dem in IIb auf- 
gestellten Grundsatz läßt sich die resultierende D. K. allgemein nach denselben 
Formeln berechnen wie bei der Einstreuung in ein homogenes Medium. Hier wird: 


e=e[l+ 4 F (xs, L,, 
{1 +46, F Ly, M,, {1 + 462 F (x2, Ly, M,, N,)}. 


Hier bedeutet von x beispielsweise die «-Komponente ae (vgl. die Nomen- 


klatur am Anfang). Durch weitere Reihenentwicklung und Vernachlässigung von 
Gliedern höherer Ordnung erhält man die Gleichung 
0 6, F (%1, M,, N,) + (#2, L,, M, N,). 
Dieselbe Gleichung gilt nach unserem Grundsatz aus IIb für jede weitere Ein- 
sprengung 64’, 62’ usw., also 
(6; + er + éy” ‘)F Ly, M,, N,) + (63 + + 6,” + AR 
F L,, Mz, N;). 
Alle Teileinsprengungen der ersten Substanz ergeben zusammen Ö,, und alle 


Teileinsprengungen der zweiten Substanz ergeben zusammen 6,. Nach Voraus- 
setzung muß dann das Aggregat genau die D. K. & besitzen, so daß mit Gl. (3b): 


L ! 1+%, (8) 
Hy )=o 


Setzt man für die x ihre Bedeutung ein, so läßt sich e ausrechnen. 
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Die Durchrechnung für drei Sonderfälle bei isotropen Substanzen ergibt: 
1. Kugelaggregat: 


e=}[(— + 2E,) + VETER + 8&6] (9) 
E, + &2 
E, = 6, 4 Oe &- 
2. Lamellenaggregat: 


mit 


+ 2m) + + 2m) = 0 


e=}(—£,+ 8E,E,) (10) 
mit 
1_ |, 


3. Kreiszylinderaggregat : 
Bet + (5 4 E,) e* — (b, + 4.8 + 6,63) €— E, = 0. (11) 


Die Formeln fiir das Kugel- und das Lamellenaggregat stimmen wieder mit 
den von Bruggemann abgeleiteten Formeln iiberein. 


€) Gleichungen zur Berechnung eines Vielkristalles 
Wenn man in der Gl. (8) 6, gleich Null setzt, stellt der Rest der Gleichung die 
Formel zur Berechnung des Vielkristalles dar. (Siehe Bemerkung zu Ila, 3.) 
Ich will wieder einige Sonderfälle durchrechnen, zunächst mit völliger Aniso- 
tropie: 
1. Kugelaggregat: 
Ay %; 


(12) 


%,, %, und x, sind nach Definition (vgl. III b) Funktionen von e,. Eine Auflösung 
nach ¢, ist möglich. Da ¢, die mittlere D. K. im Außenraum der Elemente des 
Mischkristalles ist, ist es gleich der resultierenden D. K. 


2. Lamellenaggregat, Lamellennormale parallel zur x-Achse: 


+ + = 9. (13) 
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3. Kreiszylinderaggregat, Zylinderachse parallel zur x-Achse 


Hy 
(14) 


Besonders wichtig ist der Fall des dielektrisch einachsigen Kristalles. Die 
Formeln vereinfachen sich, wenn man folgende Bezeichnungen einfiihrt: 


&=& & =e, (1—1). 
4. Kugelaggregat 
e= }e,(1+ /9—8r). (15) 
5. Lamellenaggregat, z-Achse senkrecht zur Lamellenoberfläche: 
L=M=0,N=1 
e = fe,[c—1 + 9—2). (16) 


6. Lamellenaggregat, z-Achse parallel zur Lamellenoberfläche: 


L=f, B= 
e=}e,[—1+ (17) 
7. Kreiszylinderaggregat, z-Achse parallel zur Zylinderachse: 
N=0 
+ + 2]. (18) 
8. Kreiszylinderaggregat, z-Achse senkrecht zur Zylinderachse: 
L=0, 9=N=} 
Xp 


Die Auflösung nach e wird kompliziert. 


V. Diskussion an Hand von Beispielen 


Um mit den aufgestellten Formeln die D.K. einer Mischung zu berechnen, 
muß man sich zunächst darüber Klarheit verschaffen, welche der Formeln an- 
gewendet werden muß, und in wieweit eine Annäherung des wirklichen Aggregates 
durch Kugeln, Lamellen oder andere Ellipsoide gestattet ist. An Hand von einigen 
Beispielen sollen die auftretenden Probleme erörtert werden: 

1. Rushman®) führte D. K.-Messungen an verschieden stark gesinterten 
BaTiO,-Proben durch. Durch Variation der Sinterungstemperatur und der 
Sinterungsdauer wurden verschiedene Porositäten erreicht. Wie aus den Angaben 
von Rushman hervorgeht, sind die Poren etwa kugelförmig. Bei nicht zu großer 
Porosität sind die einzelnen Hohlräume nicht miteinander verbunden, es liegt also 
die Einsprengung eines kugelförmigen Aggregates von niedriger D.K. in ein 
homogenes Medium von hoher D. K. vor. Für diesen Fall gilt die Formel (5). 

Bis zu 6; = 0,3 gibt die Formel sehr gut die D. K. der porösen Substanz an, 
wie sie von Rushman gemessen wurde. Für 6; > 0,3 beginnen die einzelnen Poren 
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zusammenzuhangen. Es liegt dann ein komplizierteres Aggregat aus Kugeln und 
mehr oder weniger langgestreckten Ellipsoiden vor. Eine einfache Berechnung 
ist nicht mehr möglich. 


2. Bei eigenen Versuchen wurde BaTiO,-Pulver mit Paraffin gemischt und die 
resultierende D. K. bestimmt. Die Elemente des Pulvers besaßen körnige Struktur 
verschiedener Gestalt. Es liegt wieder die Einsprengung eines Dielektrikums in 
ein homogenes Medium vor. Jetzt besitzt jedoch das eingesprengte Aggregat 
keine Ellipsoidgestalt mehr. Das Experiment zeigt, daß innerhalb der Fehler- 
grenze für die Einsprengung des BaTiO, in Paraffin die Kugeleinsprengformel (5) 
gilt, wenigstens auf 3%. (Die Anisotropie des BaTiO, braucht dabei, wie in IVa 
ausgeführt wurde, nicht berücksichtigt zu werden.) Eine endgültige Entscheidung, 
wie weit die Einsprengung eines Körneraggregates, dessen Elemente nicht stark 
von der Kugelform abweichen, mit der Kugeleinsprengformel berechnet werden 
darf, konnte mit den zur Verfügung stehenden Mitteln noch nicht getroffen werden. 


3. Versucht man, die Einsprengformeln (5) für gepreßtes BaTiO,-Pulver in 
Luft anzuwenden, so erhält man keine gute Übereinstimmung mit dem Experi- 
ment. Gepreßtes BaTiO, besitzt eine höhere D. K. als nach den Formeln zu er- 
warten wäre. Die Gültigkeit von Formel (5) für diesen Fall ist auch gar nicht zu 
erwarten, denn die einzelnen BaTiO,-Elemente sind nicht mehr voneinander 
isoliert, sie berühren sich mit größeren Flächen und bilden miteinander dielektrisch 
leitende Ketten. Im Grenzfall bei sehr starkem Druck bildet das BaTiO, ebenfalls 
eine geschlossene Phase. Die Einsprengformeln gelten dann sicher nicht mehr. 


4. Besonders interessante Ergebnisse erhält man bei Anwendung der Formeln 
aus [Vc auf den BaTiO,-Vielkristall. Zunächst muß man sich wieder darüber klar 
werden, wie weit sie zulässig sind. Wenn die hohe D.K. auf eine analoge Weise 
zustande käme wie die hohe Permeabilität ferromagnetischer Vielkristalle, würden 
die Formeln aus IVe sicher nicht angewendet werden dürfen. 


Da jedoch eine 90°-Drehung des Polarisationsvektors in Form einer Wand- 


verschiebung oder eines Sprunges beim BaTiO, mit einer großen Änderung der 
elastischen Energie des betreffenden Bezirkes verbunden ist (bedeutend größer als 
bei ferromagnetischen Stoffen), spielen Wandverschiebungen und 90°-Sprünge 
bei der Bildung der D.K. eines Vielkristalles keine oder nur eine sehr geringe 
Rolle. Zur Anregung von 180°-Sprüngen ist eine bestimmte Mindestfeldstärke 
erforderlich, so daß 180°-Sprünge bei kleinen Feldstärken keinen Beitrag zur Bil- 
dung der D.K. leisten. Die Anfangs-D. K. (mit verschwindend kleiner Wechsel- 
spannungsamplitude gemessene D.K.) des BaTiO,-Vielkristalles ist daher nur 
durch die Anfangs-D. K. des Einkristalles bestimmt. 


Nimmt man als D. K. des Einkristalles die von W. J. Merz *) gemessenen Werte 
bei 25°C: 
¢=200 —— = 0,955, 


so erhält man als resultierende D.K. bei Annahme der folgenden Aggregate: 


1. Kugelaggregat: 
€ = 2450 aus (15) 


6) W. J. Merz, ,,The electric behavior of BaTiO, Single-Domain-Crystals“, Physic. 
Rev. 76, 1220 (1949). 
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2. Lamellenaggregat mit der z-Achse senkrecht zur Lamellenoberfläche: 


€ = 1260 aus (16) 
3. Lamellenaggregat mit der z-Achse parallel zur Lamellenoberfläche: 
€ = 2870 aus (17) 


Man sieht aus diesen Beispielen, daß die Form der Kristallite ganz wesentlich 
die makroskopisch gemessene D. K. beeinflußt. Auf Grund dieser Tatsache lassen 
sich wahrscheinlich verschiedene Phänomene der seignette-elektrischen Stoffe 
theoretisch begründen, für die man bisher nur ungenügende oder keine Erklärungen 
besaß. 

VI. Ausblick 


1. Texturen. Wenn die Ellipsoide nicht regellos verteilt sind, sondern gegenüber 
der Feldrichtung Vorzugsorientierungen besitzen, gilt die Mittelungsformel (3a) 
nicht mehr. Es sei 


cota=p  cos*B=q  cos*y=r. 


An Stelle der Einsprengformel (3a) erhält man die Formel: 


Aus dieser Formel kann man sich ohne Schwierigkeit die Mischungsformeln für 


die einzelnen Fälle berechnen. Für die Mischung zweier inhomogener Phasen er- 
hält man z. B.: 


2. Einsprengung von Ellipsoiden, deren D, K.-Hauptachsen nicht mit den 
Hauptachsen der Ellipsoide zusammenfallen. 
An die Stelle der Formeln (2) treten die folgenden, von deren Ableitung ich 
hier der Platzersparnis wegen absehe: 


Ey, = Ex, (1 + %2L) + By + By 
By, = + By (1 + %yM) + (22) 


Eu = Ey, + Ey, +E,(1+x,N). 


Die Formeln (22) lassen sich für Spezialfälle nach E,,, E,, und E,, einfach auf- 
lösen. Für diese Fälle können einfache Mischungsformeln aufgestellt werden. 

3. Die hier abgeleiteten Mischungsformeln können ohne Änderung auf alle die- 
jenigen stationären Vorgänge übertragen werden, die analogen Differential- 
gleichungen und Grenzbedingungen wie E und D unterworfen sind. Dazu ge- 
hören: Elektrische Strömung, Wärmefluß, magnetischer Fluß. Die entsprechenden 
Konstanten, für die die Mischungsregeln gelten, sind: 


Elektrische Strömung: spezifische Leitfähigkeit o 
Wärmefluß: Wärmeleitzahl x 
Magnetischer Fluß: Permeabilität u 
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Bei der Anwendung der Mischungsregeln auf den elektrischen FluB und den 
Wärmefluß müssen die Aggregate noch Zusatzvoraussetzungen erfüllen. Bei der 
elektrischen Strömung dürfen an der Oberfläche der Elemente keine Übergangs- 
potentiale, beim Wärmefluß keine Wärmeübergangswiderstände auftreten. 

Über die Ausdehnung derselben Methode auf Probleme der Elastizität hoffe 
ich bald zu berichten. 


Für die experimentellen Untersuchungen konnte ich Apparate der Helmholtz- 


Gesellschaft verwenden, die Herrn Prof. Wolf von dieser zur Verfügung ge- 


stellt sind, wofür ich auch hier bestens danke. 


Karlsruhe, Institut für theoretische Physik der T. H. 


(Bei der Redaktion eingegangen am 21. Februar 1952.) 
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Zum ersten Randwertproblem der Maxwellschen Gleichungen 
Von E. Ledinegg und P. Urban 
(Mit 6 Abbildungen) 


Inhaltsiibersicht 


Die Anregung eines schwingungsfähigen elektromagnetischen Systemes, 
welches innerhalb eines bestimmten Raumteiles V gelegen ist, kann prinzipiell 
auf zwei Weisen erfolgen; entweder befindet sich der felderzeugende Generator 
innerhalb V und steht in unmittelbarer Wechselwirkung mit dem schwingenden 
System oder der Generator ist außerhalb V und die Anregung erfolgt mittelbar 
durch irgendein Hilfssystem. Im letzten Fall müssen längs F (wenn wir unter F 
die Hüllfläche von V verstehen wollen) entsprechende Koppel- oder Anregungs- 
flächen vorhanden sein, welche den Energiefluß vom Generator zum betrachteten 
System ermöglichen. 

Rein formal unterscheiden sich die beiden Fälle dadurch, daß bei Annahme 
einer „inneren‘ Anregung das elektrische Feld aus einer nichthomogenen Differen- 
tialgleichung zweiter Ordnung mit homogener Randbedingung zu berechnen ist, 
während im Falle der „äußeren‘‘ Anregung wir es mit einer homogenen Differen- 
tialgleichung mit nichthomogenen Randbedingungen zu tun haben. Im Gegen- 
satz zum erstgenannten Probleme, welches prinzipiell nichts Neues bietet und 
Gegenstand zahlreicher Untersuchungen war, ist das Problem mit äußerer An- 
regung, obgleich in der Hochfrequenztechnik eine Reihe derartiger Systeme ver- 
wendet werden (z. B. Hohlrohrresonatoren mit Hohlleitern als Energieträger), 
theoretisch kaum behandelt worden. Dies liegt vor allem daran, daß die Hoch- 
frequenzpraxis an exakten Lösungen mit kompliziertem analytischem Aufbau 
wenig interessiert ist und daher geeignete Näherungslösungen bevorzugt. Infolge 
dieses praktischen Bedürfnisses blieb die Frage offen, welche Bestimmungsstücke 
des elektromagnetischen Feldes im Sinne des ersten Randwertproblemes längs 
einer Hüllfläche frei vorgebbar sind. 

Die vorliegende Arbeit beschäftigt sich mit dieser, ein prinzipielles Interesse 
nicht entbehrender Problemstellung. Zunächst wird ein koordinatengebundenes 
Randwertproblem behandelt, um dann das praktisch interessierende Anregungs- 
problem in Angriff zu nehmen. Schließlich wird der Übergang zur koordinatenfrei 
formulierten ersten Randwertaufgabe vollzogen. Die Lösung des äußeren An- 
regungsproblemes wird auch in Form einer Integralgleichung angegeben. 


Verwendete Bezeichnungen 
&, u = Dielektrizitätskonstante bzw. Permeabilität in praktischen Einheiten. 
€, § = Amplituden der elektrischen bzw. der magnetischen Feldstärke. 
k*= em? = Kreiswellenzahl. 
Mit hochgestelltem * werden zweidimensionale Differentialoperatoren 
gezeichnet. 
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§ 1. Einleitung 


Die erste Randwertaufgabe einer vorgelegten Diff.-Gleichung 2. Ordnung 
D(u) = 0 besteht bekanntlich darin, eine reguläre Lösung u von D(u) = 0 inner- 
halb eines Gebietes V aufzusuchen, welche an der Begrenzungsfläche F von V 
(wir betrachten ein dreidimensionales Kontinuum) vorgeschriebene Werte an- 
nehmen soll. Im folgenden wollen wir dieses Problem für das System der Max- 
wellschen Gleichungen näher untersuchen. Wir nehmen in üblicher Weise die 
Zeitabhängigkeit in der Form e’®t an und erhalten etwa durch Substitution von 
9 die mit den Maxwellschen Gleichungen äquivalente Diff.-Gleichung 2. Ordnung 


rot rot E— 0 - (1) 


mit k? = eu @*. Im Gültigkeitsbereich V von (1) werden Quellen beliebiger Art 
ausgeschlossen. Im Sinne unserer Problemstellung ist zu ermitteln, welche Be- 
stimmungsstücke von € längs einer geschlossenen Fläche F frei vorgebbar sind, 
um eine eindeutig bestimmte zweimal stetig differenzierbare Lösung von (1) zu 
erhalten, welche am Rande das gewünschte Verhalten zeigt. Es ist sofort ersicht- 
lich, daß längs F nicht etwa der ganze Vektor (gemeint sind damit alle 3 Kom- 
ponenten von €) beliebig vorgebbar sein kann, denn wegen div € = 0 müssen z.B. 
die Beziehungen 


$E,df=0 bzw. rot, Edj= 0!) 
F 


gelten. Einen tieferen Einblick gewinnt man, wenn man statt des Systemes (1) 
das dazu äquivalente System 


AE+kRE=0 (2a) 
div € = 0 (2b) 
betrachtet. 
Legt man hier ein cartesisches Koordinatensystem zugrunde, so erhält man: 
AE,+kE,=0 (3a) 
_ 
(3b) 
oa, 
Aus der Theorie der skalaren Wellengleichung 
AU+kU=0 (4) 


wissen wir, daß durch Angabe von U längs einer geschlossenen Fläche F eine 
eindeutig bestimmte zweimal stetig differenzierbare Lösung von (4) existiert, 
sofern k? nicht gerade mit einem Eigenwert von 

AU+RU=0 mit U=0 langsF (5a) 
zusammenfällt. Überträgt man dieses Ergebnis auf (3) bzw. (2), so erhält man den 
folgenden Sachverhalt: Während das System (3a) 


AE E=0 (5b) 


unter Vorgabe von € längs F im allgemeinen?) eine eindeutige Lösung zuläßt, 
ist wegen der zusätzlichen Bedingung (3b), welche durch die or, ein- 


1) Mit der Vorgabe von € längs F ist auch wegen rotn € = Pe re) &d3 die Nor- 
malkomponente von rot & bestimmt. 


2) Nämlich wenn k? nicht Eigenwert des Systems AG + == mis = 0 längs 
ist. 
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deutige Lösung von (5b) nicht befriedigt wird, das System (1) überbestimmt. 


Weiter gilt: 


(A) Gibt man längs F nur zwei Komponenten von € (bezüglich eines cartesi- 
schen Systems x,y,2) und längs einer ebenen Randkurve auf F (welche als 


Schnitt von z = const. mit F entsteht) 
den ganzen Vektor € vor, so erhält 
man eine eindeutig bestimmte Lösung 
von (1)?). 

Obgleich der angeführte Satz, wel- 
chen wir später kurz beweisen, eine 
theoretisch befriedigende Lösung liefert, 
ist seine Anwendung für Probleme der 
Hochfrequenztechnik in dieser Form 
wertlos. Von größerem Interesse ist 
vielmehr die Frage, welche Bestim- 
mungsstücke von € längs einer An- 
regungsfläche AF eines elektromagne- 
tischen Hohlraumresonators (HR) frei 
wählbar sind, um zu einer eindeutigen 
Feldlösung innerhalb des HR zu ge- 
langen (siehe Abb. 1)). 

Einen ersten Hinweis zur Beant- 
wortung dieser Frage bekommt man 


F 


af 


Abb. 1. Schematische Ansicht eines elek- 

tromagnetischen Hohlraumresonators mit 

den Anregungsflächen AF, welche durch 

Gt + 0 charakterisiert sind. Auf den übri- 

gen Teilen der metallischen Berandung 
wird ©; = 0 angenommen 


aus der Reihendarstellung des Hohl- 

raumfeldes durch die vektoriellen Eigenfunktionen des Resonators. Bekanntlich 
lautet diese Entwicklung), wenn wir mit € das Feld der erzwungenen Schwin- 
gungen und mit ©,(v = 1, 2,...) die Eigenfunktionen von (1) (diese sind be- 
kanntlich durch [n, €] = 0 am Rande bestimmt) bezeichnen: 


¢ = Cy &, (x, %), (6a) 
wobei c, durch 
1 1 
Cy = [Gar &,) [Ga Pi In, &] Tot df (6b) 


gegeben ist. 
In (6b) tritt nur die Tangentialkomponente von €*) in bezug auf die Anre- 
gungsfläche AF auf und bestimmt in eindeutiger Weise das Hohlraumfeld. Es ist 


3) Siehe E. Ledinegg, Über Rand- und Sprungwertprobleme der Maxwellschen 
Gleichungen, S.-B. Akad. Wiss. Wien, Mathem.-naturw. Klasse, Abt. IIa, 156, H. 7 u. 8 
(1947). Wird im folgenden unter [4] zitiert. 

4) Bemerkenswerterweise hat auch diese Problemstellung nicht jene weitgehende 
praktische Bedeutung, welche man eigentlich erwarten sollte. Dies liegt vor allem daran, 
daß in Praxi Hohlraumresonatoren meistens durch Lecherleitungen oder Hohlrohre 
angeregt werden, welche, abgesehen von dem an der Anregungsstelle auftretendem Streu- 
feld, einen wohldefinierten Feldverlauf aufweisen. Die Berechnung des Streufeldes er- 
folgt näherungsweise, so daß die strenge Formulierung des „ersten Randwertproblemes“ 
auch hier mehr theoretisches Interesse besitzt. 

5) Siehe etwa E. Ledinegg u. P. Urban, Zur Theorie der erzwungenen Schwin- 
gungen elektrodynamischer Systeme. Acta Physica Austriaca 5, 4 (1952). 

%) Bekanntlich gilt: € = n (n, €) — [n [n, E]] = En + GE, d. h. also [n, €] liegt in 
der Tagentialebene und steht senkrecht auf &. Ferner ist |&| = |[n, €]. 
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also naheliegend anzunehmen, daß durch die Vorgabe von €; längs einer Anregungs- 
fläche AF und durch &; = 0 längs der übrigen metallisiert gedachten Begrenzungs- 
fläche F unser Randwertproblem eindeutig bestimmt wird. Diese, wie sich später 
ergeben wird, gerechtfertigte Vermutung läßt sich jedoch aus (6) nicht strenge 
herleiten, da wir weder über das Konvergenzverhalten der Reihe, noch über ihr 
Verhalten bei Annäherung von € nach AF allgemeine Aussagen machen können’). 
Es könnte also durch die getroffene Randbedingung möglicherweise eine Über- 
bestimmung unseres Problemes vorliegen. 

Wir werden im folgenden, von dem oben angeführten Satz (A) ausgehend, die 
Widerspruchsfreiheit der Randbedingung 


(B) {G, = 0, längs F, 
\G; + 0 beliebig vorgebbar längs F,, wobei F, und F, komplementäre 
Teilflächen von F darstellen, näher zu begründen haben. 


§ 2. Existenz des ersten Randwertproblemes von roetrE— KE=0 
und Konstruktion der Lösung 


Das System (3) lautet ausführlich geschrieben: 


AE, (x, 9,2) + E,(x, = 0 (7a) 
AE, (x, y.2) + k?E,(x, y, z) = 0 (7b) 
AE, (x, y, z) + k? E,(x, y, 2) = 0 (Te) 
ok, OH, , 0H, 


Nach dem Vorangegangenen sind die Komponenten E, = E,(x, y,2) und E, = 
E, (x, y, 2) aus (7a) und (7b) durch Vorgabe von E,= E, und E, = E}, wobei 
; und E* die Werte der betreffenden Komponenten längs F bedeuten sollen, bis 
auf eventuelle Eigenlösungen, welche zu den diskret auftretenden Eigenwerten 
+ gehören, eindeutig bestimmt. In (7d) ist damit auch die rechte Seite bekannt 
und ermöglicht daher die Berechnung von E,; wir finden als allgemeine Lösung: 


E,(x, y,2) = J p(x, y, t)dt + W(x, y). (8) 


Unter W ist eine willkürliche Funktion von x, y zu verstehen, welche noch geeignet 


zu bestimmen ist. Dazu gehen wir mit (8) in Gl. (7c) ein und erhalten nach kurzer 
Rechnung: 


[ (Ap + keg) dz + +e W =0. (9a) 
6 z= 
Da 
öE, , öE, 
Ap + 
wegen (7a) und (7b) und 


ba „ (r,8=1,2,3; % = 2) 


”) Falls die Reihe konvergiert, ist die Konvergenz sicherlich nicht gleichmäßig. 
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identisch verschwindet, bleibt als Bestimmungsgleichung für W: 


cp 
x 
A*W+kW (9b) 
Wird W am Rande (Schnittkurve von F mit z = 0) beliebig vorgegeben, so ist 
aus (9b) W eindeutig bestimmbar und damit ist auch E, (x, y, z) wegen (8) gegeben. 
Umgekehrt folgt aus (8) für z = 0 und f(x, y) = 0, wobei letztere Gleichung 
die Schnittkurve C von F mit z = 0 charakterisieren soll, die Beziehung: 


Ey (x, ¥,2= 0) = y). (9c) 


Das heißt: Gibt man E, = Ey, E, = E, längs F und außerdem E, = E3 längs C 
beliebig vor, so ist damit das System (7) oder das dazu äquivalente System (1) 
eindeutig bestimmt. Damit ist Satz (A) 
bewiesen. 

Wir versuchen nun im weiteren die 
Randwertbestimmung (A) durch jene von 
(B) zu ersetzen. Dazu denken wir uns eine 
beliebige Fläche F, gegeben, welche etwa bei 
z=( durch eine Ebene 2 abgeschlossen 
sei. F, + 2 bildet dann eine geschlossene 
Fläche (siehe Abb. 3). Längs F, geben wir 
wie bisher die Komponenten E, und E, 
d.h. also © = € vor. Langs F, wollen 
wir & = 0 fordern. Führt man in jedem 
Punkt von F, ein cartesisches Koordinaten- 
system (2%, %,, %) so ein, daß die x,-Achse 
mit der Flächennormalen n zusammenfällt, 
so besteht dann zwischen E,, E,, E, und Abb. 2. Die frei vorgebbaren Be- 
E,, E,, E, der folgende Zusammenhang: stimmungsstücke von € längs F 


| EY = 0 = + 04. + 0, (10a) 
[BS = 0 = 014, BS +0,,E (10b) 
E} = + £8. (10e) 


Die Gln. (10a) und (10b) verknüpfen also die Randbedingung mit der Randwert- 
bestimmung (A). Daraus erhält man 


E,=pl(&n)4,} (1) oder auch A ‚ (12) 
ES = p(&,n) Ay = 


wobei A,, A,, A, drei bestimmte Flächenfunktionen bedeuten, während unter 
p(£, 7) in (11) eine willkürliche Funktion der Flächenkoordinaten £, 7 zu verstehen 
ist. 

Da einerseits (12) mit der Randwertbestimmung © = 0 (längs F,) äqui- 
valent ist, andererseits die Bedingungen (12) weder mit der Forderung einer will- 
kürlichen Vorgabe von E, und E, längs F, gemäß (A) in Widerspruch stehen, 
noch darüber hinausgehen, dürfen wir annehmen, daß unsere spezielle Randwert- 
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bestimmung (B) (E&; = 0 längs F,, & = €} längs &) unser Problem nicht über- 
bestimmt ®). 

Die Forderung der Vorgabe von € längs der Randkurve C (z = 0) ist hier 
automatisch erfüllt, da aus Stetigkeitsgründen in allen Punkten der Randkurve, 
©, die Normalprojektion von € auf & darstellt. 

Wir können nunmehr zeigen, daß (B) zur eindeutigen Determinierung des 
Lösungsfeldes gerade ausreicht, daß also auch keine Unterbestimmung vorliegt. 
Bekanntlich besitzt das homogene 
System 

rot retE — 


mit &=0 längs der gesamten 
Hüllfläche F (in unserem speziellen 
Fall längs F und 2) ein unend- 
liches diskretes reelles Eigenwert- 
spektrum k, mit den dazugehörigen 
Eigenfunktionen €,?). 
Wären also zwei Felder &) und 
Abb. 3. Schematische Skizze (Aufriß) des pa ee 
unter (B) formulierten. Randwertproblemes +k, annehmen wollen, so ist 
mit einer ebenen Abschlußfläche = dann das daraus gebildete Differenz- 
feld EU — & Eigenfunktion 
von (1) mit [CH =0 am Rande zum Eigenwert k? + k; d.h. also & muß iden- 
tisch verschwinden !P). 

Wir haben noch den Übergang zur allgemeinen Randwertbestimmung gemäß 
(B) auszuführen und nehmen dazu eine geschlossene Fläche F an. Diese werde 
durch eine geschlossene Flächenkurve C in zwei Gebiete F, und F, zerlegt. Inner- 
halb F, wählen wir & = 0 als Randbedingung. In F, konstruieren wir uns ein 
beliebig zu verfeinerndes Koordinatennetz und greifen ein ebenes Flächenelement 
AF, , heraus und geben dort & = Er” iF vor (Abb. 4). Im übrigen sei in F, 
ebenfalls &; = 0. Die sicher vorhandene eindeutige und reguläre Feldlösung 
wollen wir mit © = E»%) bezeichnen. Dann ist wegen der Linearität unseres 
Problemes aber auch 


— & 
®) Man könnte auch folgendermaßen schließen: Die oben angeführte Randwert- 

bestimmung ist für die Diff.-Gl. AE + k* € = 0 unterbestimmt, da längs F, und z = 0 
zur eindeutigen Bestimmung des Feldes der ganze Vektor © vorzugeben wäre. Die durch 
(11) in das Lösungsfeld eintretende willkürliche Funktion wird es daher im allgemeinen 
gestatten, noch die Divergenzbedingung div € = 0 zu erfüllen. Wir bemerken jedoch aus- 
drücklich, daß dieser sehr plausible Sachverhalt noch durch einen strengen Beweis zu 
ergänzen wäre. Z. B. könnte man an (12) anknüpfend ein Iterarionsverfahren einleiten. 
Von einer beliebigen vorgegebenen Funktion E}(£,») ausgehend, welche man an Stelle 
von E} einführt, gewinnt man gemäß (A) ein Feld "E,, "E,,”E,. Führt man nun ‘EY 
in (12) ein, so erhält man auf analoge Weise ein Feld ”E,, ”E,, ”E,. Durch Fortsetzung 
dieses Verfahrens entsteht eine Folge, mit deren Konvergenz oder Existenz eines Häu- 
fungspunktes der strenge Beweis erbracht ist. 

®) Der diesbezügliche Beweis ist z.B. für allgemeine zylindrische Begrenzungs- 
flächen in: E. Ledinegg, Ann. Physik 41, 7 (1942). 

10) Die ausreichende Bestimmung eines Lösungsfeldes von (1) gemäß der Rand- 
bedingung (B) geht auch aus Gl. (20) hervor. 
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eine Lösung von 
rot rot E— k? € = 0, 


wobei längs F,, €; = 0 ausfällt und im Bereich F, in jedem Flachenelement AF, ,, E; 


nach unserer Konstruktion den Wert &”"* 
annimmt. Durch den Übergang AF, ,—0 
erhalten wir dann das gewiinschte Ergebnis: 


Das elektromagnetische Feld {€, 5}, wobei 
€ der Diff.-Gl. (1) genügt, wird im Inneren 
einer geschlossenen teilweise metallisierten 
Fläche F (€ = 0, längs F,) in eindeutiger Weise 
durch Vorgabe von © = € + 0 an den nicht- 
metallisierten Teilen (F,) der Fläche bestimmt. 


Der noch interessierende Übergang F, — 0, 
welcher dann zum ersten Randwertproblem 
in koordinatenfreier Form herüberleitet, soll - 

im letzten Abschnitt besprochen werden. Im = bal 
folgenden geben wir noch die Lésung von (B) ll Wh beliebig gestalteter 
mittels einer Integralgleichung an. Anregungsfläche 


$3. Zurückführung der Randwertbestimmung (B) auf die Auflösung einer 
Integralgleichung der „Flächenstromdichte“ 


Analog zur skalaren Potentialtheorie läßt sich das elektromagnetische Feld 
in einem abgeschlossenen Raumteil durch seine Randwerte und den Ableitungen 
von & am Rand explicite ausdrücken. Die diesbezügliche Integraldarstellung 
ist schon seit langem bekannt!) und lautet für das &-Feld: 

— innerhalb V 
$ {[n, rot E] + [In, E] Ve] + (n, ©) Vp}df = | — 27€ am Rande F (13) 
0 außerhalb V. 
Wir beweisen im Nachstehenden kurz Gl. (13), indem wir zum Unterschied der 
in Ya), 1b) und !!e) geführten Ableitungen rein formal vorgehen. Wir benützen dazu 
die leicht zu beweisende Identität 


rot {p rot E — [Vp, E]} = prot rot E + EAP— 2 (EV) Ve + V(Vy,&) (14) 


und integrieren über das betrachtete Volumen V, welches von der geschlossenen 
Fläche F begrenzt wird. Wir erhalten: 


f rot {p rot E— [Vp, E]} dr = je rot rot © + © Ag) dr 
v 


+ ©) — 2 Vp, €) Vo} dr. (15a) 


Wegen 
[rot {g rot El}dr = $fIn, rot [m (Mp, El} af 


und 


= © df, 


11) a) Lamor, Lond. Math. Society Proc. (2) 1 (1903); b) Kottler, Ann. Physik 
[4] 70, 405 (1923); ec) Bucholz, Arch. Elektrotechn. 37, 22 (1943). 
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sowie 

(EV) = n) Vo df, 

v 
wobei hier V& = 0 vorausgesetzt wird, geht (15) zunächst in 

je rot rot € + EAp)dr = $ {In, rot E] p— [n [Vp, E]]} df 
über. Führt man noch in der obigen Gleichung die Identität 
[Ve In, €]] — In (Vp, E]] = n (Ve, &)— Vo (n, ©) 


ein, so erhält man 
je rot rot E— E Ag) dr = $ {In, rot E] + [[n, €] Ve] + (n, ©) Ve} df. (15h) 


Die bisher unbestimmt angesetzte skalare Funktion g wählen wir nunmehr als 
ein partikuläres Integral von Ag + k? gp = 0 in der Form e/**/r und erhalten auf 
der linken Seite von (15b) identisch Null, während auf der rechten Seite noch ein 
Oberflächenintegral $ {-- -}df, welches die singuläre Stelle r = 0 ausschließt, 
hinzutritt. Fugel 

Der Grenzübergang 

lim (n,€) Vodf 

r>0 Fgugel 
liefert dann unmittelbar — 42 €, bzw. — 22 € bzw. 0 je nachdem der gewählte 
Aufpunkt innerhalb V oder auf der Oberfläche F bzw. außerhalb V gelegen ist. 

Jedes E-Feld läßt sich also innerhalb V bzw. auf F in der Form (13) darstellen. 
Es gilt auch im gewissen Sinne die Umkehrung: 

Ist eine reguläre Lösung © von (1) in einem Gebiete @ bekannt, welches 
die Berandung F ganz enthält, so liefert G1. (13) in eindeutiger Weise die reguläre 
Lösung € von (1) im gesamten von F eingeschlossenen Volumen V; dabei fällt € 
innerhalb G mit € zusammen. 

Im übrigen besitzt Gl. (13) zur Lösung der ersten Randwertaufgabe in dieser 
Form nur theoretisches Interesse, da im Integranden neben [n, €] auch noch 
[n, rot €] und (n, €) vorkommen 

Für das Folgende ist noch die Darstellung von rot € von Bedeutung. Indem 
wir zwischen Aufpunktdifferentiationen (Index P) und Quellpunktdifferentia- 
tionen (Index Q) unterscheiden und zur Abkürzung 


In, rot &] = j® (n, ©) = q® 

[n, E] = 5" (n, rot €) = q™ 
12) Wir bemerken, daß Gl. (13) durch Zusammenziehung von F über eine Unstetig- 
keitsfläche der Maxwellschen Gln. zu einer expliziten Lösung des Sprungwertproblems 


führt, welches andererseits einen Ausgangspunkt der Beugungstheorie bildet. Zum 
Unterschied des Randwertintegrales (13) können im Sprungwertintegral 


—4a€ = | + div* j” $ n,) Vp ds, 


welches außerhalb F eine überall reguläre Lösung von (1) darstellt, die Sprungwerte 


j) = [n, €], j® = [n, rot €], wobei E und rot & die Diskontinuitäten der Feldstärke 
oberhalb und unterhalb von F charakterisieren, beliebig vorgegeben werden. Siehe dies- 
bezüglich !!e) und [4]. 
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setzen, erhalten wir nach kurzer Rechnung 
+ + Var} df = Veyl— j App 
+ 5% Ve) Voy} af. 


Apg+kop=0 
und der aus dem Stokeschen Satze abzuleitender Integralformel 


Ve) Vow af = {We A] n) Voy af + ds) Veg, 


wobei X ein beliebiges genügend oft differenzierbares Vektorfeld charakterisiert, 
bekommt man unmittelbar: 
—4z rot & 
rot & 
Gl. (16) läßt sich nun zur Auflösung der Diff.-Gl. (1) mit der Randwertbestim- 
mung (B) und auch zur Lösung des homogenen Problemes (€; = 0 am Rande) 


verwenden. Beginnen wir mit letzterem; da längs F, [n, €] und (n, rot &)'*) ver- 
schwinden muß, bleibt hier 


— 2rzrotp& = $ [[n, rot E] Va q] df. (17a) 


Wegen 


$ {[{n, rot E] Ve] + In, E]p + (rot E,n) Ve} df. (16) 


Gl. (17) läßt sich unmittelbar in eine homogene Integralgleichung zur Bestimmung 
von j® = [n, rot €] umschreiben. Unter Beachtung von 


[np [np, rotp E]] = —rotp € 
22 [np, 7] = $ df. 
Berechnet man aus (17b) 7 15), so findet man dann aus 
—4n6= $ [j, df 
F 


das zugehörige &-Feld innerhalb V. 
Analog hat man im Falle (B) vorzugehen. Im Gebiete F, ist wie früher 
[n, €] = 7 = O und (n, rot &) = q® = 0, 
während innerhalb F, [n, €] = 7 + 0 vorzugeben ist. 
Mit der Kenntnis von j") ist auch q als bekannt anzusehen, denn es gilt: 


1 
— (n, rot &) = — lim $ (ds, n, j"). 
)=— lim ™) 


erhält man nämlich 


Man kann auch g® als Divergenzausdruck darstellen. Mittels des Gaußschen 
Satzes (für ein zweidimensionales Kontinuum) erhält man: 


1 
div 7) = lim — ds, n, j"). 


13) Man hätte (16) auch aus rot rot (rot ©) — n? rot € = 0 mit direkter Anwendung 
von Gl. (13) erhalten können. 
14) Man beachte, daß längs F, $ stets senkrecht auf © stehen muß. 
15) Es sind hier nur für k? = k} reelle singularitätenfreie Lösungen vorhanden. 
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Wir finden also 
= — divx (Is 
Aus (16) folgt dann: 


$ {ke div* vr} df 2n In, 72] | Ve) df. (19a) 


Man bekommt zur Bestimmung von 57 eine lineare, nicht homogene Integral- 
gleichung, deren linke Seite durch die Vorgabe von [n, €] = 7 wegen (18) völlig 
bestimmt erscheint. 

Das &-Feld innerhalb V findet man wieder aus (16); es gilt: 


= $ Vp] df + $ {k2 — divX Vo} df. 


Wir können also die Integration von (1) mit der Randwertbestimmung (2) auf 
die Integration einer zweidimensionalen vektoriellen Integralgleichung vom 
Typus (19a) zuriickfiihren. Fiihrt man noch ein orthogonales Koordinatensystem 
mit den normierten Grundvektoren e,, e, auf F ein und setzt zur Abkürzung 


so erhält man an Stelle von (19) das folgende System: 
= — + $ {his (P,Q) 7 + ky, (P,Q) df 
— then (P,Q) 9? + (P,Q) 72} af. 


Dabei bedeuten 7, 5 die Komponenten von j® in bezug auf das zugrunde 
gelegte Koordinatensystem. Die Komponenten des Greenschen Tensors von (20) 


sind durch k,,= (ee, Agg) gegeben. Es gelten die Symmetrierelationen 
kr,r (P,Q) = kyr (Q, P), kr,s (P,Q) = La (Q, P). 


(20) 


§ 4. Das erste Randwertproblem in koordinatenfreier Form 


Wir haben schon am Ende von $2 bemerkt, daß durch den Grenzübergang 
F,—0 die Randwertbestimmung (B) in die erste Randwertaufgabe übergeführt 
wird. Ein unendlich kleines Gebiet in F bleibt aber notwendig weiter metalli- 
siert, d.h. außer der Existenz mindestens eines Punktes P auf F, in welchem 
©, = 0 ausfällt, muß noch die Limesrelation 


im =0 (21a) 
4r+0 Ar 
gelten. Mit Ar haben wird den Abstand von P zu einem benachbarten Flächen- 
punkt bezeichnet. Wie man leicht einsieht, ist (21a) auch mit 


m $rot„Edf=0 (21b) 


li 
AF>0S 


äquivalent. Damit erhält das erste Randwertproblem von (1) zunächst die fol- 
gende Gestalt: 

„Wählt man &} derart, daß mindestens in einem Punkt P von F & = 0 aus- 
fällt und gilt außerdem (21a) ober (21b), so ist bei sonst beliebiger Vorgabe &; =& 
längs F, die Lösung von (1) innerhalb V eindeutig und regulär. Für k® sind dabei 
die Eigenwerte k, des homogenen Problemes auszuschließen.‘ 
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Von den, die freie Vorgabe von & belastenden Zusatzbedingungen: &{ = 0 
in P und Gl. (21a) bzw. (21b), können wir uns befreien, wenn man ausdrücklich 
die Stetigkeit von @{ und die Existenz sowie die Stetigkeit der Flächenableitung 
in jedem Punkte von F fordert. Es gilt nämlich der folgende Satz'*): (/) „Ist auf 
der Randfläche F eines einfach zusammenhängenden Bereiches ein stetiges 
Flächenvektorfeld 5 gegeben, so gibt es mindestens eine Stelle P auf F, in welcher j 
verschwindet.‘ 

Ohne auf mengentheoretische Feinheiten einzugehen, sei ein kurzer Beweis 
skizziert. Prinzipiell hat man nur die Stetigkeit des Feldes vorauszusetzen. Wir 
greifen einen Flächenpunkt P, heraus und legen durch diesen eine Vektorlinie Lo. 


Abb. 5. Die Skizze veranschau- 


licht, daß die durch L, und L, ande 


müssen durch Fortsetzung des 


begrenzte Fläche F, nicht be- 
liebig klein ausfallen kann, da 
sonst das Vektorfeld in der 
Umgebung des von P, aus- 
gehenden Kurvenzuges nicht 
stetig wäre (entgegengesetzte 
Richtung der Flächenvektoren) 


Einschachtelungsverfahrens 
sich um einen Punkt P oder 
auf ein Kurvenstück C zu- 
sammenziehen. Jede e-Um- 
gebung von P oder eines 
Punktes von C enthält daher 
mindestens zwei Vektoren ent- 


gegengesetzter Richtung 


Verschwindet 7 in keinem Punkt von Z,, so muß die Vektorlinie jedesfalls geschlos- 
sen sein. Wäre nämlich ein Anfangs- bzw. Endpunkt P, auf Z, vorhanden, so 
widerspräche dies der geforderten Stetigkeit, da nach Voraussetzung längs Z,, 
j + 0 ausfällt. Ist ferner L, stückweise rektifizierbar, so muß die Gesamtlänge der 
Vektorlinie endlich sein. Dies folgt ebenfalls aus der angenommenen Stetigkeit 
durch Ausführung eines dem Bolzano-Weierstraßschen Satze über den Häu- 
fungspunkt einer ebenen Punktmenge nachgebildeten Verfahren. 

L, teilt F in zwei Gebiete, F, und F,. Wir greifen einen Punkt P, in F, heraus 
und legen wiederum eine Vektorlinie LZ, durch, P,. Diese darf aus Stetigkeits- 
gründen Z, nicht kreuzen, verläuft also ganz in F,. Mit der Voraussetzung 7 + 0 
folgt wieder die Geschlossenheit von L,. Das Gebiet F, zerfällt dennoch durch L, 


18) L. E. J. Brouwer, Proc. Roy. Acad. (Amsterdam) 11, 850 (1909). Leider war 
uns die Brouwersche Arbeit nicht zugänglich, sondern wir entnehmen den angeführten 
Satz aus einer kurzen Notiz von: H. F. Mathis, J. appl. Physics. 
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in F, und F,. Wir betrachten nun etwa das Gebiet F, und stellen zunächst fest, 
daß der Flächeninhalt von F, sicherlich um einen endlichen Betrag kleiner als der 
von F, ausfällt. Dies folgt unmittelbar aus der Stetigkeit des Feldes, wie anschau- 
lich aus Abb. (5) zu entnehmen ist. 

Durch Fortsetzung des angedeuteten Schachtelverfahrens kommt man schließ- 
lich zu den beiden in Abb. 6 angedeuteten Möglichkeiten; entweder ziehen sich 
die Vektorlinien ZL; um einen Punkt P oder um ein Kurvenstück C zusammen. 
Damit ist in beiden Fällen ein Widerspruch zur Stetigkeitsannahme konstruiert !”), 
Es muß also im Sinne unseres Satzes mindestens eine Nullstelle von j vorhanden 
sein. 

Das angeführte Theorem gestattet eine Verschärfung, wenn nicht nur die 
Stetigkeit des Feldes j,sondern auch die Existenz und Stetigkeit der ,,Flichen- 
ableitung‘‘ gefordert wird, und in diesem Falle gilt: 

II) „Ist j ein stetig differenzierbares Flächenvektorenfeld auf F, so verschwin- 
det an den Nullstellen von 7 auch dessen Flächenableitung. Dabei wollen wir 
unter den Flächenableitungen einen kovarianten Prozeß im Riemannschen 
Sinne verstehen. Für die Nullstellen von j folgt daraus auch das Bestehen der 
Limesbeziehung (21a)‘‘!®). 

Identifizieren wir nun unter j das Tangentialfeld von €, so liefern die Sätze I) 
und II) offenbar die folgende Formulierung des ersten Randwertproblems: 

Die Maxwellschen Gleichungen bzw. die Diff.-Gl. (1) 


rot rot E— = 0 


liefern unter der Voraussetzung k? + k,, wenn am Rande F eines abgeschlossenen 
Bereiches V 

als stetig differenzierbares Feld, in sonst beliebiger Weise vorgegeben wird, eine 
eindeutige Lösung € = E(x, x, x) mit stetiger zweiter Ableitung. Die Tangen- 
tialkomponente von € fällt dabei längs F mit €{ zusammen. 

17) Ziehen sich die Kurven Z, etwa auf P zusammen, so liegen in einer e-Umgebung 
von P zwei verschieden gerichtete Vektoren von endlichem Absolutbetrag. 

18) Der Beweis des Satzes ergibt sich aus der Betrachtung eines etwa isoliert ange- 
nommenen P-Punktes (j = 0) des Feldes. Die Flächenableitungen von j müssen dann 
in einer e-Umgebung entgegengesetzt gerichtet sein, woraus sich wie früher aus dem 
Widerspruch zur Stetigkeitsannahme unsere Behauptung ergibt. 


Graz, Institut für theoretische Physik der Universität. 


(Bei der Redaktion eingegangen am 8. Februar 1952.) 
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Theorie der Beugung am Zylinder unter Berücksichtigung 
der Kriechwelle 


Von W. Franz und K. Deppermann 
(Mit 3 Abbildungen) 


Inhaltsiibersicht 


Messungen der Beugung von cm-Wellen an Metallzylindern von A. Limbach 
zeigen eigenartige Maxima und Minima. Diese werden als Interferenzen zwischen 
der geometrisch reflektierten Welle und einer um die Riickseite des Zylinders 
herumgekrochenen Welle gedeutet, und mittels der Integralgleichungen der 
Beugungstheorie von A. W. Maue quantitativ berechnet. Die Ubereinstimmung 
mit dem Experiment ist gut. 


Einleitung 


Die Beugung akustischer und elektromagnetischer Wellen an einem ideal 
reflektierenden Zylinder vom Radius a behandelt man nach dem Verfahren von 
Mie durch einen Reihenansatz, welcher nur für große Wellenlängen (A sehr groß 
gegen «) gut konvergiert, für kleinere Wellen dagegen sehr mühsam auszuwerten 
ist. Für sehr kleine A/a stimmen die Ergebnisse der Kirchhoffschen Beugungs- 
theorie bzw. der geometrischen Optik mit dem Experiment gut überein, dagegen 
zeigen Messungen im mittleren Gebiet (A wenig kleiner als a) Eigentümlichkeiten, 
welche durch die asymptotischen Formeln nicht wiedergegebem werden und aus 
den Entwicklungen von Mie schwer ersichtlich sind. Es ist daher sinnvoll, für 
dieses Gebiet nach einer Näherungslösung Umschau zu halten, die die wahre 
Lösung bis auf einen für die Praxis erträglichen Fehler approximiert und die 
rechnerisch leichter zu handhaben ist als die Reihen. Dieses soll im folgenden mit 
Hilfe der Integralgleichungen der Beugungstheorie, die wir von A. W. Maue') 
übernehmen, durchgeführt werden. Für eine skalare Welle u lauten diese für ein 
Beugungshindernis mit beliebiger Oberfläche 


1 /ea(P, 
F 
für die Randbedingung a = 0 und 


Gu(P) _ Oue(P) 1 2G(P,Q) eu(Q) 
| én, ang “la (2) 


für die Randbedingung u = 0. 
1) A. W. Maue, Z. Physik 126, 606 (1949). 
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Der Aufpunkt P und der veränderliche Integrationspunkt Q liegen beide in 
der Oberfläche des Beugungshindernisses. ö/ön ist die Ableitung nach der vom 
Hindernis nach außen weisenden Normalenrichtung, u, bedeutet die Primär- 
erregung; 


G(P,Q =" 


(mit @ = Abstand zwischen P und Q) ist die Greensche Funktion des un-, 
endlichen Raumes; dfg ist das Flächenelement im Punkt Q; zu integrieren ist über 
die gesamte Oberfläche des Hindernisses. Für ka>1 (k= 2z/A) ist die 
Kirchhoffsche Beugungstheorie anwendbar. Sie vernachlässigt jedoch Glieder, 
welche zwar für ka—oo verschwinden, aber im Bereich mittlerer ka eine aus- 
schlaggebende Rolle spielen, und ihre physikalische Ursache darin haben, daß die 
Erregung sich über die Ränder der beleuchteten Fläche hinweg ausbreitet und — 


wenn auch infolge Abstrahlung exponentiell gedämpft — auf der Rückseite des . 


Zylinders weiterkriecht. Dieser Vorgang kann mittels der Integralgleichungen 
verfolgt werden, indem die Integrale nach der Methode der stationären Phase 
näherungsweise ausgewertet werden. Vorauszusetzen ist dabei, daß zwar ka, 
nicht aber (ka/24)* groß gegen 1 ist. 


$1. Lösungsansatz für Randbedingung =0 


Wir führen ein kartesisches Koordinatensystem x, y, z und dazu ein Zylinder- 
koordinatensystem r, y, z ein. Die Koordinaten des Aufpunktes werden mit 
großen, die des Integrationspunktes mit den entsprechenden-kleinen Buchstaben 


_oberer Schattenrand 
Lichtgebiet O<@sır 
Einfallsrichtung 
Schattengebiet p 
n<®<2n 
J unterer Schattenrand 
xt 


Abb. 1. Bezeichnungen 


bezeichnet. Die Primärerregung, die eine in der y-Richtung einfallende Planwelle 
(s. Abb. 1) auf der Zylinderoberflache hervorruft, kénnen wir dann schreiben: 
u,(P) = Ay e-ikasin® (4) 
Für die Lösungsfunktion u der Integralgleichung (1) machen wir den Produkt- 
ansatz aus einer langsam veränderlichen Amplitudenfunktion mit einer rascher 
veränderlichen Phasenfunktion. Dabei setzen wir in der beleuchteten Zylinder- 
hälfte plausiblerweise als Phase der Lösungsfunktion die der eingestrahlten Welle 
an. Dies ergibt für w: 
u(P) = A(®) e-ikasin® o<®<a. (5) 
Für den Ansatz im Schattenbereich lassen wir uns von der Vorstellung leiten, 
daß vom oberen und unteren Schattenrand her je eine Welle mit exponentiell 
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gedämpfter Amplitude auf der Oberfläche des Zylinders in das Schattengebiet 
hineinkriecht und dann unendlich oft den Zylinder umläuft. Wir nennen dabei 
den Exponentialfaktor, mit dem sich die Amplitude beim Fortschreiten um den 
Winkel ® multipliziert, A(®), mit der Funktionalgleichung 

A(® + 4) = A@) - Aly). (6) 
Für die vom oberen Schattenrand herkommende Kriechwelle erhalten wir 
7 
Entsprechend ergibt sich für die vom unteren Schattenrand ausgehende Kriech- 
welle 


u,(P) =C 


ika(2n-®) vt 

(8) 
—A(2a)e 

Darin bedeutet C eine später zu bestimmende Konstante, die aus Symmetrie- 
gründen für %, und u, dieselbe ist. Als Gesamterregung eines Aufpunktes im 
Schattengebiet erhalten wir mithin: 

1 — A(2 2) 
Wir werden später sehen, daß die durch (9) dargestellte Kriechwelle die zu (1) 
gehörende homogene Integralgleichung bei geeigneter Verfügung über die Dämp- 
fung der Amplitude löst, also für jeden Oberflächenbereich des Zylinders gilt. 
Die Integralgleichung (1) lautet in Zylinderkoordianten 


u,(P)=C 


2x +00 a (eke 
=2(P) + dp f (10) 

Da u(Q) wegen der Zylindersymmetrie des Problems von z nicht abhängt, 
kann die Integration iiber z an der Greenschen Funktion ausgefiihrt werden mit 


dem Ergebnis 


ein? | u(Q) HY (2 ka 


sin? ) ao. (11) 


§ 2. Lésung der Integralgleichung im Lichtgebiet 


Die Integralgleichung (11) wird mit dem Ansatz (5) fiir Aufpunkte im Licht- 
gebiet 


4) [ |sin® Aw H® (2 ka dp. (12) 
0 


Die Kriechwelle ist, hierin unterdrückt, da sie, wie oben bemerkt, als Lösung der 
homogenen Integralgleichung angesetzt ist und daher nachträglich zu einer par- 
tikulären Lösung von (12) hinzugefügt werden kann. Freilich gilt dies nicht un- 
mittelbar an der Schattengrenze, doch erscheint es im Hinblick auf die Ver- 
hältnisse im Schatten (s. $3b) nicht unplausibel, nur die Zustrahlung von der 
Lichtseite in Rechnung zu setzen (auf die Erstellung einer im ganzen Über- 
gangsgebiet brauchbaren Lösung wollen wir demnächst zurückkommen). Gl. (12) 
ist eine Integralgleichung für die Amplitude, die wir durch ein Iterations- 
verfahren lösen. Als Ausgangsnäherung wählen wir die Inhomogenität A(®) = 2 A, 
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(dies entspricht dem Kirchhoffschen Näherungsansatz). Durch sukzessives 
Einsetzen erhalten wir als vollständige im ganzen Lichtgebiet gültige Lösung: 
n-® 


[ sin (o+= +). H® (2kasin 


-® 
sin 


-i2 kasin_ cos (o+2-+ r) 
2 2 


A(@®) = 24, 


sin HW (2%. sin dr do 


=; in?! _itkasin cos (0+) (13) 


— cos (e+ ++) +sinZ cos 
2 2 2 2 2 


AD (2ka sin (2kasin (2ka sin!2!) dr do dy + ore 


= 24, {1 + 19 + 4+ 1M +... +}, 
worin die Form des allgemeinen Gliedes ohne weiteres ersichtlich ist. In der Mitte 
der beleucheten Fläche genügt es, die ersten beiden Glieder zu berücksichtigen. 
Wir können für große k a das Integral nach der Methode der stationären Phase 
näherungsweise auswerten mit dem Ergebnis: 


7 

mit 
1 

sin? ' 
In der Umgebung der Schattengrenze (® = 0 oder 2) ist Gl. (14) nicht brauchbar. 
Für ® = 0 (x) können wir jedoch (13) näherungsweise auswerten, was für den 
Anschluß der Lösung des Schattengebietes an die des Lichtgebietes wesentlich ist. 
Wir führen die Rechnung für den unteren Schattenrand ® = 0 durch; wegen der 
Symmetrieverhältnisse ergibt sich für den oberen Schattenrand dasselbe. Durch 
asymptotische Entwicklung der Hankelfunktion und Anwendung der Methode 
der stationären Phase auf die Integrale erhalten wir im einzelnen: 


In — Lane (16) 


Durch die Substitution 


8\t. 173 
[6 4 = -— 16 4 
(14) dt 3 164-4 dt 


und Verlegung der Integration über ¢ auf die reelle Achse ergibt sich 


F@) = f 2 cos (zcos®) H( (z) dz = (15) 
0 


1 


3° (17) 
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Für /® ergibt sich: 


Mit der Substitution 
8 


81% . A 
und Verlegung der Integration in das Gebiet reeller 2 und u folgt: 
1 oo +1 su. 1 1 
J] dd J] duYl— wie 2)= 0,0385 158. (19 
iad J (19) 
Für /® ergibt sich: 
70) = ein (ka)! 
= 82} Y 


e 


dt do dy. 
Substitution : 
8 
(20) 
"Ver +484 +5 — 
Durch die 
u = sinz 
erhalten wir: 
Vv (1 + cos 2 2) 
I® = — —- d d 
+ 1608 + + 9 + cos 22)! 
TG) Vr(1+ cos 2 z) 
=~ ju fa ("+ + 8» + 
2 
1 Vv 
terrain” 
und analog: co 


Ve 
Ie) <—— - dv. 21b 
as | + 89+ 
ve 


Die Integrale (21a) und (21 b) werden numerisch ausgewertet. Dazu ist es zweck- 
mäßig, die Substitution » = einzuführen. Es ergibt sich damit: 


Y (22a) 
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bzw. 
oo 


5 
y2 J (y? + 4/2 


Die numerische Berechnung dieser Integrale ergibt: 


19 > — (0,01 104 + 0,00 006) a 


19 < — (0,00 974 -L 0,00 006) a 
1® = — 0,00 735 + 0,00 050. (23) 


Da weiter 
IM< | J®|2 = 0,00 124 


ist, und da (14) für ® = 0 eine alternierende Reihe ist, folgt: 
A(0) = A, (1,390 + 0,002). (24) 
Die durch die Vernachlässigung weiterer Reihenglieder entstehende Fehler hat 
die Größenordnung (ca) ‚so daß wir insgesamt erhalten: 
A(0) = A, [1.390 + 0,02 + o( a )| (25) 
(ka)! 
$3. Lösung der Integralgleichung im Schattengebiet 


Die Integralgleichung (11) wird mit den Ansätzen (5) und (9) für Aufpunkte 
im Schattengebiet, wenn wir 4) durch die ersten beiden Glieder seiner asym- 
ptotischen Entwicklung ersetzen: 
é 4 A(2 — @) ®) 
1 


ini) af | sin? 5? Ay -ika(sing- zsin ( 4 37 


16 kasin 


=) 4, e-ikasin® 


27 


ka [| _m) ei pika(22-9) (26) 
[ sin? | AY a) e +4(22 —p)e 
4a 2 1 — A(2 2) ef 
i2ka sin 3% 
| 16ka sin” 


Durch Nullsetzen der ersten Ableitung des Exponenten nach @ überzeugt 
man sich leicht davon, daß das erste Integral, welches über das Lichtgebiet er- 
streckt wird, genau eine Stelle stationärer Phase besitzt, und zwar beig, = 22 —9, 
also an derjenigen Stelle der beleuchteten Fläche, deren Schatten auf den Auf- 
punkt fällt; es ist daher physikalisch einleuchtend, daß der von g, herrührende 
Beitrag zur Amplitude (wie sich im folgenden zeigen wird) die Primäramplitude 
näherungsweise auslöscht (Abschattung). Die Phase des zweiten Integranden 
wird nur in der Umgebung des Aufpunktes stationär und bestimmt das Fort- 
kriechen der Erregung längs der Oberfläche. 
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a) Berechnung der Abschattung 
Nach der Methode der stationären Phase läßt sich der zweite Summand der 
rechten Seite für große k a näherungsweise ersetzen durch: 


+00 ika 
Bi T 1 
ikart /7 (k a)? cos? 


Die Auswertung ergibt bei Beschränkung auf Glieder der Größenordnung 1/k a 
~ t 
(28) 
Da der Aufpunkt entfernt von den Schattenrandern liegen sollte, erhalten wir unter 
Berücksichtigung von (15) angenähert ; 
=... 2A, h- eikasin® (29) 


was in der Tat zur Abschattung (Beseitigung der Inhomogenität) in (26) führt. 
Damit bleibt die homogene Integralgleichung: 


Ad — a) eika(®-m) + —®) eika (22-9) 


— (3iala VE sin? ‘(A@— 7) 4 A(2n — eika(ax-9)) 
8 30 
fi dq. (30) 
16 ka'sin® 


b) Lösung der homogenen Integralgleichung 
Wegen der Symmetrieverhältnisse genügt es, statt der Integralgleichung (30) 
die Integralgleichung 


ika(p+2 sin 


9-® 


zu lösen, welche die Amplitude der vom unteren Schattenrand herkommenden 
Kriechwelle bestimmt. Von dem Integral haben wir, wie oben bemerkt, nur Um- 
gebungsbeiträge zu erwarten, und zwar zeigt die Entwicklung des Exponenten 
um ®, daß die Phase nur im Bereich  < ® stationär wird. Physikalisch bedeutet 
das, daß die Umgebungszustrahlung von der Lichtseite herkommt. Sie ergibt 
sich aus (31) näherungsweise zu: 


= 24 D 
A(®) = Ves J A@—2))re dt. (32) 
Die Integralgleichung läßt sich mit dem Ansatz 
/6 3, 
Ald)=e (=) (33) 
befriedigen. Für x ergibt sich die tränszendente Bestimmungsgleichung 
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Mit der Substitution 
24 

— 
T=ei 


1= iy2 f tt dt. (34) 
% 


Von den unendlich vielen Lösungen für x kommen nur solche in Frage, für die 
a ei*/® einen positiven Realteil hat, andernfalls würde ja, wie aus (33) hervorgeht, 
die Amplitude im Schattengebiet nicht gedämpft sein. Weiter dürfen wir bei der 
Ausbreitung auf der schwach gekrümmten stetigen Zylinderfläche einen relativ 
glatten Verlauf erwarten, weshalb große Beträge von x sich verbieten. Die beiden 
Lösungen von (34) mit den kleinsten Beträgen sind: 
= 2,33 246 e-i”/® 
= 5,35 295 ei 2/3, 


Die weiteren Lösungen lassen sich mit großer Genauigkeit aus der folgenden 
asymptotischen Formel entnehmen: 


3(4 ) (2n—1) 
6(—1)"—1/4 1 5 15(—1)"—1/4\4 1 

108 G 729 (=) 
ei=/®.x, ist rein imaginär, alle weiteren x-Werte haben Beträge größer als 7. Damit 
erweist sich x, allein als brauchbar, und (33) geht über in 
~2,33246 
A(®)=e (3) (36) 
Bei genauerer Rechnung ergibt sich an Stelle von (34) die folgende Bestimmungs- 
gleichung für «x: 


erhält man daraus 


(35a) 


3/8 € € 
= + t-# —— — 
J (1 dt (37) 
(Fe 


Für «, folgt hieraus 
0,419 


(kat (38) 


= 2,33 246 7/3 + 
c) Rest der Primärerregung im Schatten 
Nach (29) wird die Primärintensität durch die Abschattung nicht vollständig 
unterdrückt; wir müssen daher auch im Schattengebiet die Partikulärlösung einer 
inhomogenen Integralgleichung berücksichtigen, deren Lösung in erster Näherung 
die Inhomogentität selbst ist, nämlich: 
A(®) = <®< 2m. (39) 
Insgesamt haben wir im Schattengebiet: 


u(P)=C* 


2Aot ikasin® 
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1 — A(2 ef (40) - 
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§ 4. Diskussion des Ubergangsgebietes 


Als Ubergangsgebiet bezeichnen wir die Umgebung der Schattenränder. Dort 
müssen die für Licht- und Schattengebiet abgeleiteten Lösungen aneinander ange- 
schlossen werden. Die für das Schattengebiet abgeleitete Lösung ist allerdings 
wegen der Singularität von A(®) in einem kleinen Winkelgebiet (1/(ka)? an der 
Schattengrenze unbrauchbar. Die Kleinheit dieses Winkelbereiches berechtigt 
uns jedoch, die homogene Lösung des Schattengebietes stetig an die Lösung des 
Lichtgebietes anzuschließen. Dadurch bestimmt sich die Konstante C des An- 
satzes (9) nach (25) zu: 


C = 1,390 A, (! +): (41) 


Der Fehler, den wir dabei machen, muß, da ja die wahre Lösung eine stetige Funk- 
tion ist, mit wachsendem ka gegen null gehen. Eine genauere Aussage über den 
Fehler läßt sich jedoch nur schwierig gewinnen, da die Lösung im Übergangsgebiet 
sehr unbequem zu handhaben ist. 


$5. Das skalare Beugungsproblem mit der Randbedingung u= 0 


Bei der Lösung dieser Aufgabe gehen wir ganz analog dem vorigen Falle vor. 
Aus der Primärerregung einer in der y-Richtung einfallenden Welle ergibt sich 
zunächst: 


= B, sin® e-ikasin®, (42) 
P 
Als Lösung der Integralgleichung (2) im Lichtgebiet des Zylinders setzen wir an 
B@) e-ikasin® O<® <a. (43) 


Im Schattengebiet setzen wir an: 
öu(P) _ p ei "+ B2a—®) 


On, B(2 2) (44) 
Die Integralgleichung (2) wird für die Zylinderoberfläche 
éu(P) -ikasin® 
2sin® e (45) 
k ® ! 


Im Lichtgebiet ergibt sich für Punkte, die nicht in der Nähe der Schattengrenze 
liegen, bei analoger ng wie in § 2 


BM) = 2 By sin® + (46) 
In der Schattengrenze wird bis auf Glieder ~ (ka)"*: 
~ ~ 4T () e' me 84 913 in/é 
BO) = Bin) = — = 4 
(0) = Be) = (47) 
Im Schattengebiet ergibt sich analog zu $3 
Bei le ( 
B(®) = e (=) (48) 


(38) 
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mit folgender Bestimmungsgleichung für ß: 


3 oo 
ePt-P dt. (49) 


Dies ist konjugiert komplex zu der Bestimmungsgleichung (34) für x; die Lé- 
sungen vom kleinsten Betrage ergeben sich somit zu (s. (35)): 


= of; Bz = a9. (50) 
Von diesen beiden Lésungen ist nur die zweite brauchbar, da nur sie auf eine ge- 
dämpfte Welle führt. Damit wird: 
B®) 5,35 295 € (#) (51) 
Die durch den Rest der Primarerregung im Schattengebiet bestimmte inhomogene 
Integralgleichung wird näherungsweise gelöst durch: 


2 Boi 
BO) = — (52) 


Die gesamte Lisung im Schattengebiet wird damit 


2u(P) 5,35 295 e~ (<*) ! ka)ı® 2) 5,35 295 e ( ira) (on ®) 
1— B(22) (53) 
2 Boi e- ikasin® 
kasin? ® 


Durch den Anschluß im Ubergangsgebiet bestimmt sich die Konstante D zu 


0,84 913 e 7/6 


D = B(0) = 


o 


$6. Beugung elektromagnetischer Wellen am ideal leitenden Zylinder 


Für eine Primärwelle, deren magnetisches Feld Hz, parallel zur Zylinderachse 
schwingt, identifizieren wir die im vorstehenden verwendete Wellenfunktion u 
mit H, und bestimmen die Ausstrahlung nach einem Aufpunkte P(R,®,Z) 
außerhalb des Zylinders aus!): 

ia (R, ®, a, 9) 
4 ca 


6 


u(Q) dy. (54) 


Fiir den anderen Polarisationsfall (H, | Zylinderachse) setzen wir > =H und 


haben für die Ausstrahlung, sofern der Aufpunkt weit vom Zylinder entfernt ist: 


2a 
b ia öu(Q) 2 
éR ény dp. (55) 
0 
Wir werten die Integrale (54) und (55) nach der Methode der stationären Phase 
aus, wobei wir R> a voraussetzen. Dann können wir zunächst die asymptotische 
Entwicklung der Hankelfunktion einführen und 


YR—a— 2a R cos (vy—®) ~ R—a cos (p—®) 
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setzen. Wir erhalten: 
=> 2a 
Hi, — | cikR S cos (¢—®) e-ika cos (p-®) u(Q) dy (56a) 


anf Q 
= e- a J e- ika cos (9-®) (56 b) 


Die Auswertung der Integrale im Limes großer ka liefert Beiträge von vier 
Stellen stationärer Phase, nämlich von je einer Glanzstelle im Licht- und Schatten- 
gebiet, sowie von den beiden Stellen, von welchen die Kriechwellen tangentiell 
zum Beobachter abstrahlen. 


a) Glanzstelle im Lichtgebiet, 9, = + + = 


+00 a ® 2 
Hs, = sin( + 5) Ali f 5 


-» | 


=—2 sin (2 +$) (—2ka sin(7 + R 


16 ka sin (7+) | 
| | 
b) Glanzstelle im Schattengebiet, 
kasin’ (7 + 2) | 2 Rein (+5 
i kasin (* + 2) | 2 Rsin (7+) 


€) Beitrag der vom oberen Schattenrand herkommenden Kriechwelle, 7, =P —— 


ak 1,390 A, 


» — /_— 
Hz, | 1 — A(2 2) 


f a,e = )'r- ika(sint-r) 
Hz, OER : 116 er 1 ++. | (59a) 
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- a, 7/6 ka\t +in/i2 


b e 
Ho, = — 3, Baa) 


Man kann abschätzen, daß der relative Fehler der angegebenen Formeln 
proportional 1/ka ist, sofern die Beobachtung nicht aus der Umgebung des Winkels 
G®= + erfolgt. Die Abstrahlung in dieser Richtung kann nur bei genauerer 


Behandlung des Ubergangsgebiets berechnet werden. Wir beabsichtigen, hierauf 
an anderer Stelle zuriickzukommen. 


$7. Vergleich mit dem Experiment 


Wir vergleichen die Theorie mit unveröffentlichten Messungen von A. Lim- 
bach?) aus dem Jahre 1944 für Rückstrahlung von 9,1-em-Wellen an Metall- 
zylindern verschiedener Radien. Für den Fall H || Zylinderachse ergibt sich aus 
(57a), (59a), (60a) folgende Formel: 


H’(Z,R)Y2kR | 4,-22(ka)d a) 5,142 ka) 
Ay 1— A(2 2) e'*"?* 
amplitude 
0 1 2 3 4 5 6 7 


Abb. 2. Amplitude der gebeugten Welle fiir magnetischen Vektor parallel Zylinderachse. 
Fein ausgezogen: geometrische Optik, stark ausgezogen: theoretische Kurve, ©: Meß- 
punkte von A. Limbach 


2) Wir danken Herrn Baurat Dipl.-Ing. A. Limbach für die freundliche Erlaubnis, 
die in unserer Hand befindlichen Photokopien seiner Ergebnisse zu benützen. Wie uns 
Herr Limbach mitteilt, sind die übrigen Unterlagen seiner damaligen Meßreihen infolge 
der Kriegsereignisse leider verloren gegangen. 


d) Beitrag der vom unteren Schattenrand herkommenden Kriechwelle, ,=® + > 
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Der erste Summand ist der Anteil von der Glanzstelle im Lichtgebiet, der zweite 
rührt von den Kriechwellen her. Der Anteil von der Glanzstelle im Schattengebiet 
kann im vorliegenden Falle vernachlässigt werden. Abb. 2 zeigt das theoretische Er- 
gebnis als Funktion von ka zusammen mit den Meßpunkten von Herrn Limbach. 


Der Gang im großen wird durch die mit Vka anwachsende geometrisch reflek- 


Gangoifferenz (2+n)a 
le 


geometrisch reflektierte Welle 
Kriechwelle 


Abb. 3. Laufweg und Gangdifferenz zwischen geometrisch reflektierter Welle und 
Kriechwelle 


tierte Welle bestimmt. Die dieser Parabel überlagerten Maxima und Minima 
rühren von der Interferenz mit den Kriechwellen her; ihre Periode entspricht ange- 
nähert der geometrischen Gandifferenz (2 + 2) a zwischen den beiderlei Wellen 
(s. Abb. 3). Die Meßpunkte stimmen bis herunter zu ka » 1 befriedigend mit der 
theoretischen Kurve überein. 


Die vorstehende Theorie geht im wesentlichen auf das Jahr 1945 zurück, 
und wurde durch die Meßergebnisse von Herrn Limbach angeregt. Die dabei 
zutage getretenen eigenartigen Maxima und Minima wurden damals von dem einen 
von uns (F.) als Interferenzen zwischen geometrisch reflektierter Welle und Kriech- 
welle gedeutet, und die Theorie der Kriechwellen mittels der etwa gleichzeitig 
mündlich mitgeteilten Integralgleichungsmethode von A. W. Maue (veröffent- 
licht 1949 1. c.) entworfen. Diese Theorie wurde jetzt von dem anderen von uns 
in der vorliegenden Form durchgeführt. Die Anwendung auf die Beugung an der 
Kugel wird in Kürze in einer weiteren Arbeit dargestellt werden. 


Münster i. W., Institut für theoretische Physik der Universität. 


(Bei der Redaktion eingegangen am 17. März 1952.) 
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Uber An- und Abklingen der Photoleitfahigkeit 
von (dS-Einkristallen 


Von J. Faßbender und B. Seraphin 
(Mit 13 Abbildungen) 


Inhaltsiibersicht 


An sechs CdS-Einkristallen wurde das An- und Abklingen der Leitfähigkeit 
bei Belichtung mit Rechteckimpulsen verschiedener Intensität und Wellenlänge 
oszillographisch gemessen. Das Anklingen folgt überwiegend einer hyperbolischen 
Tangensfunktion. Die Abklingkurven setzen sich aus einem hyperbolischen 
Anfangsteil und einem exponentiellen Restverlauf zusammen. Mit wachsender 
Belichtungsintensität nimmt die Steilheit des Anklingens sowie der hyperbolische 
Anfangsteil des Abklingens zu, die Neigung des exponentiellen Restverlaufs er- 
wies sich als von der Belichtungsintensität unabhängig. Abhängigkeiten von der 
Wellenlänge und der Saugspannung wurden nicht gefunden. 

Die Ergebnisse lassen sich durch Annahme eines aus einem mono- und einem 
bimolekularen Vorgang zusammengesetzten Rekombinationsmechanismus quan- 
titativ deuten. Für den monomolekularen Rekombinationskoeffizienten werden 
Werte zwischen 30 und 200 see”, für den bimolekularen Werte zwischen 0,5 - 10712 
und 1,2 - 10-12 cm? sec”! gemessen. 


$ 1. Einleitung 


Zur Analyse des Leitungsmechanismus photohalbleitender Substanzen werden 
allgemein zwei Verfahren angewandt: Halleffektsmessungen und Leitfähigkeits- 
messungen bei nichtstationärer Belichtung. Halleffektsmessungen liefern bei 
richtiger Meßmethodik weitgehend hypothesenfrei direkte Aussagen über die 
Elektronenkonzentrationen und Beweglichkeiten, von denen die Leitfähigkeit 
direkt bestimmt wird. Dagegen ist es vor allem bei komplizierterer Reaktions- 
kinetik schwierig, aus den experimentellen Daten weitergehende Rückschlüsse 
auf die Elektronenbilanz zu ziehen, da nur Gleichgewichtszustände gemessen 
werden. Hier führen prinzipiell Messungen des nichtstationären Verhaltens, wie 
An- und Abklingmessungen der Leitfähigkeit, Glow-Kurvenmessungen und ähn- 
liche weiter. Allerdings bedarf es zu ihrer Interpretation schon ziemlich detaillierter 
Modellvorstellungen. Die erhaltenen Resultate hängen stark davon ab, welches 
spezielle Modell man ihrer Deutung zugrunde legt. 

Zur Beschreibung der Leitungsvorgänge in Halbleitern hat sich das Bänder- 
modell in letzter Zeit allgemein durchgesetzt; auch für die Photoleitung ist es 
verschiedentlich mit Erfolg verwandt worden, da es eine einfache Darstellung 
der in der klassischen Terminologie als Sekundärströme bezeichneten Vorgänge 
erlaubt. Die Grundzüge des Modells lassen sich etwa folgendermaßen skizzieren: 
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In dem Halbleiter sind mehrere Arten von Ladungsträgern (Elektronen, 
Löcher, Ionen) vorhanden, die einen Stromtransport bewirken können. Von den 
Ionen abgesehen, werden sie entweder durch thermische oder durch photoelektri- 
sche Ablösung erzeugt und durch Rekombination oder Anlagerung für den Strom- 
transport ausgeschaltet. Ablösung und Anlagerung kann aus dem Grundband 
oder aus Gitterfehlstellen erfolgen. Im stationären Fall stellt sich ein Gleich- 
gewicht zwischen Erzeugung und Vernichtung ein. Werden die Ladungsträger 
nun durch ein äußeres Feld in Bewegung gesetzt, so resultiert ein Strom. Die Ab- 
wanderung der Ladungsträger kann Konzentrationsverschiebungen verursachen, 
die zu Feldverzerrungen und Änderungen der Reaktionskinetik führen. Eine be- 
sondere Rolle spielen noch die Elektroden mit dem Sprung der Fermigrenze und 
den dadurch entstehenden Doppelschichten. 

Die mathematische Formulierung dieses Modells macht keine Schwierigkeiten, 
ist aber in voller Allgemeinheit bis heute noch nicht versucht worder. Man gelangt 
zu einem partiellen Differentialgleichungssystem, das als Unbekannte die Träger- 
konzentrationen und die Feldstärke, als Parameter die Dichte und energetische 
Lage der Störterme, die zugehörigen Übergangswahrscheinlichkeiten, Beweg- 
lichkeiten der Ladungsträger und die Bestrahlungsstärke enthält und streng nicht 
zu lösen ist. 

Die Vereinfachung des Systems zum Zwecke der näherungsweisen Lösung ist 
im wesentlichen in zwei Richtungen möglich. Die Art der Prozesse, deren Be- 
trachtung im Vordergrund steht, entscheidet über die Art der Vereinfachung: 

Stehen Leitungsvorgänge im Vordergrund, so beschränkt man sich meist auf 
den stationären Fall. Dann verschwinden die zeitlichen Ableitungen und das 
Differentialgleichungssystem reduziert sich auf ein gewöhnliches, das nur noch 
von der Ortskoordinate abhängig ist!) ?)®). 

Betrachtet man dagegen vorwiegend die Reaktionskinetik, so sieht man die 
Trägerkonzentrationen als ortsunabhängig an‘-*). Dann erhält man die 
Feldstärke durch Division der Spannung mit dem Elektrodenabstand, der Strom 
wird proportional der Trägerkonzentration und es bleibt nur noch ein Differen- 
tialgleichungssystem für die Trägerkonzentrationen zu lösen. Systeme dieser 
Art haben grundsätzlich die Form 


Die Indizes j und k laufen über alle betrachteten Termgruppen mit den Term- 
dichten N, bzw. N,. n, bzw. n, als Unbekannte des Systems bedeuten die Konzen- 
trationen in der j-ten bzw. k-ten Termgruppe, zwischen denen die Übergangs- 
wahrscheinlichkeiten y,, und ö,, bestehen. a, und 6, bedeuten Gewinn und Ver- 
lust an Trägern in der k-ten Termgruppe durch Einstrahlung. Auch dieses System 
ist nur näherungsweise lösbar und enthält dann noch eine Reihe experimentell 


1) K. Lehovec, Z. Naturforschg. 1, 258 (1946). 

2) J. Faßbender, Ann. Physik (6) 5, 33 (1949). 

3) F. Stöckmann, Z. Physik 128, 185 (1950). 

4) D. Blochinzev, Physik. Z. Sowjet-Union 12, 586 (1937). 

®) W. de Groot, Physica 6, 275 (1939). 

6) F. Möglich, R. Rompe, Physik. Z. 41, 236 (1940). 

‘) M. Schön, Verh. dtsch. physik. Ges. (3) 18, 70 (1939). 

») J. Faßbender, F. Möglich, R. Rompe, Ann. Physik (6) 3, 327 (1948). 
®») 1. Broser, R. Warminsky, Ann. Physik (6) 7, 289 (1950). 
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meist nicht bestimmbarer Parameter, die es nachträglich ermöglichen, die theore- 
tische Voraussage an den experimentellen Befund anzuschmiegen. Je nach dem, 
in welcher Richtung man das Differentialgleichungssystem (1) vereinfacht, d.h. 
welche Termgruppen man vernachlässigen zu können glaubt, lassen sich so aus den 
gleichen Experimenten sehr verschiedene Schlüsse ziehen. Die Modelle werden 
dabei häufig aus der Lumineszenzphysik übernommen, wo sich oft unmittelbare 
Beweise für die Existenz der einen oder anderen Termgruppe erbringen lassen. 
In welchem Umfang sich die an stark aktivierten Luminophoren gewonnenen 
Ergebnisse auf die viel weniger verunreinigten Einkristalle, die Gegenstand von 
Leitfähigkeitsuntersuchungen sind, übertragen lassen, ist in vielen Fällen fraglich. 
Gleichzeitige Untersuchungen von Leitfähigkeit und Lumineszenz sind nur ver- 
einzelt ausgeführt worden”). 

Die große Anzahl frei wählbarer Parameter, die Systemen der Art (1) den 
Charakter einer Interpolationsformel verleihen, lassen das Bedürfnis entstehen, 
weitgehend hypothesenfreie Aussagen zu gewinnen!!) oder die Zahl der Para- 
meter auf ein Minimum herabzudrücken. 

Man kann etwa die Frage stellen, wieviel Parameter mindestens notwendig 
sind, um die einfachsten experimentellen Erfahrungen über die nichtstationären 
Vorgänge zunächst rein formal darzustellen. Die vorliegende Arbeit verfolgt 
dieses Ziel speziell im Falle des CdS und leitet eine Berechtigung für dieses Ver- 
fahren aus der Tatsache ab, daß CdS verhältnismäßig rein und gut einkristallin 
herstellbar ist. Es soll deshalb versucht werden, mit einer einzigen Differential- 
gleichung, die nur zwei Parameter enthält, auszukommen. Es zeigt sich, daß der 
Erfahrungskomplex, der mit diesem primitiven Ansatz bereits beschrieben werden 
kann, erstaunlich umfangreich ist. Die physikalische Interpretation dieser beiden 
Parameter ist dabei zunächst zurückgestellt. 


$2. Theoretischer Überblick 


Derartige Differentialgleichungen sind bisher in drei Typen verwandt worden: 


—an (2) 
(3) 
(4) 


(a Anzahl der pro cm und sec im Kristall absorbierten Quanten, x und ß Rekom- 
binationskoeffizienten). 

Dabei entspricht (2) einem monomolekularen, (3) einem bimolekularen und (4) 
einem gemischt mono- und bimolekularen Rekombinationsmechanismus. 


I. Der stationäre Fall 


Im Fall stationärer Belichtung (2 = 0) folgt bei Verwendung der in (2), (3) 
oder (4) vorgeschlagenen Mechanismen fiir die sich im Gleichgewicht zwischen Er- 
1) R. H. Bube, Physic. Rev. 88, 393 (1951). 
1) Lee Gildart and A. W. Ewald, Physic. Rev. 88, 359 (1951). 
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zeugung und Vernichtung einstellende Elektronenkonzentration 


My = — (2a) 
Vs (3a) 
2a 


Trägt man diese Beziehung zwischen absorbierter Belichtungsintensität und 
stationärer Elektronenkonzentration in doppeltlogarithmischen Koordinaten auf, 
so liefert der monomolekulare Mechanismus eine Gerade mit der Steigung 1, der 
bimolekulare eine Gerade mit der Steigung 0,5. 

Die auf (4a) fußende Beziehung des gemischten Mechanismus liefert in der 
doppeltlogarithmischen Darstellungsweise eine gekrümmte Kurve, die bei kleinen 
Intensitäten sich einer Geraden mit der Steigung 1, bei großen Intensitäten sich 
einer Geraden mit der Steigung 0,5 asymptotisch anschmiegt. Die Lage des Über- 
gangsgebietes in bezug auf die Intensität hängt von dem Verhältnis der Kon- 
stanten x und ß ab, die den Charakter des vorliegenden Mechanismus bezüglich 
seiner Mischung aus einem mono- und einem bimolekularen bestimmen. Definiert 
man als „Knick“ der Kurve (Abb. 3) den Punkt mit der Steigung 0,75, so rückt 
dieser Punkt mit wachsendem Anteil des monomolekularen Mechanismus zu 
größeren Intensitäten a’ hin, und zwar gemäß ‘ 


a TEE (5) 


Die Krümmung der Kurve an der Stelle des ,,Knicks“ bleibt von einer Ver- 
schiebung längs der Intensitätsachse unberührt. Untersucht man den Abstand 
der Punkte, in denen sich die Kurve bis auf 10% der Steigung ihrer Asymptoten- 
geraden genähert hat, so zeigt sich, daß diese Punkte ohne Rücksicht auf die durch 
(5) beschriebene Lage des Knicks stets um zweieinhalb Größenordnungen der In- 
tensität auseinanderliegen (Abb. 2). 


Il. Der nichtstationäre Fall 
Bei nichtstationärer Belichtung mit Rechtecklichtimpulsen der absorbierten 
Intensität a erhält man für das Anklingen der Leitfähigkeit nach Einsetzen der 
Belichtung bzw. für das Abklingen nach Abschalten der Belichtung beim mono- 
molekularen Mechanismus 


= in (Anklingen) (2b) 

= = eat (Abklingen) (2c) 
beim bimolekularen 

= — Tg (Ya ß: t) (Anklingen) (3b) 

= (Abklingen). (3e) 
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Der gemischte Mechanismus liefert kompliziertere Resultate, die deshalb etwas 
eingehender diskutiert werden sollen. Als allgemeine Lösung von (4) findet man 


t = 
“0 2p(1—Ce*') 
mit der Abkiirzung 
A=Yx:+4aß (7) 


und C als Integrationskonstante. 
Fiir die Anklingkurve (n (0) = 0) wird dies 


2a 


n(t) =~ 8 

oder unter Beachtung von (4a) 
a+ A (9) 


Mo a+ A-Ctg At 


Für die praktische Auswertung ist es nützlich, A durch die experimentell 
leicht meßbare Halbwertszeit t auszudrücken; es ist dann 


n 1 
(10) 
No 
2 A 
t= ea: 
Mit 
1 
Ar Tq x 
folgt dann schlieBlich 
_ 
oder 
& | 
1] 
A e4*—8 ( ) 


Bei bekanntem x und r läßt sich hieraus A berechnen. 
Die Gleichungen lassen sich erheblich vereinfachen, wenn 


(12) 

ist. Dann geht wegen 

hod 47 < 1 

A Vat+4aB 
(9) über in 

2-2 t; 

No 2 
und wegen (12) in 

= Iq(Vap-t). (3b) 


Unter der Voraussetzung (12) klingt also auch der gemischte Mechanismus im 
wesentlichen bimolekular an. A läßt sich dann einfach aus 
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bestimmen. Bei verschiedenen absorbierten Intensitäten a, und a, verhält sich 
dann 


Umgekehrt kann man die Gültigkeit von (3b) als Kriterium für das Zutreffen 
von (12) ansehen. 

Für das Abklingen der Leitfähigkeit findet man mit den Randbedingungen 
n(0) = %, a=0 


n 1 
14 
m +1 er 
Für kleine Zeiten (! <=) geht dies iiber in die bimolekulare Abklingkurve 
n 1 
my t+’ 
und fiir groBe Zeiten (t > +) in die monomolekulare 
a. (2c) 
Die Lage des Knickpunktes t, = < hängt nur von x, nicht aber von ß oder der 


Bestrahlungsintensität ab. 
Für die Anfangssteigung des Abklingens erhält man 


= — nf. (15) 


§ 3. Apparative Anordnung 


Der Kristall wurde mit Impulsen von etwa 50 Hz belichtet, die in iiblicher 
Weise durch einen rotierenden Sektor erzeugt wurden, bei gleicher Dauer von 
Belichtung und Dunkelpause. Die auf diese Weise unausweichliche Abrundung 
der Impulse wurde mit einer Alkali-Photozelle gepriift. Es ergab sich ein Ab- 
rundungsgebiet von 10% sec Dauer, einem Zeitintervall also, das etwa fiinfzigmal 
kleiner als die fiir den Anklingvorgang und sogar tausendmal kleiner als die fiir 
den Abklingvorgang maßgeblichen Zeiten ist. Die erzeugten Impulse können also 
für die hier vorgetragenen Messungen als scharf rechteckig bezeichnet werden. 

Als Lichtquelle wurde eine Quecksilberhöchstdrucklampe HBO 200 verwendet, 
in den Strahlengang konnten bei Bedarf Monochromatfilter für die Linien 5461 AE 
(Ausläuferabsorption) und 4358 AE (Grundgitterabsorption) eingeschaltet werden. 
Der CdS-Kristall wurde über einen regelbaren Arbeitswiderstand (1...10&Q) 
mit einer stufenweise veränderlichen Spannung (7...35 V) versehen; die an dem 
Arbeitswiderstand auftretenden stromproportionalen Spannungen wurden auf die 
Eingangsbuchsen des Oszillographen gegeben. 

Um die Kurvenform des An- und Abklingens genauer analysieren zu können, 
braucht man vor allem für den Anklingvorgang eine sehr hohe Zeitauflösung. 
Bei den üblichen Oszillographen nimmt das eigentliche An- und Abklingen nur 
einen kleinen Teil der Schirmbreite ein, wenn man eine volle Periode schreibt. 
Synchronisiert man mit einem höheren Vielfachen der Lichtimpulsfrequenz, so 


man 
i a 
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erhält man viele sich überlagernde Kurvenzüge, deren Auswertung nicht sehr 
genau ist. Deshalb wurde folgender Weg gewählt: 


Der verwendete Oszillograph'?) besaß ein impulsgesteuertes Kippgerät. Da- 
durch kommt man von der sonst bestehenden Kopplung zwischen Schreibge- 
schwindigkeit und Periodendauer frei und kann irgendeinen beliebigen Bruchteil 
der Periode über den ganzen Schirm schreiben. Bei den Anklingmessungen wurde 
etwa eine hundertstel Periode auf den ganzen Schirm geschrieben. Die Größe 
der Zeitauflösung wird charakterisiert durch den zeitlichen Abstand zweier Punkte 
des Schirmbildes, der, sofern nichts anderes angegeben wird, bei den Ankling- 
messungen 5,5 - 10-* sec, bei den Abklingmessungen 0,65 - 10”? sec betrug. Die 


s 8 Kompensationsspannung 


zum Wehnelt- Zylinder 


zZ 


Abb. 1. Schema der Versuchsanordnung. L Lichtquellen, S Drehscheibe, Z, gemessener 
Kristall, Z, Steuerzelle, V Verstärker, K Kippgerät, Z Zeitmarke, B Braunsches Rohr 


für die Betätigung des Kippgerätes notwendigen Steuerimpulse wurden photo- 
elektrisch mit einer zweiten CdS-Zelle erzeugt, die über den gleichen Sektor eben- 
falls mit Rechteckimpulsen belichtet wurde. Durch eine mechanische Schwenk- 
vorrichtung konnte der Phasenwinkel zwischen den beiden Rechteckimpulsen 
beliebig verändert und so die ganze Periode über den Bildschirm geschoben werden. 


Um Nichtlinearitäten des Oszillographen in Richtung der Zeitachse zu eli- 
minieren, steuerte ein Zeitmarkengeber den Strahl hell-dunkel. Die Frequenz 
dieser Hell-Dunkel-Steuerung konnte zusammen mit der Schreibgeschwindigkeit 
des Strahles stufenweise verändert werden und lieferte auf diese Weise die oben 
beschriebene veränderliche Zeitauflösung der Vorgänge. Zeitmarke und Kipp- 
gerät wurden in jeder Periode durch den Steuerimpuls synchron angestoßen. 


Die Ablenkspannung in Ordinatenrichtung wurde durch eine Kompensations- 
spannung gemessen, indem jeder Kurvenpunkt durch eine meßbare Gegenspannung 
an den Ablenkplatten auf eine Nullinie zurückgebracht wurde. Auf diese Weise 
fallen die durch das Rohr verursachten Nichtlinearitäten in Ordinatenrichtung 
völlig heraus. 


Der stark gegengekoppelte Gegentaktverstärker besaß eine völlig geradlinige 
Frequenzkurve von 16 Hz bis 3 MHz und steuerte den Schirmdurchmesser völlig 
verzerrungsfrei aus. Das Meßverfahren liefert ohne Umweg über das Schirm- 
bildphoto in 5 Minuten etwa 40 zeitlich äquidistante Meßpunkte. 


. 12) Für die Konstruktion des Oszillographen sind wir Herrn Palitsch zu Dank ver- 
pflichtet. 
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sehr § 4. Diskussion der MeBergebnisse 
An sechs CdS-Zellen der üblichen Bauart!3) wurden An- und Abklingmessungen 
5 Da- sowie Messungen des Photostroms in Abhängigkeit von der Bestrahlungsintensität 
fen durchgeführt. Die auf den Kristall auftreffende Strahlungsdichte betrug in der 
htei 
urde 
TOBe 
nkte 
ling- 
Die 
ner 
ohr 
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ık- 
en 
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p- 
S- Abb. 2. Schirmbild einer An- und einer Abklingkurve. Aus Gründen der Darstellung 
ig wurde die Zeitauflösung geringer gewählt alsim Text angegeben sowie der Durchmesser 
se der Punkte vergrößert 
mit „7 cm“ bezeichneten Entfernung von der Lichtquelle 10: Quanten/cm? sec 
e für die Linie 4358 AE. Variationen der Bestrahlungsintensität wurden durch 
| Veränderung des Abstandes zwischen Zelle und Blende erreicht. Der Dunkel- 
‘ strom lag bei allen sechs Zeilen mindestens 5 Größenordnungen unter dem 
| Photostrom. 
Die in den Reaktionsgleichungen (2)—(4) auftretende Größe a bezieht sich 
nicht auf die auftreffende Intensität, sondern auf die im Kristall wirksam zur An- 
18) J. Faßbender, l.c. 
Ann. Physik, 6. Folge, Bd. 10 26 
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regung verbrauchte Quantenzahl pro cm? und sec. Von der auftreffenden Inten- 
sität ist also zunächst der reflektierte Anteil abzuziehen und der Rest a, auf eine 
Eindringtiefe zu verteilen, die durch den Absorptionskoeffizienten 4. angegeben 
wird. Korrekt wäre hier der Ansatz a, e~“* für die in der Tiefe x vorliegende Quan- 
tendichte. Da jedoch der Ansatz einer einzigen Reaktionsgleichung der Typen 
(2) bis (4) ohnehin nur eine nullte Näherung darstellen kann, wurde es als aus- 
reichend erachtet, den eindringenden Anteil a, homogen über eine Schicht von 
10 cm Dicke zu verteilen, entsprechend dem Absorptionskoeffizienten u = 
1,5 - 105cm~ nach Messungen von Mollwo) für das Gebiet der Grundgitter- 
absorption. Der Reflektionskoeffizient wurde nach Messungen von Mühlschlegel®s) 
mit 0,3 angenommen. Sofern im folgenden also von der Bestrahlungsintensität a 
die Rede sein wird, handelt es sich stets um die Dichte der wirksam zur Anregung 
verbrauchten Quanten, die in einer Schicht von 10cm Dicke als gleichbleibend 
angenommen wird. 


I. Der stationäre Fall 


Die Messungen des Photostroms unter dem Einfluß stationärer Belichtung, 
die über drei Größenordnungen variiert wurde, lieferten die bekannten Ergeb- 
nisse!®) 17) 18), In doppeltlogarithmischer Auftragung des Photostroms über die 
Intensität ergab sich für große Intensitäten genähert eine Gerade mit der Stei- 
gung 0,5, die in einem „Knick“, der sich über eine knappe Größenordnung der 
Intensität erstreckt, in eine genäherte Gerade der Steigung 1 übergeht. 

Gilt für den Leitungsvorgang die Gleichung 


j=enb-€ (16) 


und vernachlässigt man ortsabhängige Effekte, so lassen sich die gemessenen 
Kurven der Photoleitfähigkeit mit den nach (4a) berechneten Kurven für die 
stationäre Elcktronenkonzentration n, vergleichen, sobald man bestimmte An- 
nahmen über die Beweglichkeit b macht. Bei den sechs gemessenen Kristallen 
läßt sich mit Beweglichkeitswerten zwischen 50 und 200 cm?/Volt sec eine an- 
nähernde Übereinstimmung erzielen. Die Werte sind im Mittel um einen Faktor 2 
größer als die von Faßbender und Lehmann”) und Bube®°) gefundenen. Dies 
ist nicht verwunderlich, da bei Halbleitern die Beweglichkeit stark vom Störstellen- 
gehalt abhängt. Wie weit durch eine Konzentrationsabhängigkeit der Beweglich- 
keit zusätzliche Abweichungen vom theoretisch zu erwartenden Verlauf der Leit- 
fähigkeitskurve vorgetäuscht werden, ist nicht ohne weiteres zu entscheiden, da 
eingehende Messungen hierüber noch nicht vorliegen. Nach unveröffentlichten 
Hall-Messungen von Diedrich scheinen solche Abhängigkeiten vorzuliegen; 
allerdings sind die Messungen ohne Potentialelektroden ausgeführt worden und 
deshalb mit den unvermeidlichen Unsicherheiten durch Potentialsprünge an den 
Elektroden usw. behaftet. 

Die unter der Annahme einer Beweglichkeit der obigen Größe und Eiger- 
schaft aus (4c) berechneten Photostromkurven ähneln in ihrer Form stark den ge- 


14) Mollwo, Reichsber. f. Ph 1,1 

15) B. Mühlschlegel, Ann. 

16) J. FaBbender, I. c. 

17) J. Broser u. R. Warminsk 

16) E. A. Niekisch, Ann. Physik os, 291 (1951). 

19) J. Faßbender u. H. Lehmann, Ann. Physik 6) 6, 215 (1949). 
20) R. H. Bube, |. c. 


hysik (6) 9, 29 (1951). 
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messenen. Der Knick liegt befriedigend in der durch die Beziehung (5) angegebenen 
Größenordnung, ist bei der gemessenen Kurve jedoch wesentlich schärfer als aus 
der n,-Kurve zu erwarten war (Ubergangsbereich von errechneten zweieinhalb 
Größenordnungen gegenüber dem gemessenen Übergangsbereich von einer Größen- 
ordnung)?!). Als asymptotische Werte für die Steigung der gemessenen Photo- 
stromkurven erhält man jedoch nicht 1 und 0,5, sondern meist etwas kleinere 
Werte. Immerhin sind gelegentlich auch Kurven gefunden worden, die zumindest 
asymptotisch eine recht gute Übereinstimmung mit der Theorie aufweisen*). 
Die Abweichungen sind zweifellos eine Folge der eingangs erwähnten Vernach- 
lässigungen, die bei der Aufstellung der Reaktionsgleichung (4) gemacht werden. 
Jedenfalls läßt sich sagen, daß sich die Meßkurve verhältnismäßig gut an die des 
gemischten Mechanismus anschließt und in der Tat sind ja auch erstmalig aus 
solchen Messungen Schlüsse auf die Existenz eines zusätzlichen monomolekularen 
Reaktionsmechanismus gezogen worden. 


ot 
B — 

10 10" 0” 
Abb. 3. Die stationäre Elektronenkonzentration n, als Funktion der Bestrahlungsinten- 
sität a. ———— nach (2a) errechnet (monomolekular); —-—-+—- nach (3a) 


errechnet (bimolekular); - - nach (4a) errechnet (gemischt); °°°°° mit einer 
Beweglichkeit 6 = 100 em?/Volt sec m Photostrommessung berechnet (Kristall 
r. 3) 


II. Nichtstationäre Vorgänge 
Das An- und Abklingen der Leitfähigkeit wurde untersucht bezüglich 
1. Der Kurvenform, d. h. des der Rekombination zugrunde liegenden Mechanis- 
mus, 


21) E. A. Niekisch, 1. c. 
2) J. Broser u. R. Warminsky,l.c. 
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2. Der Abhängigkeit von der Intensität der Bestrahlung. 
3. Der Abhängigkeit von der Wellenlänge der Bestrahlung. 
4. Der Abhängigkeit von der an die Zelle gelegten Saugspannung. 


II,1. Die Kurvenform 


Das Anklingen läßt sich in seinem überwiegenden Teil streng durch die hyper- 
bolische Tangensfunktion 


%q(c-t) (17) 
No 
beschreiben. 


b 
02} 


5 10 15 20:55:10 sec 
Abb. 4. Das Anklingen der relativen Elektronenkonzentration (Kristall Nr. 6) 


Die Abklingkurve setzt sich aus zwei Teilen zusammen: Einem steilen, hyper- 
bolischen Beginn, der während einiger 10-* sec wirksam ist, und einem Verlauf 
streng nach einer Exponentialfunktion mit einer Abklingkonstante a zwischen 
10? und 200 sec je nach Kristall. 

Der Zeitpunkt des Überganges zwischen beiden Bereichen wird gut durch die 
Beziehung 


(18) 


wiedergegeben, so daB das hyperbolische Anfangsstiick um so ausgedehnter ist, 
je flacher das nachfolgende exponentielle Abklingen vor sich geht. Bei sehr steilem 
exponentiellen Abfall ist es nur noch in hoher Zeitauflésung zu erkennen. 

In den Bereichen hoher Elektronenkonzentration, d.h. im Ausläufer der An- 
kling- und im Anfangsteil der Abklingkurve zeigen sich bei der Mehrzahl der unter- 
suchten Kristalle Abweichungen von dem soeben beschriebenen Verhalten: 

a) Bei etwa 65—75proz. Auffüllung der stationären Elektronenkonzentration 
weicht die Anklingkurve im Sinne eines Sättigungseffektes von der hyperbolischen 
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Tangensfunktion ab, die relative Konzentration wächst langsamer an. Der Einsatz- 
punkt dieses Sättigungseffektes ist von Kristall zu Kristall verschieden und ver- 
schwindet mitunter ganz, als Folge der Vorbehandlung. Er ist von der eingestrahl- 
ten Intensität und somit von der Zeit, die seit dem Einschalten der Belichtung 
verflossen ist, völlig unabhängig. 
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Abb. 5. Abklingen der relativen Elektronenkonzentration in einfach-logarithmischer 
Auftragung. Zum Nachweis der Hyperbelnatur ist der Anfangsteil in höherer Zeitauf- 
lösung rechts oben doppeltlogarithmisch aufgetragen. (Kristall Nr. 5) 


b) Die hyperbolischen Anfangsteile der Abklingkurven besitzen in der Mehrzahl 
der untersuchten Fälle Hyperbelexponenten, die etwas kleiner sind als 1 mit 


der unteren Grenze 0,8. 
Das überwiegend bimolekulare Anklingen 


und das überwiegend monomolekulare Ab- 


klingen läßt sich nur mit einem gemischten 
Mechanismus deuten, bei dem außerdem 
die Bedingung (12) erfüllt ist. Bis auf einen 
Faktor 3 ist diese identisch mit der Lage 
des Knickpunktes der stationären Kurve, 
d.h. die Intensität muß so hoch gewählt 
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08} 
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werden, daß man sich im oberen Zweig der a 
Kurve der Abb. 3 befindet. 3 
Die beiden Teile der Abklingkurve wer- 
den durch GI. (14) gut wiedergegeben, auch 
der Zeitpunkt des Überganges zwischen 
dem hyperbolischen Anfangsteil und dem 
exponentiellen Restverlauf liegt in der rich- € J . 
tigen Größenordnung und wird nur von höhere a zeigt. (Kristall 
der Abklingkanstanten x, nicht dagegen von nr 
der Anregungsintensität oder von der Rekombinationskonstanten ß bestimmt. 
Erstaunlich ist hierbei die Gültigkeit der Differentialgleichung (4), die doch 
bestimmt nur eine Annäherung darstellt, für einen so großen Bereich von x. Bei 


+ 
5 10: 0.65-10°sec 
Abb. 6. Abklingkurve eines Kristalls 
mit sehr steilem exponentiellen Abfall. 
Hyperbolischer Abfall nur noch 
während etwa 10? sec Dauer, wie 
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einem extrem trägen Kristall mit einer Abklingzeit von Minuten ließ sich die ganze 
Meßkurve durch Gl. (14) wiedergeben unter Verwendung eines x in der Größen- 


08} 


067 


02} ---- 


i i t 
5 10 15 20-5510 sec 
Abb. 7. Anklingkurve mit Sättigung 


ordnung 10? sec-! und der Lage des Ubergangsbereiches bei etwa 2 Minuten. Ge- 
messen wurde die Kurve mit Stoppuhr und Milliamperemeter®). Damit er- 
streckt sich der Giiltigkeits- 
N % bereich über 4 Größenordnun- 
gen von «. 


ol 11,2. Die Abhängigkeit der Vor- 
gänge von der Bestrahlungs- 
EN intensität 
y=-095 Die Steilheit der Ankling- 
: kurven wächst gleichsinnig mit 
der Intensität der Bestrahlung, 
hohe Bestrahlungsintensität lie- 
Abb. 8. Die hyperbolischen Anfangsteile der Ab- Schnellem Auffüllen der statio- 
klingkurven mit einem von Eins abweichenden nären Elektronenkonzentration, 
Hyperbelexponenten. (Kristalle 2 und 4) geringe Intensität flache Kur- 
ven mit längerer Anklingzeit. 
Formal läßt sich dies Verhalten gut mit einer Veränderung des Faktors c in der 
hyperbolischen Tangensfunktion (17) beschreiben, und zwar verhalten sich die 
Quadrate der c wie die Intensitäten a 


y = -080 


(20) 


oflich, Für die Überlassung dieser Meßkurve sind wir Herrn Muscheid zu Dank ver- 
pflichtet. 
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Da die Bedingung (12) bei den Messungen ausreichend erfüllt war, läßt sich (20) 
aus der Gl. (13) verstehen. 


Der hyperbolische Teil der Abklingkurve ist um so steiler, je größer die ein- 
gestrahlte Intensität ist. Da der exponentielle Teil des Abklingens von der 


nt Photostrom (willkiirl Einheiten ) 
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Abb. 10. Anklingkurven als Funktion der Bestrahlungsintensität (Kristall Nr. 3) 


Intensität unabhängig ist, ergibt sich somit eine Parallelverschiebung der Ab- 
klingkurven als Funktion der Bestrahlungsintensität. 

Dies Verhalten des hyperbolischen Anfangsteils wird durch (14) gut wieder- 
gegeben, da mit steigender Anregungsintensität auch n, wächst und die dadurch 
bewirkte Erhöhung der Hyperbelsteilheit eine Parallelverschiebung des nach- 
folgenden exponentiellen Teil nach sich zieht. 
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11,3. Die Abhängigkeit von der Wellenlänge der Einstrahlung 


Die An- und Abklingvorgänge zeigen keine Abhängigkeit von der Wellenlänge 
der Einstrahlung. Monochromatische Bestrahlung im Bereich der Ausläufer- 


Abb. 11. 
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Abklingkurven als Funktion der Bestrahlungsintensitat (Kristall Nr. 2 und 5) 
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Abb. 12. Anklingkurve in Abhängigkeit von der Saugspannung. (Kristall Nr. 1) 


absor 
der w 
ist. 
kling 
sehr 
> 
span 
[ 
gene 
venf 
tats 
wied 
Abh 
Well 
stra! 
Zelle 
spar 
men 
sung 
hier 
gun 
dun 
t weis 
1 
— 
1 sis? 
kon 
Tei 
suc 
~~ 
— 
| 
Au 


nge 
fer- 


J. Faßbender u. B. Seraphin: Photoleitfähigkeit von CdS-Hinkristallen 389 


absorption (5461 AE) liefert gegeniiber der Einstrahlung im Gebiet der Grund- 
gitterabsorption nur diejenige Veränderung der Kurven, die für die Erniedrigung 
der wirksamen Quantendichte durch Verkleinerung der Absorption charakteristisch 
ist. Sowohl die Anklingkurven als auch der hyperbolische Anfangsteil der Ab- 
klingkurven verlaufen also gemäß der sehr herabgesetzten wirksamen Intensität 
sehr viel flacher, bleiben aber in ihrer charakteristischen Form erhalten. 


11,4. Die Abhängigkeit von der Saugspannung 


Sowohl An- als auch Abklingkurven sind von der_an die Zelle gelegten Saug- 
spannung unabhängig. 


Derin (4) vorgeschla- me 

gene Ansatz, der Kur- ast ty 
venform und Intensi- 4 é 
tätsabhängigkeiten gut 4 
Abhängigkeiten vonder 
Wellenlänge der Ein- 
strahlung und der an die Saugspannung 
Zelle gelegten Saug- 02 + 75 Volt 
spannung. Im Zusam- x 16 Volt 
menhang mit den Mes- + 35 Volt 

1 2 t 
sungen ergibt sich auch 5 0 15 065-10 


hieraus eine Berechti- ii ae in Abhängigkeit der $ 

um di 2 . 18. ngkurve in ngigkeit von der Saug- 
55 = en spannung. (Kristall Nr. 2). Die Meßkurven der Abb. 12 
ung des näherungs- nd 13 geben zugleich eine Darstellung der bei der Messung 
weisen Ansatzes (4). etwa erzielten Reproduzierbarkeit 


$5. Bestimmung der Konstanten 


I. Bestimmung von « 
x läßt sich direkt aus der Abklingkurve entnehmen. Für große Zeiten ("> =) 


geht das Abklingen in einen streng monomolekularen Prozeß mit der Abkling- 
konstanten « über. « ergibt sich somit einfach aus der Neigung des exponentiellen 
Teiles der Abklingkurve in einfachlogarithmischer Auftragung. Fiir die unter- 
suchten Kristalle ergaben sich die Werte 


Kristall Nr. «+ 102 sec-t KristallNr. & + 1072 sec7! 
1 0,71 5 0,35 
2 0,42 6 0,54 
3 1,70 Träger Kristall | 10-4 
4 0,78 (Muscheid) 


II. Bestimmung von ß 


Mit bekanntem « läßt sich die Gültigkeit der Voraussetzung (12) nachprüfen. 
Aus der Halbwertszeit + der Anklingkurve läßt sich mittels (11) der numerische 
Wert der Konstanten A berechnen. Man findet Werte in der Größenordnung 108 
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Da «? höchstens von der Größenordnung 10% ist, muß (7) fast ausschließlich vom 
Glied 4 a ö herrühren. Bei bekanntem a läßt sich dann ß aus der Beziehung 


$= (Vaß-r) 
bestimmen. Um die Genauigkeit zu erhöhen, wurde für jeden Kristall die An- 


klingkonstante A für mehrere Intensitäten bestimmt. Über eine Streuung von 
10% ergaben sich hierbei für 8 die Mittelwerte. 


Kristall Nr. B+ 1012 cm? sec-1 KristallNr. B - 1012 cm? sec 
1 0,6 4 | 0,5 
2 1,2 5 1,1 
3 0,8 6 0,8 


Daß Faßbender und Lehmann*) bei ihrer 8-Bestimmung aus Wechsel- 
lichtmessungen auf Werte in der gleichen Größenordnung kamen, trotz Zugrunde- 
legung eines rein bimolekularen Rekombinationsmechanismus, ist plausibel, denn 
die dort verwendeten Belichtungsfrequenzen lagen in der Größenordnung von 
einigen Kilohertz, ließen den Kristall also niemals bis in den monomolekularen 
Teil hinein abklingen. 


III. Bestimmung von n, 


Zur Abschätzung der stationären Elektronenkonzentration n,, die sich als 
Folge einer bestimmten Bestrahlungsstärke a im Leitfähigkeitszustand befindet, 
bieten sich drei Möglichkeiten an, von denen jedoch nur die dritte von den ersten 
beiden unabhängig ist, da die erste und zweite über die Differentialgleichung (4) 
miteinander verknüpft sind und die Übereinstimmung ihrer Ergebnisse nur als 
Kontrolle zu werten ist. 

III,1: Die Beziehung (15) liefert eine Möglichkeit, das Produkt n,ß aus dem 
Anfangsteil der Abklingkurve abzulesen. Bei bekanntem ß ergibt sich z. B. für 
den Kristall Nr. 1 bei Bestrahlung in 7 cm Abstand 


N, = 0,7 - 1015 cm”®, 


III,2: Aus der Gl. (4a) ergibt sich andererseits unter Verwendung der be- 
kannten Größen & und ß für die Intensität a, die der Entfernung von 7 cm ent- 
spricht, die stationäre Elektronenkonzentration n, zu 


My = 2+ 1015 


III,3: Eine hiervon unabhängige Möglichkeit der Abschätzung von n, bieten 
die Photostrommessungen als Folge stationärer Belichtung. Die Gültigkeit der 
Leitungsgleichung (16) erfordert, daß man den gemessenen Photostrom an die 
berechneten n,-Werte anschließen kann unter Verwendung einer Beweglichkeit, 
die in der Größenordnung der anderweitig gemessenen Werte liegt. Einsetzen 
einer Beweglichkeit in der Größenordnung 10? cm? Volt— sec-!, die den Messungen 
von Faßbender und Lehmann*), Bube**) und Diedrich?) entspricht, 


2) J. Faßbender u. H. Lehmann,l.c. 
25) J. Faßbender u. H. Lehmann, |. c. 
2) R. H. Bube, l.c. 

27) Unveröffentlichte Diplomarbeit. 
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führt wieder zu dem Wert 
N, = 1015 em”? 


für die stationäre Elektronenkonzentration bei Photostrommessungen im Ab- 
stande von 7 cm. 


$6. Diskussion 


Erstaunlich ist die gute Übereinstimmung zwischen Theorie und Experiment, 
wenn man sich vor Augen hält, daß es sich bei der Differentialgleichung (4) be- 
stimmt nur um eine nullte Näherung handeln kann. Vernachlässigt sind die Orts- 
abhängigkeit der Trägerkonzentration — insbesondere Elektrodeneffekte — und 
die Konzentrationsabhängigkeit der Beweglichkeit. Die bei Anregung im Gebiet 
der Grundgitterabsorption sehr schnell exponentiell abfallende Anregungsdichte 
wurde durch eine von gleicher Anregungsdichte erfüllte Schicht von 10° cm 
Dicke ersetzt. Der Einfluß der Störterme wurde summarisch durch ein lineares 
und ein quadratisches Glied berücksichtigt (direkte Band-Band-Übergänge der 
Rekombination sollten quantenmechanisch sehr unwahrscheinlich sein). Des 
weiteren wurde angenommen, daß das Gleichgewicht der Randschicht sich in 
einer Zeit einstellt, die kurz ist gegenüber den Auf- und Abbauzeiten des Elektronen- 
plasmas. Andernfalls würde sich den gemessenen An- und Abklingkurven die Auf- 
und Abbaueffekte der Randschicht überlagern. Die Berechtigung für diese An- 
nahme wurde aus der guten Übereinstimmung mit dem Ansatz (4) entnommen. 


Angesichts dieser Vernachlässigung ist der von diesem primitiven Ansatz 
wiedergegebene experimentelle Tatbestand erstaunlich umfangreich. Umgekehrt 
läßt sich daraus entnehmen, daß man aus Experimenten der vorstehenden Art 
nur mit großer Vorsicht weitergehende Schlüsse auf Einzelheiten der Reaktions- 
kinetik ziehen darf. 

Über das Anklingen der Leitfähigkeit bei CdS lagen exakte Messungen bisher 
nicht vor. Messungen von Lee Gildart und A. E. Ewald) benutzen die Stei- 
gung des ersten Teiles der An- und Abklingkurven, um die Beweglichkeit der 
Elektronen des Leitungsbandes zu bestimmen. 


R. W. Smith**) führt Anklingmessungen an CdS mit dem Ziele aus, eine 
„Lebensdauer“ der Elektronen im Leitfähigkeitsband aus der Halbwertszeit des 
Anklingens zu bestimmen. Mit einer Quantendichte von 5 - 1013 erzielt er Halb- 
wertszeiten in der Größe 10"? sec. Schätzt man aus diesen Angaben ß etwa nach 
(3b) ab, so gelangt man zu guter Übereinstimmung mit der bei unseren Kristallen 
gefundenen Größenordnung. Die von Smith auf eine stationäre Belichtung auf- 
gesetzten Rechteckimpulse verursachten eine um so steilere Anklingkurve, d.h. 
eine um so geringere Halbwertszeit, je stärker die stationäre Belichtung war. Diese 
Beobachtung ist unmittelbar aus dem in $2 entwickelten zu verstehen, sobald 
man etwa in (3b) a durch n, ersetzt, etwa vermittels (4a). In gleicher Weise läßt 
sich die Verschiebung der dort abgebildeten Abklingkurve mit zunehmender 
stationärer Belichtung deuten, da von Smith das gleiche Verhalten des Abklingens 
wie in unserer Abb. 11 beobachtet wird, allerdings ohne auf Kurvenform oder 
ähnliches einzugehen. Für die Beweglichkeit findet Smith aus Hallmessungen 
Werte zwischen 40 und 200 cm? Volt”! sec—!. Die von Smith benutzten Kristalle 


28) Lee Gildart and A. E. Ewald, l.c. 
2) R. W. Smith, RCA Review Vol. XII, Nr. 8, 350 (1951). 
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müssen in stärkerem Maße verunreinigt sein als die unseren, da Smith Potential- 
abfälle von 90% der gesamten Saugspannung in unmittelbarer Elektrodennähe 
findet. Die in unserer Arbeit verwendeten Kristalle stammen aus einer Fertigung, 
für die bereits früher durch Institutsmessungen an reichlichem Material sicher- 
gestellt war, daß der Potentialabfall an den Elektroden 20% der Gesamtspannung 
nicht überstieg und das Potential dazwischen streng linear verlief. 


Das Abknicken der Anklingkurve bei höheren Konzentrationen ist schon von 
De Groot beim Anklingen der Phosphoreszenz beobachtet worden **). Er suchte es 
als Überlagerung verschiedener Anklingkurven mit abnehmender Anregungs- 
intensität zu deuten. Infolge des hohen Absorptionskoeffizienten im Gebiet der 
Grundgitterabsorption nimmt die absorbierte Quantenzahl mit der Eindringtiefe 
sehr schnell ab. Wenn sich die unterschiedlichen Elektronenkonzentrationen nicht 
durch Diffusion oder Leitung ausgleichen, überlagern sich beim Anklingen die 
verschieden stark angeregten Bereiche und täuschen ein Art Sättigungseffekt vor. 


Für das Anklingen der Leitfähigkeit trifft diese Deutung mit Sicherheit nicht 
zu, da einmal der Effekt auch bei Anregung im Gebiet schwacher Absorption auch 
vorhanden ist, wo die absorbierte Quantenzahl in der ganzen Kristalldicke prak- 
tisch die gleiche ist, zum anderen der Effekt sich beim gleichen Kristall durch die 
Vorbehandlung ändern läßt und bei manchen gar nicht auftritt. Aus dem gleichen 
Grunde kann es sich auch nicht um einen durch die Apparatur vorgetäuschten 
Effekt handeln. Beim Abknicken der Kurve vor dem Erreichen des Sättigungs- 
wertes bleibt zu beachten, daß die Trägerkonzentration asymptotisch durchaus 
den gleichen Wert erreicht, wie bei strenger Befolgung der hyperbolischen Tangens- 
funktion. Dies geht daraus hervor, daß zur Angleichung der gerechneten an die 
gemessene Kurve die Halbwertszeit r benutzt wird. Die Halbwertszeit 7 bestimmt 
sich aber aus dem asymptotischen Wert der Funktion für große Zeiten. Der Sät- 
tigungseffekt kann also nicht durch zeitunabhängige Effekte gedeutet werden, 
wie z. B. durch eine Konzentrationsabhängigkeit der Beweglichkeit, sondern nur 
durch vorübergehend mitbeteiligte Prozesse, die nicht in den Rahmen der Näherung 
fallen. 


Bei der Abklingkurve treten genau wie beim Anklingen bei hohen Konzen- 
trationen Abweichungen vom theoretischen Verlauf auf. Der hyperbolische An- 
fangsteil entspricht nicht exakt einem Hyperbelexponenten —1, sondern liegt 
zwischen —0,8 und —1,0. Ähnliche Beobachtungen sind schon früher beim Leucht- 
abklingen der Kristallphosphore gemacht worden. 


Bisher sind über die Natur der Prozesse, denen der monomolekulare, bzw. 
der bimolekulare Prozeß entspricht, keine besonderen Annahmen gemacht worden. 
Um zu einer Aussage über die physikalische Natur dieser Prozesse zu gelangen, ist 
eine Zuordnung zu den Gliedern der Differentialgleichung (4) erforderlich. Auf 
Grund eines vor allem von der Lumineszenzphysik gelieferten umfangreichen Tat- 
sachenmaterials wird allgemein der bimolekulare Prozeß einem leuchtenden, der 
monomolekulare Prozeß einem strahlungslosen Übergang zugeordnet, so daß sich 
die Differentialgleichung für das Abklingen in der symbolischen Form schreiben 


läßt 
(Breve (ai) leuchtend + (21) 


#) W. de Groot, Physica 7, 432 (1940). 
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Unter dieser Annahme würde ß der Ubergangswahrscheinlichkeit des leuch- 
tenden, bimolekularen Prozesses entsprechen. Fiir diese Annahme spricht die 
größenordnungsmäßige Übereinstimmung mit der in der Lumineszenzphysik 
bekannten Konstanten ß. 


Der monomolekulare Prozeß ist seiner Natur nach noch nicht näher bekannt. 
Wegen seiner Spannungsunabhängigkeit läßt er sich nicht als Rekombination 
der Träger mit den Elektroden (Absaugen) deuten. Die gesamte, beim Abklingen 
fließende Elektrizitätsmenge ist zu groß, um die Stromleitung während des ex- 
ponentiellen Teils durch eine Nachlieferung aus Haftstellen zu bestreiten. Wahr- 
scheinlich handelt es sich um eine, durch Verunreinigung induzierte monomolekulare 
Übergangsmöglichkeit ®!) ®). Diese Annahme wird gestützt durch die Möglichkeit, 
beim gleichen Kristall x durch die Vorbehandlung reversibel zu verändern. 


Es bleibe jedoch dahingestellt, ob der durch ein lineares Glied beschreibbare 
Rekombinationsanteil wirklich durch einen selbständigen Prozeß verursacht wird, 
der proportional der Elektronenkonzentration verläuft, oder ob die Linearität 
nur eine formale Konsequenz der Vereinfachung des Differentialgleichungssystems 
(1) darstellt. Mit Sicherheit kann jedoch der pseudolineare Effekt ausgeschlossen 
werden, der trotz rein bimolekularer Recombination entstehen würde, falls ein 
sehr großer Bestand n, an Dunkelelektronen durch Photoelektronen größen- 
ordnungsmäßig gleicher Konzentration erhöht wird, gemäß der Gleichung 


=a —B(m + n) n. (22) 


Diese Deutung scheidet aus, weil der Dunkelstrom bei sämtlichen gemessenen 
Kristallen mindestens 5 Größenordnungen unter dem Photostrom lag. 


Als Folge dieser Zuordnung der Prozesse sollte die Lumineszenz für größere 
Zeiten monomolekular abklingen. Die Leitfähigkeitsmessung zeigt, daß für große t 
der Elektronenbestand ausschließlich durch den monomolekularen Prozeß be- 
stimmt wird. Er läßt sich in der Form 


n(t) = (23) 


schreiben, wobei n{ durch Extrapolation des exponentiellen Teils auf t = 0 zu 
bestimmen wire. 


Der leuchtend iibergehende Anteil, der in diesem Gebiet nach Vorhergehendem 
sehr viel kleiner ist als der strahlungslose, beträgt definitionsgemäß 


= — 


dn = #2 


also 


31) E. J. Adirowitsch, Ber. Akad. Wiss. USSR 68, 111 (1948); 68, 635 (1948); 
Nachr. Akad. Wiss. USSR, Physik. Ser. 18, 101 (1949). In deutscher Sprache erschienen 
in Abhandlungen aus der sowjetischen Physik, Folge 1, S. 71, Verlag Kultur und Fort- 
schritt 1951. 

3) E. A. Niekisch u. R. Rompe, Z. physik. Chem. 198, 200 (1951). 
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Die Leuchtintensität sollte also exponentiell abklingen. Solche Kurven sind ven 
Garlick und Wilkins®) an ZnS:Ag 10-* Ni sowie an ZnSCdS:Cu auch ge- 
messen worden. 

Es steht jedoch außer Frage, daß die zuletzt angeschnittenen Probleme der 
Zuordnung physikalischer Prozesse zu den formalen Rekombinationsgliedern 
sowie des Zusammenhanges zwischen Leitfähigkeit und Phosphoreszenz einer 
weitaus gründlicheren Bearbeitung an Hand eines umfangreicheren experimen- 
tellen Materials bedürfen. 


"Herrn Prof. Dr. F. Möglich danken wir für viele fördernde Diskussionen, 


®) Garlick u. Wilkins, Proc. Roy. Soc. London A 184, 408 (1945). 


Berlin-Buch, Institut für Festkörperforschung der Deutschen Akademie 
der Wissenschaften zu Berlin. 


(Bei der Redaktion eingegangen am 9. April 1952.) 
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dern 
a Uber elektrische und akustische Einschwingvorgange 
nen- 
Von H. Dänzer 
nen, (Mit 9 Abbildungen) 
Inhaltsiibersicht 
; Die Einschwingvorgänge auf rückgekoppelten elektrischen und akustischen 
nn Leitungen werden mathematisch untersucht. Die nicht harmonischen Eigentöne 


der mit induktiver Belastung versehenen Leitung laufen bis zu einer durch die 
Dämpfung bestimmten Ordnungszahlgrenze an, verschwinden aber bei nicht 
linearer Charakteristik wieder. Stationär bleibt nur die Eigenfrequenz bestehen, 
für die die Rückkopplungsbedingung am besten erfüllt ist. 


Den Anlaß zu der vorliegenden Untersuchung bildeten die schönen Experimente 
von F. Trendelenburg und Mitarbeiter!) über die Einschwingvorgänge bei 
Orgelpfeifen. Die Verfasser haben bei ihren Untersuchungen an Orgeln aus der 
Barockzeit auf die merkwürdigen Klangphänomene hingewiesen, die dem einge- 
schwungenen Zustand vorausgehen und die für den Klangcharakter dieser Re- 
gister und die Klangschönheit der betreffenden Pfeifen so wesentlich sind. In 
folgendem soll die Grundlage für eine mathematische Theorie dieser Vorgänge 
| gegeben werden. Aus mehr prinzipiellen Gründen wollen wir uns dabei auf das 

elektrische Ersatzschaltbild Abb. 1 stützen, das die wesentlichen Züge des akusti- 
| schen Phänomens mit umfaßt und den Vorteil größerer Übersichtlichkeit besitzt. 

In Analogie zu den Verhältnissen bei den Pfeifen ist in Abb. 1 der Schwingkreis der 
Rückkopplungsschaltung durch eine am anderen Ende geschlossene Leitung ersetzt. 

Um die Analogie der Schwingungsvorgänge in einer Pfeife mit den elektrischen 
Vorgängen in Abb. 1 im einzelnen darzutun, betrachten wir zunächst die Zungen- 


Abb. 1. Prinzipschaltbild. Die Komponente i des Anodenstromes J geht in die Leitung, 
die Komponente J läuft über die Induktivitat L 


1) F,Trendelenburg, F. Thienhausu. E. Franz, Akust. Z. 1, 59 (1936); 8, 7 (1938). 
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pfeife Abb. 2. Die Zunge Z läßt in das Innere der Pfeife einen Luftstrom J ein- 
treten, dessen Stärke von der jeweiligen Zungenelongation abhängt und im übrigen 
Sättigungscharakter hat. Die funktionelle Abhängigkeit zwischen dem Luft- 
strom J und dem Schalldruck 2, im Innern der Pfeife entspricht also durchaus der 
Charakteristik einer Elektronenröhre. Da der Einlaß der Luft in die Zungenpfeife 
von dem Druck p, auf die Zunge abhängt, entspricht dieser der Gitterspannung. 
Durchläuft der Druck ein Maximum, so strömt bei hochabgestimmter Zunge eine 
maximale Luftmenge in die Pfeife ein, ganz analog wie in der elektrischen Schaltung 
Abb. 1, bei der ein negatives Spannungsmaximum an der Klemme K, über die 
Gittersteuerung einen starken Elektronenzufluß zu X, bewirkt, wenn die Rück- 
kopplung richtig gepolt ist. Die potentielle Energie der akustischen oder elek- 
trischen Schwingung wird so im Moment ihres 
Maximums künstlich noch verstärkt, so daß 
Verluste, die sich im Verlauf der Schwingung 
einstellen, ausgeglichen werden können, was 
zur Erhaltung einer stabilen, ungedämpften 
Schwingungsform notwendig ist. 

Ähnlich liegen die Verhältnisse bei den 
Lippenpfeifen der Orgel, die in Abb. 3 dargestellt 


i sind. In dem Lippenspalt Z der Pfeife schwingt 

. eine bestimmte Luftmasse im Takte der Schwin- 

Pi gung hin und her. Diese Luftmasse repräsen- 
tiert die „Induktivität“ des Schwingkreises. 
"ZI ze Pi An die „Induktivität“ Z schließt sich im 

7 L 77, Pfeifeninneren der eigentliche Schwingkreis an, 
der in der elektrischen Schaltung durch das 

Lechersche Drahtsystem dargestellt wird und 

in der Akustik dem Pfeifenrohr entspricht. 

Abb. 2. Abb. 3 Da am offenen Ende der Pfeife der Druck in 
Zungenpfeife Lippenpfeife . erster Näherung gleich Null gesetzt werden 


kann und die elektrische Leitungsspannung dem 
Druck entspricht, so ist in Abb. 1 das Lecher'sche Drahtsystem, um die Analogie 
vollständig zu machen, am anderen Ende geschlossen. Die aus dem Kernspalt 
der Lippenpfeife austretende Luft vermischt sich mehr oder weniger mit der im 
Lippenspalt hin- und  herschwingenden Luftmasse. Für die Unterhaltung der 
Schwingung ist es auch hier notwendig, daß die Luftzufuhr ins Innere der Pfeife 
in dem Moment am größten ist, wo der Druck p, im unteren Ende des Pfeifen- 
rohres sein Maximum durchläuft. Diese Forderung ist dadurch verwirklicht, 
daß die Luftmasse im Lippenspalt eine freie Schwingung ausführt, bei der die 
Elongation und die wirkende Kraft bekanntlich eine Phasenverschiebung von 
180° gegeneinander besitzen. (Dieser Phasenverschiebung wird in der bekannten 
Dreipunktschaltung durch die Überkreuzung der Leitungen, die nach der Kathode 
bzw. dem Gitter führen Rechnung getragen.) Die Stromverzweigung der drei 
Ströme 


ImJ+i 


ist im elektrischen und akustischen Fall (Abb. 1 u. 3) dargestellt. J ist der ,,Anoden- 
strom‘, J fließt über die Induktivität Z zur ,,Erde“ bzw. in den AuBenraum ab, 
i ist der Leitungsstrom, der ins Lechersystem bzw. in das Pfeifenrohr übertritt. 
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Wir werden zunächst annehmen, daß der Ruhepunkt des Anodenstromes in 
der Mitte der Charakteristik liegt, so daß für mäßige negative oder positive Gitter- 
spannung die Röhrencharakteristik als gradlinig angesehen werden kann. Die 
Wirkung der Krümmung bzw. des Knickes der Charakteristik soll weiter unten 
untersucht werden. Wir fassen das Schwinggebilde als „komplexen Wechsel- 
stromwiderstand“ auf, an dessen Klemme K, der Anodenstrom J eintritt und bei 
K, austritt. 

Wir betrachten eine zeitliche Verteilung des Anodenstromes /(t’), die sich in 
der Vergangenheit abgespielt haben möge. Der gegenwärtig betrachtete Zeit- 
punkt, in dem der Strom und die Spannung V an den Klemmen K des komplexen 
Wechselstromwiderstandes betrachtet werden, wird mit ¢ bezeichnet. Die Zeit, 
um welche ein Ereignis zurückliegt, ist also 


T=t-t. 
Die drei Zeiten sind in der Zeitachse Abb. 4 eingetragen. Das betrachtete 


Ereignis, welches 7’ Sekunden zurückliegt, sei der Durchgang der Elektrizitäts- 
menge I(t’) dt’ durch den komplexen 


Widerstand R(w). Dabei entsteht an t2-00 \ t’ T 
dessen Klemmen ein Spannungsver- er + + ir 
lauf, dessen Ausläufer noch zur Zeit ¢ t=t=0 at 

wirksam sind. Uberlagert man alle Abb. 4. Zeitachse 


Spannungsverläufe, die zu verschie- 
denen Zeitpunkten dt’ in der Vergangenheit ausgelést worden sind, im Zeit- 
punkt ¢, so erhält man die an den Klemmen K gerade vorhandene Spannung V (t): 


t 
J I(t’) f(t—t’) dt’. (1) 


Hier bedeutet f(t—t') den Spannungsverlauf, den die Einheit der Elektrizitäts- 
menge, welche im Zeitmoment t’ übergeht, ausgelöst hätte. 

Bedenkt man, daß vermöge der Rückkopplung der Anodenstrom /(t’) eine 
eindeutige Funktion der Schwingkreisspannung V (t’) ist, so stellt Gl. (1) eine 
Integralgleichung dar, zu deren weiterer Behandlung der „Kern“ f(t — t’) bekannt 
sein muß, den wir als Fourierintegral darstellen wollen. Zu diesem Zweck be- 
trachten wir zunächst einen einzelnen Stromimpuls J mit der Dauer dt’ =r 
für sich isoliert. Die einzelne Stromzacke läßt sich im Sinne des Fourierschen 
Integraltheorems in ein Frequenzkontinuum auflösen: 


IT) = A(w) dan. (2) 


Für die Amplitudenverteilung A(w) ergibt sich hieraus 


A(w) = —— J I{T)ei®T AT. 3 
J (3) 
Der betrachtete Stromimpuls ist nämlich nur innerhalb des Zeitabschnittes r 
von Null verschieden. 
Durchsetzt die Stromzacke den komplexen Widerstand R(@), so schließt sich 
daran ein Spannungsverlauf dV an. 
+00 
av = md A(w) R(w) ei®T dv. (4) 


Ann. Physik. 6. Folge, Bd. 10 27 
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Jede Frequenzkomponente der Spannung wird nach (4) durch Multiplikation der 
Stromkomponente A (w) mit R(w) erhalten. Die Zulassung negativer Frequenzen 
hat nur formale Bedeutung. Es gilt 


R(— w) Ro), 


wobei die Uberstreichung den Übergang zum komplex-konjugierten Wert bedeutet. 

Das Integral (4) konvergiert nur, wenn R(w) für unendlich hohe Frequenzen 
hinreichend schnell gegen Null konvergiert. Diese mathematische Voraussetzung 
ist ohne Beschränkung der Allgemeinheit erfüllt, weil R(w) immer einen zwischen 
seinen Klemmen liegenden Kondensator C enthält, der die Anodenkapazität der 
Röhre und die Kapazität der Drahtzuleitungen verkörpert (vgl. Abb. 6). Für kleine 
Frequenzen kann von der Existenz dieses Kondensators abgesehen werden. 

Bei der Auswertung des Integrales (3), dessen Wert wir in (4) einzusetzen ge- 
denken, brauchen wir also nur Frequenzen unterhalb einer bestimmten Schranke 
Q in Betracht zu ziehen. Für solche Frequenzen kann in (3) die Exponential- 
Funktion durch 1 ersetzt werden, wenn wir so klein wählen, daß 


Qr<l. 
Dann wird also 
A(o) = (5) 


und wir erhalten durch Einsetzen in (4) und Integration über alle Zeitelemente r 
auf der Zeitachse (Abb. 4) zwischen t’ = — oo und !' = 


+00 


Vi) = 5 I(t’) J R (ww) dey dt’. (6) 
Der Vergleich mit (1) liefert für den gesuchten ‚Kern“ 
+00 
fe—r)= f dw. (7) 
-© 


Eine weitere experimentell stets erfüllte Voraussetzung für die Möglichkeit der 
Fourierschen Integral-Darstellung (7) ist das zeitliche Abklingen dieser Funk- 
tion. Durch die stets vorhandene Dämpfung ist diese mathematische Notwendigkeit 
immer erfüllt. 

Für das Folgende ist es zweckmäßig, die Kombination i ® durch p zu ersetzen 
und R als Funktion von p zu betrachten. Für unsere mathematische Diskussion 
ist es entscheidend, für welche komplexen Werte für p die 
Funktion R(p) Pole besitzt. Polstellen mögen existieren 
fx für p,, Pr, p,... In der Reihe dieser komplexen Fre- 
(x quenzen treten auch die konjugiert komplexen Größen auf, 
so daß die Lage der Pole, wie in Abb. 5 dargestellt, gegeben 
ist. Indem wir in (7) die Integration über p aus der imagi- 


fx naren Achse herausziehen und um die Pole schlieBen, er- 
én hält man unter Einführung der „Residuen“ ß, der Pole 
Bn = dp (8) 
nach (7): 
p-Ebene +00 
Abb. 5 (9) 
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also nach (1)?) 
VQ= ZB, ae. (10) 


Setzt man also nach (10) abkürzend für die Spannung V (t) 


Vit) = = Bn Galt) ernt, (11) 


so ergibt sich fiir die ,,.komplexen Amplituden“ 


t 
= f I(t’) erat dr. (12) 


An Stelle von (12) erhalten wir durch Differentiation nach ¢ folgenden Differential- 
gleichungen fiir die 9,: 


Wir identifizierten in unserem Fall den Strom /(t) mit dem Anodenstrom der 
Rückkopplungsschaltung (Abb. 1). Dann bedeutet V (t) die Spannung an den Klem- 
men der Induktivität Z. Die Summendarstellung (10) für die Schwingkreisspan- 
nung bedeutet offenbar eine Superposition der auf die einzelnen Eigenfrequenzen 
des Schwingkreises entfallenden Spannungsverläufe, die von den zwischen t’ = — 00 
und t/ = ¢ übergegangenen Elektrizitätsmengen angestoßen werden. In der Tat 
haben wir ja die komplexen Frequenzen 


P,=io, (14) 
definiert durch die Pole des komplexen Wechselstromwiderstandes, für die dieser 
Widerstand unendlich groß wird. Diese Eigenschaft besitzt bekanntlich ein in 
der Sperrkreisschaltung befindlicher Schwingungskreis. Die komplexen Fre- 
quenzen p, bedeuten also die Eigenschwingungszustände des Schwing- 
kreises. 

Die Tatsache der Rückkopplung (vgl. Abb. 1) liefert mathematisch einen Zu- 
sammenhang zwischen der Stromstärke J (t) (Anodenstrom) und der Spannung V (t), 
von der ein bestimmter Bruchteil © über die induktive Kopplung am Gitter liegt 
und auf Grund der Röhrencharakteristik auf diese Weise den Strom / bestimmt. 
Dabei wird es freilich nicht gleichgültig sein, ob die Amplituden der Spannung 
V(t) sich auf den linearen Teil der Charakteristik beschränken oder von vorn- 
herein ein nichtlinearer Zusammenhang zwischen Spannung und Strom berück- 
sichtigt werden muß. Wir beschränken uns zunächst auf den linearen Fall, der 
durch den Ansatz gekennzeichnet werden kann 


I(t) =OS V(t). (15) 


Hierin bedeutet S die Röhrensteilheit und © ein echter Bruch, der angibt, welcher 
Teil der Spannung V(t) als Gitterspannung über den Rückkopplungsmechanismus 
wirksam wird. Ehe wir mit diesem Ansatz in (13) eingehen, ist es zweckmäßig, 


2) Ein einfaches Modell für (10) sind n hintereinander geschaltete Schwingkreise 
den Kapazitäten C, und den Eigenfrequenzen iw, = p, für die nach (8) folgt 
in, = 
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an Stelle der komplexen Amplituden, die Größen ®, durch die Definition einzu- 
führen 


D,=gp,emt, (16) 
mit denen (11) übergeht in 
+00 
Vi) By. (17 
Die Differentialgleichung (13) geht mit (16) über in 
d® 
PnP, = Il). (18) 
Unter Benutzung des linearen Ansatzes (15) 
I(t) =O8 (19) 
folgt dann 
d® 
PnP, + OS SB, ®,. (20) 
Ausgeschrieben erhalten wir das lineare System 
d®, 
= (OSB. + +OSB,D,+--- 
OSB,®, + (OSB, + p,)®, + (21) 


Die Lésung des linearen Differentialgleichungssystems (21) liefert bekanntlich 
der Ansatz 

= C, *, (21a) 
den wir in das Schema (21) eintragen. Wir erhalten dann lineare homogene Glei- 
chungen zur Bestimmung der Amplituden C,. Die Forderung nicht trivialer 
Lösungen dieses Systems führt bekanntlich zum Verschwinden seiner Determinante, 
womit eine Bestimmungsgleichung für die P gewonnen ist: 


OSB, OS py--- 
0= OSB, OSf,---|. (22) 


Berücksichtigt man, daß in unserem Fall unendlich vieler Eigenfrequenzen die 
„Säkular-Determinante“ unendlich viele Lösungen zuläßt, so nimmt die allge- 
meine Lösung des Differentialgleichungssystems (21) an Stelle von (21a) die 
Form an: 


= Pe, (23) 


Der allgemeinen Struktur nach ist dieses Ergebnis wohl bekannt. Es bedeutet, 
daß jede Eigenschwingung ,,Kopplungsschwingungen“ ausführt. Die Koppel- 
frequenzen P, sind die Lösungen der Säkulardeterminanten. Trägt man das 
Ergebnis (23) und (17) ein, so erscheint die Schwingkreisspannung V (t) ebenfalls 
in der Form 
V(t) = (24) 
v 


jedoch mit anderen Amplituden A,. 
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Es handelt sich also darum, die komplexen Frequenzen P, und die Ampli- 
tuden A, von (24) zu berechnen. Die p, sind die Pole des komplexen Wechsel- 
stromwiderstandes R(p) und als solche im allgemeinen ebenfalls komplex (vgl. 
Abb. 5): 

Pr = Pon— En: (25) 


Pp, ist rein imaginär und es bedeutet 9,, = i@,, die Eigenfrequenz des ge- 
dämpften Schwingsystemes, welche bei schwacher Dämpfung mit der ungedämpf- 
ten Eigenfrequenz praktisch übereinstimmt. Die Dämpfung des Schwingkreises, 
die zu dem Realteil —e, von p, Veranlassung gibt, denken wir uns als Parallel- 
widerstand w, (vgl. Abb. 6) berücksichtigt, der alle dämpfenden Ursachen der 
h-ten Eigenschwingung repräsentieren möge. Die Dämpfungskonstante e, in (25) 
hängt also — wie wir gleich sehen werden — mit w, unmittelbar zusammen. 

Da nach (25) in p, noch die Dämpfung des Schwingkreises enthalten ist, 
setzen wir an Stelle von (11), unter Einführung der y,, die definiert sind durch 


(26) 
= rat, (27) 


(27) bringt zum Ausdruck, daß die y, im eigentlichen Sinn als „komplexe Ampli- 
tuden‘‘ anzusprechen sind, da der Exponent in (27) jetzt rein imaginär ist. Aus 
(13) folgt für die nach (26) definierten Größen y, die Differentialgleichung 


= I(t) e-Pat—e, y,. (28) 


Nehmen wir hinzu, daß auf Grund der Réhrencharakteristik der Strom /(t) eine 
Funktion @ der nach (26) gegebenen Spannung darstellt: 


T(t) = G(T Bryn er), (28a) 
so ist die Lésung unseres Problems in den Gln. (28) und (28a) enthalten. Im nicht- 
linearen Fall, wo @ eine beliebige Funktion sein kann, ist allerdings eine mathe- 
matische Diskussion der Lösung kaum möglich. Wir beschränken uns deshalb 
zunächst auf eine lineare Abhängigkeit zwischen J und V, für die wir eine Lösung 
von (28) und (28a) suchen. Bevor wir eine Näherungslösung für die y, angeben, 
berechnen wir zunächst die Dämpfungskonstanten ¢, der einzelnen Eigenschwin- 
gungen und ziehen zu diesem Zweck die Definition (8) der ß, heran. Wir setzen hier 


1 1 
Rp) = gay = — 
Rp) (p - 
so daß aus (8) folgt 


By (29) 


* 
dp 
Das ungedämpfte System habe den komplexen Widerstand R,(p), so daß nach 
Abb. 6 zu setzen ist 


+=. (30) 


en 


Also 
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Wir bedenken, daß bei schwacher Dämpfung e,<|p,n| gesetzt werden kann 
= 0. (32) 
Andererseits folgt auf (30), der Definition für p, entsprechend, 


O = + 4 


Durch Subtraktion der beiden letzten Gleichungen ergibt sich (bei Vernachlässigung 
von Gliedern der Größenordnung e&,) 


so daß (29) unter Berücksichtigung von (31) übergeht in 


Br = En (33) 
Trägt man nun den linearen Ansatz (15) in (28) ein, so folgt unter Berück- 
sichtigung von (26) und (33) 


d 1 
(© By Yq Pend! (34) 


In der Summe fehlt das Glied mit g = h und die Exponenten der Exponential- 
funktion sind rein imaginär. Da wir bei schwacher Rückkopplung die Amplituden 
y, als „langsam veränderlich‘‘ betrachten dürfen, ist in (34) die Durchführung 
einer über mehrere Perioden sich erstreckenden Mittelung erlaubt, bei der N 
in erster Näherung verschwindet. Dann folgt also angenähert (die strenge Lösung 
für den linearen Fall geben wir weiter unten) 


1 
Auf Grund von (26) und (35) liefert ein Vergleich mit (24) in derselben Näherung 
P, = + (36) 
Ir) 
vie) l 
Je 


x=0 


Abb. 6. Leitung mit Induktivität Z, Kapazität C und Dämpfungswiderstand w, 


wo die ß, durch ihre Definition (8) gegeben sind. Zur Ermittelung der ß, betrachten 
wir das Schaltbild Abb. 6. Verstehen wir unter Z den Wellenwiderstand der 
Leitung, so igs für die komplexe Leitfähigkeit Rp) 


R-1(p) = pC + — Pi) + = (37) 
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und unter Hinzunahme der Definition der p, als Pole von R(p): 


Durch Subtraktion der beiden letzen Gleichungen ergibt sich 
1 l —i 
(89) 


worin die p, auf Grund von (38) festliegen. Die Berechnung des Residuums (8) 
kann auf Grund von (39) vorgenommen werden und liefert 
+e 1 1 
Bei vernachlässigbar kleiner Dämpfung und Kapazität C folgen aus (40) und (38) 
c 


mit der Definition 2», = + 


L 

(41a) 
14375 (1 + 

(41b) 


Die Bedingung (41 b) liefert den Imaginärteil von (36), während (41a) den Realteil 
x, von P, festlegt, der also die Form annimmt 


08-.- 


ax, = Las i 42 

Für die Spannung V(t) folgt damit nach (24) die Lösung 
V(t) = =A Ment, (43) 


(42) und (43) kénnen sowohl im elektrischen als auch akustischen Fall diskutiert 
werden. Im Falle der Orgelpfeifen muß das akustische Analogon des für (42) 
charakteristischen Widerstandsverhältnisses 

oL 

zZ 
bekannt sein. Es hängt von den konstruktiven Daten der Pfeife ab, wie weiter 
unten gezeigt wird. 


Zur Deutung von (42) und (43) betrachten wir — als „negativen Widerstand“, 


der in Abb. 6 dem positiven Widerstand w, parallel geschaltet ist. Die Entdämp- 
fung der Eigenschwingung ®,, funktioniert nach (42) natürlich nur, solange 


1 
,>0 (44) 


positiv ist. Werden die Schwingungsamplituden immer größer, so muß (nach 
Barkhausen) mit einer mittleren Steilheit gerechnet werden, die infolge der 
Krümmung der Charakteristik stets kleiner ausfällt als die Steilheit im linearen 
Teil. Bezeichnen wir die mittlere Steilheit mit S, so unterbleibt eine weitere An- 
fachung der Schwingung, wenn 


(44a) 


n 


1 
(33) 
ick- 
den 
= 
1 
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geworden ist. Da die Dämpfung der Eigenschwingungen sehr verschieden sein kann 
und bei hohen Frequenzen infolge der Abstrahlung sowohl im elektrischen als 
auch im akustischen Fall mit wachsender Ordnungszahl A schnell ansteigt, mußten 
wir jeder Eigenschwingung einen individuellen Dämpfungswiderstand w, geben, 
der mit wachsendem A abfällt. Es gibt also eine obere Schranke für die Ordnungs- 
zahl h, oberhalb derer eine Selbsterregung nicht mehr stattfinden kann. 

Bedenkt man, daß die mittlere Steilheit im zeitlichen Ablauf der Anfachung 
immer mehr absinkt, so verschiebt sich diese Schranke für dauernd nach unten, 
so daß Eigentöne, die zunächst angelaufen waren, wieder verschwinden, bis 
schließlich die Eigenschwingung mit der geringsten Dämpfung, im allgemeinen 
die Grundschwingung allein übrig bleibt und den eingeschwungenen Zustand be- 
stimmt. Da ihre Schwingungen im allgemeinen nicht genau sinusförmig erfolgen, 
tritt dann ein System von harmonischen Obertönen zum Grundton hinzu, die von 
den Eigentönen wohl zu unterscheiden sind. 

Die nicht harmonischen Eigentöne, deren Frequenzen durch (41b) gegeben 
sind, treten also nur im Verlauf der Anfachung hervor, nicht aber im eingeschwun- 
genen Zustand. Der Verfasser ist der Meinung, daß dieser Mechanismus bei der 
Erklärung der Trendelenburgschen Experimente mit herangezogen werden 
muß. Komplizierend kommt dabei hinzu — was von Trendelenburg zur 
Verständlichmachung der „Vorläufer“ allein herangezogen wird —, daß beim 
Öffnen des Pfeifenventils der Luftdruck eine gewisse Zeit braucht, um seinen 
Endwert zu erreichen (was in elektrischer Analogie so viel bedeuten würde, 
daß die Kathode der Elektronenröhre langsam angeheizt wird). Wie lange dieser 
Vorgang dauert, hängt von dem benutzten Windladensystem ab (,,Schleiflade“ 
oder „Kegellade‘‘). 

Das Verschwinden zunächst angefachter Eigentöne, überhaupt die ganze 
zeitliche Entwicklung des Klangspektrums im nichtlinearen Fall, ist natürlich in 
unseren Gln. (28) und (28a) ebenfalls enthalten, aber — wie schon bemerkt — 
mathematisch sehr schwer zu diskutieren. Im Fall schwacher Rückkopplung er- 
gibt sich jedoch eine prinzipielle Vereinfachung dadurch, daß dann die kom- 
plexen Amplituden y, ‚langsam veränderlich‘ sind. Mittelt man also über eine 


Periode der h-ten Eigenschwingung, so folgt mit p,, = it ; 
dy 1 
J I(t) dt— en yy, (45) 
Por = T, 


Aus (45) erkennt man das Abklingen ursprünglich angefachter höherer Eigentöne, 
wenn der Grundton die Oberhand gewonnen hat und praktisch für sich allein den 
Verlauf des Anodenstromes /(t) bestimmt. Während der Halbperiode der Grund- 
tonschwingung, in der nach Abb. 3 der Luftstrom in das Pfeifeninnere eintritt, 
ist nämlich J während dieser Zeit als konstant zu betrachten. Das Integral in (45) 
ist dann für höhere Eigenfrequenzen gleich Null. Aus (45) geht also in der Tat 
hervor, daß die Amplitude y, nach Maßgabe ihrer Dämpfungskonstanten ver- 
klingt: 
yv„=Be®t, 


Die Anfachung der Grundschwingung in nichtlinearen Fall kann auf Grund 


des Energieprinzips in folgender Weise behandelt werden. Auf Grund von Abb. 6 
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liefert der Energiesatz für die als Fre veränderlich zu betrachtende Energie E 


= VI-— (46) 

Darf man (bei schwacher Riickkopplung) fiir V einen praktisch sinusférmigen 

Verlauf voraussetzen, so folgt durch Mittelung von (46) über eine Periode, wenn 
wir setzen: V = Vysinwt 

dE 


anf I(t) sin w tdw t— 


(47) 


Die Schwingungsenergie E ist dem Quadrat der Spannungsamplirude V¢ propor- 
tional. Wir setzen ®) 


12 
und erhalten damit aus (47): 
dV, 
"19 sin tdo Vy, (48) 
Nach (28a) ist im Integral zu setzen 
I(t) (0 V, sin ot). (48a) 


Fiir die Anwendung von (48) betrachten wir zwei besonders einfache Beispiele fiir 
die Ermittlung der Anfachung der Spannungsamplitude V,. 

1. Die Charakteristik G (© V) sei geradlinig mit der Steilheit S ohne Sättigung, 
der Arbeitspunkt liege im unteren Knick (vgl. Abb. 7). Dann wird 


27 


J G(O V, sin x) sinadz = OSV, = 3.058 


3) Unter Zugrundelegung von Abb. 6, in der wir in erster Näherung die Dämpfung 
unberiicksichtigt lassen wollen und auch von der Existenz des Kondensators C absehen 
können, ergibt sich eine sinusférmige Abhängigkeit der Leitungsspannungsamplitude 
vom Ort: 

V=Asinkz. (I) 
An der Induktivitat Z liegt die Spannung V5, so daß (I) übergeht in 


Vo 
= oh sink2. 


Bezeichnen wir die Kapazität pro Längeneinheit der Leitung mit K und bedenken, daß 
die Schwingungsenergie E bei verschwindendem Leitungsstrom J rein elektrostatisch 
gegeben ist, so folgt 


1 ı 
E= [Md= 
0 
Zieht man (41b) heran und beachtet, daß der Wellenwiderstand der Leitung durch 


1 
2 


cZ = 
gegeben ist, so folgt durch Berechnung des Integrales mit dem Wert (41a) für 8 
Vi 
E= 18° 


Im Fall des einfachen Schwingkreises ist HR: fB-! = 2 C, was auch direkt aus (8) ent- 
nommen werden kann. 
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also nach (48) 


oder integriert 
a= (F-4). (49) 


Der Vergleich mit (36) zeigt, daß die ,,mittlere Steilheit‘‘ in unserem Beispiel gerade 
halb so groß ist wie im linearen Fall (36). Nach (49) ergibt sich eine Mindest- 
steilheit, für die eine Entdämpfung noch möglich ist. Im linearen Fall ist diese 
Steilheit durch (44a) gegeben. In unserem nicht linearen Beispiel (1) wird 


2 
Smin = Su’ (50) 


eV 


Abb. 7. Charakteristik 1 . Abb. 8. Charakteristik 2 


2. Die Charakteristik sei unendlich steil und zeige Sättigung fiir 


vgl. Abb. 8. 
Hier wird 
2n 


f Isin =I, f sinxdx = 2],. 
6 6 


Es folgt also nach (48) 
28, 
8 


a EV (51) 


wo nach (33) £ die Dämpfungskonstante e des Schwingkreises bedeutet. 

Aus (51) geht hervor, daß sich in unserem Fall des Beispiels (2) die Spannungs- 
amplitude V, exponentiell einer stabilen Endamplitude V(oo) annähert, für die Me 
verschwindet : 


V(co) = =, w. (52) 


Der exponentielle Anstieg ist nach (51) durch die Dämpfungskonstante e = B 
des Schwingsystemes bestimmt. 
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Das Ergebnis (52) ist verständlich, wenn wir bedenken, daß die Fouriersche 
Grundschwingung des Anodenstromes im Falle der Charakteristik (Abb. 8) durch 


21 , gegeben ist und daß der Wechselstromwiderstand des mit der Resonanz- 


frequenz arbeitenden Schwingkreises unendlich groß ist, so daß dann nur die 
Paralleldämpfung w übrigbleibt. 


Unser Beispiel (2) erscheint geeignet, als Modell für die Vorgänge in der 
Lippenpfeife (Abb. 3) zu dienen. 


Für eine akustische Diskussion müssen (42) und (43) in akustische Gleichungen 
transformiert werden. Im Falle der Orgelpfeifen (Abb. 3) ist die ‚„Induktivität“ 
durch die Menge der am Lippenspalt mitschwingenden Luft gegeben. Der Druck p 
im Innern der Pfeife entspricht elektrisch der Leitungsspannung, der elektrischen 
Stromstärke ordnen wir das Produkt Fv zu, wo F der jeweilige Pfeifenquerschnitt 
und v die Geschwindigkeit der Luftteilchen bedeutet. Es ist üblich, Druck und 
Geschwindigkeit durch das „Geschwindigkeitspotential‘“ qm auszudrücken. Be- 
deutet o die Luftdichte, so ist bekanntlich 


v = — grad 


Das Potential m gehorcht der Wellengleichung 


Im Pfeifenrohr gelten die Ansätze (mit b= 2) 


y=Asinkxsinwt a) 
=—Akcoskzsinwt b) (53) 
p=Aowsinkxcoswt c) 


x = 0 bedeutet hier das offene Ende der Pfeife, wo praktisch p = 0 gesetzt werden 
kann. Als Randbedingung an der Lippe fordern wir in Analogie zum Ohmschen 
Gesetz 

g=W-t, (54) 
wo 

I=v,:f=v-F 


der Schallstrom durch den Lippenspalt mit dem Querschnitt f bzw. durch den 
unteren Teil der Pfeife mit dem Querschnitt F bedeutet. 


Das Geschwindigkeitspotential m setzen wir im AuBenraum gleich Null. W be- 
deutet einen Strömungswiderstand, der durch die geometrische Form der Lippe 
bestimmt ist. Wie an anderer Stelle näher ausgeführt werden soll, wird man für 
rechteckige Spalte quadratischer Pfeifen mit der Quadratseite d nicht allzu fehl 
gehen, wenn man setzt 


W=- (55) 


Ay 
1) 
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Dabei bedeutet H die Höhe des Lippenspalts (,,Aufschnitt“). Setzt man F = d, 
I=vF, so folgt aus (53a) und (53b) durch Division 


(56) 


Ein Vergleich mit (41 b) liefert dann mit der früher benutzten Abkürzung + = 2» 


L 


2 NZ (57) 


list die Länge des Pfeifenrohres. Wie man sieht, kommt es in erster Linie auf das 
„Mensurverhältnis‘ $ der Pfeife an. 


Da wenigstens die ersten Eigenténe nach (56) noch quasi harmonisch liegen, 
kann, wenn h die Ordnungszahl des Eigentons angibt, in (42) annähernd ersetzt 
werden: 


= h?- 72 (58) 


Wählen wir der Einfachheit halber für eine bestimmte Pfeifenkonstruktion 


d 
H™ 1,5? = 100, (59) 
so wird 
Mon L? ad 2 vL_ d 
100h? , 3 T° (60) 


(42) geht damit über in 
l a 
1+ 54 (1+ 100m 5 


(61) 


Für den Grundton ist A = 1 zu setzen. Eigentöne oberhalb einer bestimmten 
Schranke für h, für die der Zähler Null oder negativ wird, können überhaupt 
nicht anspringen. S muß als Funktion der Zeit betrachtet werden, wie oben aus- 
einandergesetzt wurde. 

Die Amplituden, mit denen die Eigentöne anspringen, können einfach nur 
unter Zugrundelegung einer linearen Charakteristik ermittelt werden. Wir geben 
deshalb in folgendem noch eine strenge Lösung für den linearen Fall. 

Strom und Spannung seien also nach (15) miteinander verkoppelt und man wird 
auf Grund von (24) erwarten, daß eine Partikularlösung von (6) gegeben ist durch 


V=4AePt, (62) 
wobei wir voraussetzen können, daß der Realteil von P positiv ist: 
RP>0. (63) 


Setzt man in (6) 7@ = p und trägt den Ansatz (62) ein, so liefert die Durchführung 
der Integration über ¢’ 


R(p) er: dp = 428 (64) 
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Zur Durchführung der Integration beachten wir, daß die Pole p, des Zählers 
nach Abb. 5 alle links der imaginären Achse liegen, die Nullstelle des Nenners 
jedoch wegen der Voraussetzung (63) auf der anderen Seite liegt, vgl. Abb. 9. 
Der Integrationsweg kann also um den Pol P 
zusammengezogen werden und liefert 


V=AOSeP'*R(P). 
Damit unser Ansatz (62) eine Lésung von (6) ist, p 
muß also P der Bedingung genügen Mn® 
1 


was die Parallelschaltung des ‚negativen Wider- 
standes“ (O8)! zu R(w) bedeutet. Die allgemeine 
Lösung von (6) werden wir in Übereinstimmung 
mit (24) in der Form ansetzen p-Ebene 


V(t) = SA, Pot (66) Abb. 9 


(65) ist die strenge Bedingungsgleichung fiir die komplexen Frequenzen P, fiir die 
bereits nach (36) eine Näherungslösung vorliegt. Es ist leicht zu sehen, daß diese 
Näherungslösung auch aus (65) erhalten werden kann. Setzt man nämlich an Stelle 
von (65) analog wie bei (32) 


o(P,) + OS 
(Pon) 0, 
so folgt durch Subtraktion, fiir den Fall, daB OS —+ hinreichend klein ist, 
h 


(P,— Po) (GF), = 98 


Wr 


Also ergibt sich unter Beriicksichtigung von (31) und (32) aus (29) 
1 
Pi — Por = (OS——)Br 


in Übereinstimmung mit (36). Die strenge Lösung (65) liefert gleichzeitig eine 
Möglichkeit zur Beurteilung der Genauigkeit dieser Näherungslösung. 

Die Amplituden A, in (66) können nur auf Grund einer besonderen Annahme 
über die Ursachen der Anfachung ermittelt werden. Wir wollen annehmen, der 
ganze Schwingungsvorgang möge dadurch in Bewegung kommen, daß in die 
Rückkopplungsspule R der Abb. 1 z. Z. t’ = 0 eine scharfe Spannungsspitze Eg(t') 
hinein induziert wird, wodurch eine bestimmte Elektrizitätsmenge Q die Induk- 
tivität Z durchsetzt und dabei alle Eigenschwingungen anstößt. Unter Hinzu- 
nahme dieser Spannungsspitze, die zusammen mit der rückgekoppelten Spannung 
© V den Anodenstrom /(t’) bestimmt, geht dann die Integralgleichung (1) über in 


t 
(67) 


= a, 

(56) 
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Die Funktion f(t — ¢’) ist durch (7) gegeben. Findet der Spannungsstoß z. Z. 
t' = 0 statt und war vorher alles in Ruhe, so liefert (6) 


t t 
ViQ=S f E,W’) (68) 
0 -00 


Wir wollen der Einfachheit halber annehmen, daß der Spannungsstoß L(t’) 
nur so kurz andauert, daß während seiner Dauer /(t—t’) durch f(t) ersetzt werden 
kann. Dann folgt einfach, weil die durchgegangene Elektrizitätsmenge durch 
Tt 
Q=S8 f at 
6 


gegeben ist: 
t 
Vi)=S8 ft — dt’ —Q fit). (69) 


Die Funktion f(t — t’) ist durch (9) gegeben. 


Tragen wir (66) und (9) in (69) ein, so folgt nach Durchführung der Inte- 
gration über (¢’) 


Die Forderung (70) erfüllen wir dadurch, daß wir die P, und A, den Bedingungen 
unterwerfen: 


1 

O8 
Q A, _ 

as— (72) 


Die Pole p, sind die Wurzeln von 
H-(p) = 0, 
die 6, sind durch (41a) gegeben. 


Aus (71) können wir der Reihe nach die P, berechnen. (72) stellt für = 1, 2,3... 
ein lineares System zur Berechnung der Amplituden A, dar, womit unsere Auf- 
gabe im linearen Fall streng gelöst ist. 


Zunächst ist leicht zu sehen, daß (71) mit (65) identisch ist. Bedenkt man 
nämlich, daß alle p, links der imaginären Achse und die P rechts davon liegen, 
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so kann das Integral 
+400 
/ R(p) dp 
i 


entweder durch Umgang um die Singularität P, oder durch die Umgänge um die 
Pole p, berechnet werden. Deshalb ist nach (8) 


= 


womit in der Tat (71) in (65) übergeht. (36) besteht also zu Recht und kann auch 


aus (71) entnommen werden, wenn alle Glieder der Summe außer dem größten, 
für den der Nenner besonders klein wird, vernachlässigt werden. Dann folgt nämlich 


P,— = 


was nach (25) und (33) mit (36) identisch ist. 


Mit demselben Genauigkeitsgrad folgen dann aus (71) und (72) durch Division 
die Amplituden A, 


A, = (73) 
Unter Benutzung von (41) ergibt sich also 


Q L 
1 + 


(74) 


Die Lösung (43) ist damit restlos bekannt. x, ist durch (42) gegeben, im akusti- 
schen Fall durch (61). 


Für die Beurteilung der Hörbarkeit der Eigenfrequenzen muß berücksichtigt 
werden, daß der Energiestrom E durch den Lippenspalt der Pfeife gegeben ist 
durch 


E= (75) 
fv, ist die Amplitude des Stromes, die aus der Spannungsamplitude (74) durch 


Division mit w,, Z hervorgeht. Der abgestrahlte Energiestrom durch den Lippen- 
spalt ist also im Anschluß an (61) auf Grund von (74) 


h* vi eth t 


=K 
l1+ + 100» 
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Aus (76) ist zu entnehmen, daß bei der Anfachung höhere nicht harmonische 
Eigentöne bis zu einer gewissen Frequenzgrenze anlaufen, für die nach (61) « „Vver- 
schwindet. Wie früher auseinandergesetzt wurde, klingen jedoch diese Eigentöne 
der Reihe nach wieder ab, wenn die Grundschwingung sich ihrer stabilen 
Endamplitude nähert — in qualitativer Übereinstimmung mit den Tren- 
delenburgschen Experimenten. 


Frankfurt a.M., Institut für angewandte Physik der Universität. 


(Bei der Redaktion eingegangen am 9. April 1952.) 
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